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PREFACE 


The  information  documented  in  the  publication  is  meant  to  provide  the  ffee-flight 
ballistics  test  engineer  with  the  basic  insights  and  understandings  he  requires  to 
accomplish  his  research  duties.  However,  it  certainly  does  not  attempt  to  provide  all  the 
information  the  ballistic  range  practitioner  needs  to  fulfill  these  duties.  There  are 
numerous  outstanding  references  available  to  the  researcher.  Some  of  these  are  listed 
herein  and  those  references  will  refer  to  even  more  outstanding  publications  associated 
with  ballistic  technologies.  It  is  hoped  that  the  information  contained  herein  will  be  of 
assistance  to  the  ballistic  range  test  personnel,  especially  those  involved  in  free-flight 
spark  ranges. 

This  publication  was  drafted  intermittently  from  December  1990  through  March  1997. 
Information  was  extracted  from  numerous  papers,  reports,  and  text  books.  Hopefully, 
credit  has  been  given  where  due  but  after  practicing  over  30  years  in  a  technology  area  it 
is  difficult  to  remember  where  all  the  information  one  possesses  originally  came  from. 
Nevertheless  applicable  references  have  been  provided  throughout  this  publication  and  a 
bibliography  of  related  publications  is  also  provided  for  the  interested  reader.  Because 
much  of  the  information  was  extracted  from  other  sources  and  collected  over  a  period  of 
time  the  units  used  throughout  are  a  hodge  podge  of  English  and  Metric  and  in  many 
cases  mixed.  The  author  apologizes  for  this  inconvenience  but  has  attempted  to  at  least 
clearly  identify  the  units  used. 

Even  after  the  above  statement,  it  is  still  necessary  to  acknowledge  certain  individuals, 
i.e.,  Bob  Whyte  and  Wayne  Hathaway  of  Arrow  Tech,  Inc.,  Burlington  VT,  Dr.  C.  H. 
Murphy,  Army  Research  Laboratory,  Aberdeen  MD,  Dr.  John  D.  Nicolaides,  California 
Polytechnical  Inst,  (retired),  and  Dr.  Gary  Chapman,  University  of  California,  Berkeley. 
Although  all  of  these  individuals  have  been  referenced  they  are  the  persons  primarily 
responsible  for  much  of  the  information  provided  herein.  Very  little,  if  any,  of  the 
information  presented  in  this  publication  was  original  work  by  the  author.  Some  of  the 
derivations  included  herein  may  appear  to  differ  from  those  in  previous  publications; 
however,  this  is  only  due  to  the  author’s  style. 

There  is  always  the  desire  and  inclination  to  add  one  more  section  or  present  another 
derivation.  Also,  since  modem  computational  techniques  and  electronics  are  rapidly 
evolving  some  sections  have  required  significant  rewrites  during  the  drafting  process.  In 
fact  this  could  be  a  continuous  process.  But,  eventually  an  end  must  be  reached  and  the 
work  product  delivered  in  the  hope  that  it  adds  something  to  the  available  literature.  In 
that  hope  this  document  is  submitted. 
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L  INTRODUCTION 


1.1  HISTORICAL  BACKGROUND 

Ballistics  could  very  well  be  the  oldest  technology  known  to  man.  It  certainly  predates 
fire  making  and  may  even  predate  the  use  of  primitive  tools.  Undoubtedly  one  of  man’s 
two  legged  ancestors,  in  a  fit  of  anger,  picked  up  a  stone  and  threw  it  at  a  perceived 
enemy  or  rival.  After  this  rival  figured  out  that  he  could  do  likewise,  an  enterprising 
individual  must  have  recognized  that  whoever  could  throw  the  stone  the  farthest  had  a 
distinct  advantage.  Especially  if  the  stone  could  be  thrown  with  any  accuracy. 

This  recognition  would  have  led  to  the  development  of  crude  slings  which  can  still  be 
found  in  use  by  primitive  tribes  even  today.  Although  we’ll  never  know  for  certain 
whether  or  not  this  scenario  is  correct  we  do  know  that  strife  and  war  have  been  the 
engine  driving  the  development  of  ballistics  as  a  science  during  the  following  ages.  These 
crude  slings  evolved  to  become  bows  and  arrows  which  were  used  as  personal  weapons 
and  even  later  as  catapults  which  were  the  medieval  equivalent  of  artillery. 

However,  the  engine  really  didn’t  get  into  high  gear  until  the  Chinese  invention  of  gun 
powder  was  brought  to  Europe  by  Marco  Polo  in  the  1 3th  century.  This  invention  not 
only  made  it  possible  to  project  force  even  farther  using  guns  or  rockets  but  also  the 
projectile  could  carry  an  explosive  charge  and  hence  deliver  more  energy  to  the  target 
than  just  the  kinetic  energy  of  the  projectile  itself.  This  single  event  totally  altered 
western  civilization  as  well  as  warfare  because  castles  suddenly  became  very  vulnerable 
to  attack.  Since  the  sole  purpose  of  these  damp  musty  piles  of  stone  was  to  provide 
protection  to  the  nobility  and  ruling  classes  and  this  protection  was  no  longer  available 
the  residents  of  these  castles  moved  into  the  more  comfortable  towns  and  villages.  This 
intermingling  of  the  ruling  classes  with  the  masses  of  the  people  was  the  beginning  of  the 
breakdown  of  the  existing  social  structure.  From  this  breakdown  followed  the 
establishment  of  a  new  social  order  which  has  progressed  to  the  modem  concepts  of 
personal  and  religious  freedoms  that  we  enjoy  today. 

Nevertheless  the  original  concept  that  whoever  could  throw  the  stone  the  farthest  had 
an  advantage  over  his  enemies  was  still  valid.  Napoleon  used  this  concept  by  first  using 
sabots  in  firing  subscale  projectiles  from  large  bore  gims.  These  subscale  projectiles 
could  be  fired  much  farther  since  the  total  in  gim  weight  of  the  projectile  plus  sabot  was 
significantly  less  than  the  weight  of  the  full  scale  projectile.  This  permitted  Napoleon’s 
gunners  to  engage  the  enemy  troops  long  before  they  could  return  fire.  Although  the  use 
of  sabots  added  a  complexity  to  the  accuracy  of  which  the  projectiles  could  be  fired, 
Napoleon’s  artillery  tactics  became  legendary  and  accounted  for  much  of  his  military 
successes. 
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Very  early  in  the  development  of  ballistics  as  a  science,  ballisticians  realized  that  if  a 
projectile  could  be  fired  at  a  high  enough  velocity  then  the  curvature  of  the  trajectory  due 
to  gravity  would  equal  the  curvature  of  the  Earth.  Since  the  projectile  would  then 
essentially  fall  around  the  Earth,  orbit  would  be  achieved.  This  is  the  basic  concept  from 
which  orbital  mechanics  was  derived.  Obviously,  once  we  had  developed  the  capability  to 
fire  a  projectile  at  this  orbital  velocity  we  had  reached  the  ultimate  ability  to  throw  our 
stone  to  the  farthest  possible  distance. 

The  purpose  of  this  historical  diatribe  is  not  to  provide  a  history  lesson  of  the  world 
but  to  give  the  reader  a  basic  understanding  of  the  importance  that  ballistics  has  played  in 
the  development  of  civilization  and  science  as  we  know  it  today. ' 

1.2  EXTERIOR  BALLISTICS 

Since  we  have  now  reached  the  ultimate  of  throwing  our  stone  the  farthest  possible 
distance,  future  advances  in  ballistics  will  be  primarily  directed  in  three  general  areas. 
First  to  improve  the  efficiency  of  the  launch  systems  in  order  to  optimize  the  velocity  and 
mass  of  the  payload  that  we  want  to  deliver.  Secondly,  to  improve  our  ability  to  hit  the 
target  of  interest.  And  lastly,  to  deliver  more  energy  and  maximize  damage  to  the  target. 
These  three  areas  are  the  distinct  technologies  commonly  called  interior,  exterior,  and 
terminal  ballistics. 

Interior  ballistics  is  the  study  of  the  processes  which  occur  while  the  projectile  is  in 
the  barrel  of  the  gun.  Exterior  ballistics  is  concerned  with  the  understanding  of  the  forces 
and  moments  acting  on  the  projectile  during  its  flight.  Terminal  ballistics  is  the  science  of 
what  happens  when  the  projectile  hits  the  target.  Although  all  three  of  these  areas  will  be 
mentioned  in  varying  degrees  herein,  the  primary  discussion  will  be  assoeiated  with  the 
exterior  ballistics  of  the  projectile  or  missile. 

The  historical  event  from  which  exterior  ballistics  became  a  true  science  occurred  in 
the  mid  1 8th  century.  At  that  time  Benjamin  Robbins  used  a  ballistic  pendulum  to 
measure  the  muzzle  velocity  of  a  projectile.  Until  then  muzzle  velocities  were  estimated 
by  measuring  the  distance  the  projectile  flew  and  calculating  the  velocity  using  the  point 
mass  kinematic  relations  derived  by  Newton  and  others  a  hundred  years  earlier.  Since 
these  relations  ignored  the  air  density  they  grossly  imderestimated  the  associated  muzzle 
velocities.  This  apparent  discrepancy  between  the  measured  muzzle  velocity  as 
determined  by  Robbins  and  the  best  theoretical  estimates  of  the  time  led  to  the 
understanding  that  the  air  itself  provided  a  resistance  force  to  the  passage  of  the 
projectile.  The  concept  of  aerodynamic  drag  was  henceforth  developed  and  a  realization 
that  this  force  was  a  function  of  the  shape  of  the  projectile  formed  the  basis  for  the 
science  of  aerodynamics. 
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Since  Robbins  made  his  muzzle  velocity  measurements  many  individuals  have 
contributed  to  the  science  of  exterior  ballistics.  Much  of  the  more  recent  work 
accomplished  by  some  of  these  individuals  will  be  described  and  referenced  herein. 

1.3  FREE  -  FLIGHT  SPARK  RANGES 

A  free-flight  spark  range  is  an  enclosed  facility  containing  a  number  of  orthogonal 
spark  photography  stations  and  assorted  other  instrumentation  systems.  Various 
aerodynamic  configurations  of  interest  are  flown  through  the  facility  and  shadowgrams 
are  obtained  at  each  of  the  spark  stations.  These  shadowgrams  are  then  processed  and 
analyzed  to  provide  an  experimentally  measured  position-attitude-time  history  of  the  free- 
flight  configuration.  The  theoretical  equations  of  motion  are  then  fit  to  the  experimentally 
measured  trajectory  and  the  unknown  aerodynamic  coefficients  and  derivatives  extracted. 
These  aerodynamic  parameters  are  the  basic  data  required  to  understand,  predict,  or 
improve  the  exterior  ballistics  performance  of  the  configuration  of  interest.  This  complete 
process  will  be  described  and  discussed  in  some  detail  in  the  following  sections. 

Many  varieties  of  test  ranges  must  have  been  used  by  researchers  over  the  centuries. 
But,  the  first  spark  photograph  of  a  free-flight  object  was  obtained  by  Ernst  Mach  late  in 
the  19th  century.  This  development  was  put  to  good  use  by  A.C.  Charters  when  he 
constructed  and  operated  the  first  modem  free-flight  spark  range  ’  ^  at  Aberdeen 

Proving  Ground,  Maryland,  during  World  War  II.  That  facility  proved  very  useful  in 
solving  many  ballistic  problems  and  provided  the  impetus  in  developing  theoretical  free- 
flight  dynamic  analysis  techniques.  That  facility  is  still  a  very  valuable  research  tool  and 
the  technologies  developed  there  have  been  used  throughout  the  world. 

Since  World  War  II  ballistic  ranges  have  been  adapted  to  study  many  different 
physical  phenomena,  i.e.,  impact  mechanics,  basic  material  properties,  reentry  physics, 
ablative  properties  of  various  materials,  wake  studies,  radar  cross-section  measurements, 
planetary  entry  probes,  and  many  others.  As  a  result  ballistic  range  technology  has 
become  a  diverse  and  multifaceted  mature  science'  This  is  not  to  suggest  that  we  know 
everything  there  is  know  about  this  technology  area,  but  only  that  the  beginner  has  a  lot 
to  know  about  before  he  can  become  expert  in  the  field. 

The  primary  discussion  herein  will  be  related  to  free-flight  spark  ranges  as  applied  to 
aerodynamic  testing  which  was  Charters  primary  purpose  in  constmcting  that  first 
facility.  The  facilities  that  have  been  constmcted  since  have  varied  considerably  in  size, 
number  of  spark  photography  stations,  geometry,  and  capabilities.  Hence,  the  facility 
used  herein  to  describe  the  various  processes  and  considerations  associated  with  free- 
flight  aerodynamic  testing  is  the  Aeroballistic  Research  Facility  (ARF)'  ''  located  at  Eglin 
AFB,  Florida  and  operated  by  the  Wright  Laboratory  Armament  Directorate. 

This  atmospheric  test/research  facility,  see  Fig.  1.1,  became  operational  in  1976.  It 
contains  55  orthogonal  shadowgraph  stations,  see  Fig.  1.2,  along  with  other 
instrumentation  systems,  i.e.,  interferometry,  Schlieren,  front  lighted  LASER 
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photography,  atmospheric  measuring  instruments,  and  a  master  chronograph  system. 

This  enclosed  facility  is  about  800  ft  long  with  660  ft  ( ~  200  meters)  of  instrumented 
length  and  has  been  used  to  evaluate  the  ffee-flight  characteristics  of  various  objects  from 
subsonic  to  hypersonic  velocities.  A  list  and  locations  of  the  instrumentation  systems 
presently  contained  in  the  ARF  are  presented  in  Appendix  1 . 


Figure  1 . 1  Sketch  of  the  Aeroballistic  Research  Facility 
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Figure  1.2  Schematic  of  a  Shadowgraph  Station 
1.4  AERODYNAMIC  TESTING 

Webster’s  definition  of  aerodynamics  is  “  the  dynamics  of  gases,  especially  of 
atmospheric  interactions  with  moving  objects  Therefore,  aerodynamic  testing  is  the 
measurement  /  determination  of  these  atmospheric  interactions.  The  basic  tenet  in  wind 
tunnel  testing  is  that  instead  of  moving  the  object  (model)  the  object  is  held  and  the 
atmosphere  (air  /  gas)  is  moved  over  the  model.  The  idea  being  that  the  interactions 
occurring  between  the  atmosphere  and  model  is  independant  of  whether  or  not  the  model 
or  the  gas  is  moving.  Although  this  is  generally  true  there  are  many  issues  concerning  the 
quality  of  flow  in  wind  tunnels  which  can  affect  the  resulting  aerodynamic  data.  Some  of 
these  issues  are  the  tunnel  wall  effects,  variation  of  the  velocity  across  the  tunnel  test 
section,  flow  turbulence,  interactions  caused  by  the  moimting  system,  i.e.,  model  base 
flow  interference  effects,  condensation  of  the  flow  and  others.  Many  of  these  aspects 
associated  with  wind  tunnel  testing  are  discussed  in  various  publications,  see  Pope’s  book 
entitled  “Low-speed  Wind  Tunnel  Testing”  * 

In  contrast  to  wind  tunnel  testing,  ballistic  range  testing  is  the  study  of  the  interaction 
of  ‘  moving  ’  objects  in  a  quiet  or  still  atmosphere.  Hence  the  concerns  associated  with 
the  quailty  of  flow  are  eliminated.  Nevertheless,  in  whichever  type  of  facility  is  used 
there  are  flow  properties  which  must  be  matched  in  order  to  duplicate  the  full  scale  flight 
conditions  or  to  compare  data  from  facility  to  facility.  Two  of  the  most  important  of  these 
are  Mach  number  ,(M),  and  Reynolds  number  ,(ReL).  These  are  defined  as: 
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M  =  V/a 

1.1 

and 

ReL  =  pVL  /  p 

1.2 

where : 

V  =  velocity  of  the  model ,  or  velocity  of  the  flow  in  the  test  section 

1  n 

a  =  speed  of  sound  (local),  (  yRgT) 
p  =  density  of  flow  medium 
L  model  length 
p  =  viscosity  of  flow  medium 

Since  air  viscosity  is  insensitive  to  pressure  a  linear  relationship  for  p  exists  for  the 
60-90  degree  Fahrenheit  temperature  range.  This  relationship  for  p  of  air  is , 

p  =  3.7463x1 0'”^  +  (0.00551 5x1 0'^)(Tp  -  60)  1.3 

where  p  has  the  units  of  Ib.sec.  /  ft.  and  Tp  is  in  degrees  F.  Flow  Reynolds  numbers  are 
frequently  presented  as  per  unit  length,  or,  Re  =  pV  /  p ,  where  the  model  length  is  not 
included.  There  are  many  other  numbers  which  are  used  to  identify  and  qualify  various 
flow  regimes,  i.e.,  Knudsen  number  ~  M  /  Re,  but  it  is  not  our  purpose  here  to  discuss  alt 
of  these  various  flow  quantification  parameters. 

The  test  engineer  should  also  recognize  that  when  testing  subscale  objects,  parameters 
other  than  flow  conditions  should  also  be  matched  with  the  full  scale  object  if  possible. 
The  center-of-gravity ,  L^g  /  L  ,  and  the  inertia  ratio  ,  /  ly ,  are  two  such  model 

characteristics  that  should  be  matched  when  feasible.  Although  the  aerodynamic  moment 
data  can  be  shifted  to  any  Leg  /  L  position  this  shift,  as  will  be  shown  in  Section  VII, 
requires  a  determination  of  the  normal  force  coefficient  and  if  this  determination  is  in 
error  the  shifted  moment  data  will  also  possess  errors.  Therefore  it  is  wise  to  design  the 
subscale  models  with  L^g  /  L  positions  as  close  to  the  full  scale  vehicle  as  practical  and 
minimize  the  required  shift  in  the  moment  data.  Also  since  the  inertia  ratio  can  affect  the 
free-flight  angular  motion  pattern,  this  ratio  should  be  representative  of  the  full  scale 
vehicle. 

Here  a  precautionary  note  concerning  subscale  testing  in  any  facility  must  be 
mentioned.  Boundary  layers  do  not  scale!  Therefore,  the  boundary  layer  thickness  will  be 
overly  large  compared  to  the  full  scale  configuration.  This  overly  large  boundary  layer 
can  affect  the  effectiveness  of  the  control  surfaces  (i.e.  fins)  immersed  in  the  boundary 
layer  and  result  in  misleading  results  and  conclusions.  Also  the  boundary  layer 
transition  point  ( laminar  to  turbulent  flow)  will  not  scale.  Frequently  an  experimentalist 
will  force  the  boundary  layer  to  transition,  using  trips,  at  the  equivalent  location  of  the 
full  scale  object  in  order  to  better  simulate  the  actual  flow  conditions.  Tripping  boundary 
layers  can  be  an  art  in  itself  and  beads,  grooves,  grit,  and  other  mechanisms  are  used. 
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When  comparing  dynamic  data  (free-flight,  air  bearing,  forced  oscillation,  etc.)  two 
additional  parameters  can  be  important.  These  are  the  reduced  frequency  parameter, 
tod  /  2V,  and  the  nondimensionalized  spin  parameter,  pd  /  2V.  Omega,  to,  is  the  frequency 
of  oscillation  and  in  the  free-flight  facility  is  normally  assumed  to  be  equal  to  the 
nutational  frequency,  (1)n'  ,  and  p  is  spin  rate.  The  characteristic  length  (d)  is  usually  the 
model  diameter  but  when  testing  aircraft  type  configurations  is  frequently  defined  as  the 
aerodynamic  chord  of  a  wing  or  fin.  These  parameters  should  be  recorded,  and  the 
determined  aerodynamic  results  plotted  versus  cod  /  2V  and  pd  /  2V  to  determine  if  they 
are  affecting  the  results.  Examples  of  cases  where  they  have  affected  the  results  are 
discussed  in  Refs.  1.6  and  1.7. 

The  ultimate  purpose  of  aerodynamic  testing  is  to  measure  those  interactions  between 
the  test  model  and  the  atmosphere  in  such  a  way  as  to  quantify  those  interactions  .  This 
quantification  is  accomplished  by  measuring  the  forces  and  moments  acting  on  the  test 
object  as  a  result  of  these  interactions.  Once  these  forces  and  moments  have  been 
measured  they  can  be  used  to  predict  the  performance  of  the  object  at  conditions  other 
than  those  tested.  Hence  the  vehicle’s  flight  performance  is  known  and  a  flight  control 
system  can  be  designed.  Obviously  for  those  configurations  which  do  not  have  a  flight 
control  system,  i.e.,  dumb  bombs,  purely  ballistic  missiles,  bullets,  artillery  shells,  etc., 
the  accurate  determination  of  these  forces  and  moments  are  an  absolute  necessity  in  order 
to  predict  the  trajectory  from  the  platform  to  the  target. 
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II.  DEFINITIONS 


2.1  COEFFICIENT  AND  DERIVATIVE  NOMENCLATURE 

The  purpose  of  the  aerodynamic  nondimensionalized  nomenclature  is  to  represent  the 
aerodynamic  forces  and  moments  acting  on  a  body  in  a  manner  such  that  these  forces  and 
moments  can  be  easily  used  in  conjunction  with  the  equations  of  motion  to  predict  the 
flight  performance  of  the  body  in  question.  This  nomenclature  permits  measuring  the 
aerodynamic  parameters  using  a  subscale  body,  i.e.,  from  either  a  wind  tunnel  or  free- 
flight  spark  range,  and  using  these  measurements  directly  in  the  predictions.  It  also 
permits  direct  comparison  of  the  aerodynamic  forces  and  moments  as  measured  from 
various  scale  models  in  different  facilities. 

In  general  forces  are  represented  by  capital  letters,  i.e.  F,  and  moments  are  represented 
by  lower  case  letters,  i.e.  m.  The  nondimensionalized  force  coefficient  is  defined  as, 

Cf  =  F/(1/2)pV^A  2.1 


and  the  moment  coefficient  as. 


C„  =  m/(l/2)pV^Ad  2.2 

Where  A  and  d  are  the  model’s  reference  area  and  reference  length  respectively.  The 
reference  area  and  length  should  always  be  clearly  specified  when  presenting  or 
documenting  results.  The  aerodynamic  coefficients  are  useless  in  comparing  results  or 
predicting  the  flight  performance  if  these  reference  values  are  not  specified.  It  is  also 
imperative  to  specify  the  moment  reference  location  as  well.  This  is  normally  the  center- 
of-gravity  location  for  the  free-flight  models;  however,  not  always. 

The  forces  in  the  x,y,z  directions  are  therefore  defined  as: 

F,-C,(l/2)pV"A  1 

F,-C,(l/2)pV"A  }  2.3 

F,.C.(l/2)pV"A  J 

And  the  moments  about  the  x,y,z  axis  respectively  are  defined  as: 

1  =  C,  (1/2)  p  V^  A  d  ] 

m  =  C^(l/2)pV^Ad  \  2.4 

n  =  Cn(l/2)pV^Ad  J 
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Since  the  x  axis  is  associated  with  the  direction  of  flight  (see  Section  2.2.2)  1  is  the  rolling 
moment,  m  is  the  pitching  moment,  and  n  is  the  yawing  moment.  The  derivative 
nomenclature  of  the  force  and  moment  coefficients  is  illustrated  by  the  examples  below: 


5Cx  /  5M  =  CxM 

2.5 

dC^  /  da  = 

2.6 

where  M  is  Mach  number  and  a  is  the  angle  of  attack.  Hence,  C^m  in  Eq.  2.5  represents 
the  slope  of  the  force  coefficient  in  the  x  direction  with  respect  to  Mach  number. 

Likewise  ,  in  Eq.  2.6  represents  the  slope  of  the  pitching  moment  with  respect  to  the 
angle  of  attacL  This  moment  coefficient  derivative  ,  is  frequently  referred  to  as  the 
static  stability  derivative  and  will  be  discussed  in  detail  in  the  later  sections. 

The  aerodynamic  nomenclature  as  defined  above  is  used  throughout  the  following 
sections  and  is  the  basic  nomenclature  system  used  by  most  designers  and  researchers 
around  the  world.  However,  in  the  past  some  organizations  /  countries  did  not  include  the 
1/2  in  the  force  and  moment  definitions  of  Eqs.  2.1  and  2.2.  Therefore,  when  comparing 
data  or  using  the  results  from  other  organizations  one  should  be  alert  to  how  the  force  and 
moment  coefficients  are  defined. 

2.2  AXES  SYSTEMS 

2.2.1  Facility  Axis  System 

When  conducting  aerodynamic  tests  in  a  free-flight  spark  range,  the  time-position- 
attitude  history  of  the  test  item  is  precisely  measured  at  various  locations  along  its 
trajectory.  These  precise  position  and  attitude  measurements  are  obtained  by  reading  and 
numerically  coding  the  positions  of  the  test  item’s  shadow  with  respect  to  the  range 
reference  system  as  derived  from  photographs.  The  determination  of  these  positions  and 
attitudes  from  the  shadow  measurements  will  be  explained  in  detail  in  Section  V. 

The  range  reference  system  in  the  ARE  consists  of  four  Kevlar®  fiber  bundles 
(henceforth  called  wires)  with  reference  beads  positioned  at  ±45.72  cm  about  1.52  meter 
centers  over  the  entire  length  of  the  instrumented  range.  Two  of  these  wires  are  strung  in 
front  of  the  wall-mounted  reflective  screens  and  two  are  strung  below  the  ceiling- 
mounted  reflective  screens.  With  the  positioning  described  above,  eight  beads  on  each  of 
four  wires  are  located  in  front  of  each  of  the  131  potential  station  locations  (windows). 

At  each  of  the  50  windows  which  house  fully  operational  dual  plane  shadowgraph 
stations,  two  beads  on  each  of  two  wires  are  in  the  field  of  view  of  each  camera. 

Once  the  wires  and  reference  beads  have  been  positioned  and  calibrated,  see  Section 
3.5.3,  they  are  maintained  in  this  position  by  the  alignment  system^ '.  The  zero-zero-zero 
coordinate  of  the  ARE  axis  system  corresponds  to  the  first  bead  on  the  lower  wall  wire  at 


10 


window  1  (there  is  no  shadowgraph  station  at  window  1).  Emanating  from  this  bead  and 
pointing  downrange  is  the  X  axis,  pointing  to  the  left,  looking  downrange,  is  the  Y  axis, 
and  Z  is  vertically  up.  A  sketch  of  the  range  axis  system  is  shown  in  Fig.  2. 1 . 

2.2.2  Aerodynamic  Axis  Systems 

2.2.2.1  Fixed  Plane  Equations  of  Motions 


The  fixed  plane  axis  system  is  defined  as  shown  in  the  sketch  below. 


X,Y,Z  Earth  fixed  coordinates,  parallel  to  the  ARF  axis  system,  Fig.2.1 
x',y',z'  First  rotation  about  the  Z  (ypp)  yields  intermediate  coordinate  system 
x,y,z  Second  rotation  about  y'  (Gpp)  yields  fixed  plane  coordinate  system 

Where  u,v,w  are  the  velocities  along  the  missile  axes  x,y,z  ;  p,q,r,  are  the  angular  rates 
about  the  x,y,z  axes  ;  and,  l,m,n  are  the  moments  about  the  x,y,z  axes.  The  total  velocity 
of  the  missile  ,V,  is  therefore. 


V  =  (u^  +  v^  +  w^)'^^  2.7 

The  fixed  plane  coordinate  system  is  positive  in  the  same  direction  as  the  Earth  fixed 
system  X,Y,Z  of  the  ARF  as  shown  in  Fig.  2.1 .  Or,  x  is  positive  downrange,  y  is  positive 
to  the  left  looking  downrange  ,  and  z  is  positive  up.  The  missile  angles,  0n,  and  \j;,„  ,  are 
defined  as  shown  in  the  sketches  below. 
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ctive  Screen 
ng,Station.3 


Figure  2.1  ARF  Axis  System 


z 


Note  that  if  the  trajectory  (x)  is  parallel  to  the  X  axis  (flight  parallel  to  the  range 
alignment  system)  6^  is  equal  to  0pp .  Also,  since  the  direction  cosines  (mg,  rig, 
pg )  are  what  are  actually  measured  from  the  shadowgraph  images,  then 

0FP  =  -sin ‘(iig)  I 

and  r  2.8 

\j/FP  =  sin’  [  mg  /  (  mg^  +  Pg^  )'^  ]  J 


The  missile  angles  0^  and  xj/m  are  then  obtained  by  subtracting  the  instantaneous 
trajectory  angles  from  0fp  and  \\>fp .  It  should  be  recognized  that  during  the  fitting  process 
the  instantaneous  trajectory  will  vary  slightly  and  hence  the  missile  angles  0^  and  v(;„  will 
also  vary  with  each  iteration. 
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The  roll  orientation ,  (j) ,  is  defined  in  the  sketch  below, 


The  aerodynamic  roll  angle ,  y ,  is  therefore  defined  as, 

Y  =  (t>-(t)i  2.9 

where  (j»i  =  tan’*  ( -v  /  w) 

Normally  the  aerodynamic  roll  angle  is  calculated  in  the  analysis  routine  using  the 
following  relation, 

Y  =  tan’’ (  Vb/ Wb)  2.10 

where 

Vb  =  V  cos  <|)  +  w  sin  (|) 

Wb  =  w  cos  ([>  -  V  sin  (j) 

Note  that  when  the  reference  fin  ( or  the  spin  pin )  is  vertically  down  (|)  is  equal  to  zero  ; 
and,  the  missile’s  measured  roll  orientations  will  be  between  0  and  360  degrees.  During 
the  process  of  fitting  the  theoretical  equations  of  motion  to  the  experimental 
measurements  a  continuous  roll  orientation  profile  is  required.  Hence,  we  are  required  to 
determine  when  the  roll  passes  through  each  360  degree  revolution.  This  ‘unwinding’  of 
the  raw  roll  angles  can  be  one  of  the  more  difficult  tasks  to  accomplish  when  analyzing 
fi'ee-flight  spark  range  trajectories.  Although  tools  have  been  developed  and  incorporated 
into  the  analysis  routines  to  assist  the  analyst  in  accomplishing  this  task,  it  can  still  be 
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difficult  when  the  roll  rate  is  rapidly  changing  or  when  full  spin  is  not  imparted  by  the 
barrel’s  rifling. 


The  6DOF  fixed  plane  equations  of  motion  have  been  derived  in  Refs.  2.2  and  2.3  and 
are  listed  below.  These  equations  include  the  Coriolis  accelerations  and  recognize 
that  0  is  0n,  and  \\i  is  \|/n,  in  these  relations. 

\ 

X  =  u  cos  ^cos  y/  -  V  sin  +  w  sin  ^  cos  y/ 
y  =  u  cos  ^sin  y/  +  \  cos  y/  +  w  sin  ^sin  y/ 
z  =  -  u  sin  ^  +  w  cos6 


d  =  q 
y/  =  X  cos  9 
^  =  p  +  rtan^ 

ii  =  g  sin  ^  -  q  w  +  r  v  -  a^„  +  /  m^ 

v==  -ru-rw  tan^  -  a^^  +  F^  /  m^ 
w  =  -  g  cos^  +  r  V  tan^  +  q  u  -  a^^  +  F^  /  m3 

P  =  1/Ix 

q=  -  rMan6>  -  rp  //^  +  m/ly 
f  =  qrtan(9  +  qp  +  n/I^ 


>  2.11 


where  :  m^  =  projectile  mass 

g  =  acceleration  due  to  gravity 

The  Coriolis  effects  included  above  are  defined  as, 

a^u  =  Ucx  cos  0  cos  v  +  a^y  cos  0  sin  v|/  -  a^^  sin  0 
^cv  “  -acxSin  v|/  +  a^yCos  y 
a^w  =  a^xSin  0  cos  \\i  +  a^ySin  0  sin  \|/  +  a^zCos  0 
and 

a^^  =  -  2  Q}^  (y  sin  z  cos  sin  5^  ) 
a^y  =  2  o)^(x  sin  Ar  -  z  cos  ArCos  J  r  ) 
a.^  =  2  (o^(x  cos/Ir  sint^R  +  y  cosArCos  Jr) 

where  : 

©e  is  the  Earth’s  rotation  rate ,  7.272  x  10'^  rad  /  sec 
Xr  is  Earth  latitude  of  test  facility ,  30.5  deg  for  ARF 
6r  is  Earth  azimuth  of  test  facility ,  126.0  deg  for  ARF 
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The  above  equations  of  motion  are  applicable  for  symmetrical  spin  stabilized  or  fin 
stabilized  free-fiight  missile  configurations.  Because  of  the  symmetry  assumptions  the 
missiles  must  have  at  least  three  fins. 

2.2.2.2  Body  Fixed  Equations  of  Motion 

The  body  fixed  axis  system  is  simply  the  fixed  plane  system  defined  in  the  previous 
section  but  is  allowed  to  roll  with  the  missile.  Or, 


This  coordinate  system  ( Xb ,  yj, ,  Zj, )  is  defined  with  Xb  aligned  with  the  longitudinal  axis 
of  the  missile  and  points  out  the  nose.  The  yb  axis  points  out  the  left  wing  (or  equivalent) 
and  the  Zb  axis  points  up  with  respect  to  the  body.  The  body  fixed  coordinate  system  is 
rigidly  attached  to  the  missile  and  rotates  with  the  missile  about  the  Xb  axis  as  shown 
above. 

The  equations  of  motion  in  this  system  can  then  be  defined  as  shown  below,  see 
Refs.  2.4  and  2.5.  Again  0  =  Bm  and  ij/  =  \|/^  in  these  relations. 
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X  =  Ub  cosOcosy/  +  VbCsin^sin^cos^^/'  -  cos^z^sin^^/') 
+Wb(sin^cos^cos^  +  sin^sin^) 
y  =  Ub  cos^sin^  +  sin^siruz/'  +  cos^cos^) 

+Wb  (  sin^cos^sin^^  -  sin^cos^z/') 
z  =  -Ub  sin^  +  Vb  cos6'  sin^  +  Wb  cos^cos^ 

^  =  qb  cos^  -  tb  sin^ 

^  =  ( Qb  sin^  +rb  cos^)  /  cos^ 

^  =  Pb  +  tan^(  Qb  sin^  +  tb  cos^) 

“b  =  rb^b  -  qbWb  -  acd,  +(Fxb  !  >",)  +  gsin^ 

Vb  =  Pb Wb  -  +( /  m  J  -g  sin«>cos6> 

Wb  =  Qb^b  -  PbVb  -  acub  +(Fzb  /  )  -gcoS(!>cos6' 

Pb  =[Iy*b  +Ix,mb  -  (I.+Iy-IJI^Pbfb]  /  (1,1,-!^") 
+  [(1^"  +  Ib^b]  !  (Ixly  -  Ixy") 

Qb  “Uxl^b  +  Ixy'b  +(Ix+Iy-Iz)Ix,qbrb]  /  (Ixly  "Ixy^) 

+  [(I,(I,-Iv)-V)Pbrb]/axIy-Ix,') 

fj  =  [nb+I^(Pb^-  qt^)  +  (Ix-  Iy)Pbqb]  ! 


Note  that  the  Coriolis  accelerations  are  also  included  in  the  body  fixed  equations  of 
motion.  The  moments  of  inertia  about  the  x,y,z  axes  are  Ix ,  ly ,  and  respectively.  There 
is  only  one  product  of  inertia,  Ixy ,  because  it  is  assumed  that  there  is  one  plane  of 
symmetry,  i.e. ,  Ixy  =  lyz  •  If  the  model’s  principal  axis  is  aligned  with  the  body  axis  then 
Ixy  is  also  equal  to  zero.  This  could  be  a  symmetric  body  with  wings  or  a  body  with  an 
oval  shaped  cross  section  with  or  without  wings. 


2.12 
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2.3  AERODYNAMIC  MODELS 


2.3.1  Body  Fixed 


As  previously  stated  the  interaction  between  the  atmosphere  and  the  missile  results  in^. 
forces  and  moments  acting  on  the  missile.  These  forces  and  moments  are  illustrated  in  the 
sketch  below  for  an  asymmetric  configuration. 


2 


Using  the  aerodynamic  nomenclature  as  defined  in  Section  2. 1  the  forces  can  be  written 
as, 


Fx  =  -(l/2)p  V  A  Cx  +  ma  g  sin  fipp 
Fy  =  (1/2)  p  V^A  Cy  -  ma  g  cos  0pp  sin  (j) 
Fz  =  (1/2)  p  V^A  Cz  -  ma  g  cos  9pp  cos  (|) 


y  2.13 


The  aerodynamic  coefficients  Cx ,  Cy ,  and  Cz  are  expanded  as  functions  of  the  total 
angle  of  attack ,  a ,  or  the  components  of  the  total  angle  of  attack  ,a  and  p,  as  shown  in 
the  sketch  above.  These  expanded  coefficients  can  also  be  assumed  to  be  functions  of 
Mach  number  and  the  aerodynamic  roll  angle.  The  nonlinearities  with  respect  to  a  or  a 
and  P  have  been  assumed  to  be  functions  of  their  respective  sines,  i.e., 

sin  a  =  ( Wb^  +  Vb^  )*^  /  V  =8 
sin  a  =  Wb  /  V 
sin  P  =  Vb  /  V 
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The  resulting  expansions  for  Cx  ,  Cy ,  and  Cz  are  then, 

Cx  =  Cxo  +  Cxa  (  Wb  /  V  )  +  Cxa2  (  Wb  /  V  )^  +  Cx|32  (  ^b  /  V  )^  +  Cxm  (  H  ’  Mref  )  2.15 

CY  =  CY0  +  CYp(Vb/V)  +  CYp3(Vb/V)'  +  CYp„(pd/2V)(Wb/V) 

+  CYya  s^(Wb/V)cos(Ny)  2.16 

Cz  =  Czo  +  Cza(Wb/V)  +  Cza2(Wb/V)'  +  Cza3(Wb/V)' 

+  Czpp(pd/2V)(Vb/V)-Czya  e'(Vb/V)cos(NY)  2.17 

and  N  =  number  of  planes  of  symmetry.  Substituting  these  expansions  for  Cx ,  Cy ,  and 
Cz  into  Eq.  2.13  yields  the  associated  force  definitions. 

Again  using  the  moment  definitions  of  Section  2.1, 

1  =  (1/2)  p  A  d  C, 

m  =  (l/2)pV^AdC„, 

n  =  (l/2)pV^AdCn 

and  similarly  expanding  Cj ,  C^  ,  and  C^  as  we  did  for  the  associated  forces  ,we  define 
C,  =  C,o  +  C,p  (pd  /  2V )  +  C,p  ( Vb  /  V )  +  C,,a  2  s'  sin  (Ny)  2.18 

C„,  =  C„,0+  Qa(Wb/V)  +  C„,^(Wb/V)'  +  C,„3(Wb/V)' 

+  a  3  (  Wb  /  V  )  sin  (Ny)  +  C^q  (  Qb  d  /  2V) 

+  C,qa2  ( Wb  /  V )'  (qbd  /  2V  )  +  c,,a  3  s'  ( Vb  /  V  )  sin  (Ny)  2.19 

C„  =  C„o  +  C„p(Vb/V)  +  C„p3(Vb/V)'  +  C„,(rbd/2V) 

-C„ya3s'(Vb/V)sin(Ny)  +  C„,p2(Vb/V)'(rbd/2V) 

+  C„ya3  8'(Wb/V)sin(Ny)  2.20 

Again  substituting  Eqs.  2.18  through  2.20  into  the  moment  definitions  above  one  arrives 
at  the  general  moment  expansions. 

It  should  be  noted  that  the  general  force  and  moment  expansions  (Eqs.  2. 1 5  to  2.20) 
are  modified  and  extended  as  required.  When  analyzing  experimentally  measured 
trajectory  data  the  analyst  should  strive  to  fit  the  measurements  to  the  experimental 
accuracy  of  the  data.  This  may  necessitate  modifying  the  force  and  moment  expansions 
above.  Knowing  what  expansions  and  terms  to  extend  or  modify  can  be  difficult  at  best 
and  usually  requires  considerable  experience  and  insight. 
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2.3.2  Fixed  Plane 


Unlike  the  body  fixed  axis  system  which  rotates  with  the  body,  this  system  does  not 
rotate  but  rides  along  with  the  missile  and  is  defined  as  previously  discussed  in  Section 
2.2.2.  However,  since  the  missile  is  symmetric  with  at  least  3  planes  of  symmetry,  certain 
of  the  aerodynamic  force  and  moment  coefficients  are  equal. 


The  forces  and  moments  are  expanded  similarly  as  the  body  fixed  forces  and  moments 
in  the  previous  section.  However,  the  nonlinearities  for  this  case  are  assumed  to  be  only  a 
function  of  the  sine  of  the  total  angle-of-attack  ,  a  ,  and  not  the  components  of  the  total 
angle-of-attack  as  was  for  the  body  fixed  case.  Again, 

sin  a  =  (  v^  +  w^)^^  /  V  =  8  s 

where  : 


Note  that 


sin  0n,  =  w  /  ( u^  +  w^)*'^ 
sin  v|/^  =  -  V  /  V 


\  2.22 

J 


sin^  a  ^  sin^  0n,  +  sin^  v|/n, 
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The  force  equations  are  therfore  written  as  : 


Fx  -  -  (1/2)  p  A  [  Cxo  +  Cx  cr  2  Cx  a  4  ^xm  ( ) 

+  Cxya2S^cos(NY)  2.23 

FY  =  (l/2)pV'A[-CNa(v/V)-CNa3e'(v/V)-CNa5e'(v/V) 
-CNy«3e'(v/V)cos(Ny)  +  CYpa(pd/2V)(w/V) 

+  CYpa3s'(pd/2V)(w/V)  +  CYya3e'(w/V)sin(Ny) 

+  CN5  5Asin^  -  Cn8  6bCos(1)]  2.24 

Fz  =  (l/2)pV'A[-CNa(w/V)-CN«3e'(w/V)-CNa5e\w/V) 
-CNy53s'(w/V)cos(Ny)  +  CYpa(pd/2V)(v/V) 

+  CYpa  3  s'(  pd  /  2V )  ( V  /  V )  -  Cy^  « 3  e'C  v  /  V )  sin  (Ny) 

-CNg5AC0S(l)  -  CNsSssincl)]  2.25 

and  the  moment  equations  are  : 

l  =  (l/2)pV'Ad[C,o  +  (pd/2V)C,p  +  (pd/2V)C,pa2e" 

+  (  pd  /  2V  )  C,pM  (  Mi  -  Mref  )  +  C„a  2  e'  sin  (Ny)  ]  2.26 

m  =  (l/2)pV2Ad[C„a(w/V)  +  C^a3s"(w/V)  +  C,a5S^w/V) 

+  C„,  a  M  ( Mi  -  Mref  )  (  w  /  V )  +  « 3  (  w  /  V )  cos  (Ny) 

+  C„q(qd/2V)  +  C„,«2  £"(qd/2V)  +  C,qa4e\qd/2V) 

+  C„a(v/V)  +  C„pa(pd/2V)(v/V)  +  C„p«3s"(pd/2V)(v/V) 

+  C„pa5e^pd/2V)(v/V)  +  C„,a3e'(v/V)sin(Ny) 

+  Cm8SAC0S(|)  -  0^563  sin  (j)]  2.27 

n  =  (l/2)pV"Ad[-C„5(v/V)-C„53s"(v/V)-C„ase‘'(v/V) 
-C„SM(Mi-MKEF)(v/V)-C„53e"(v/V)cos(Nr) 

+  C„,(rd/2V)  +  C^52S=(rd/2V)  +  C„,54s‘(rd/2V) 

+  C,a  (w/V)  +  C„pS(P<l/2V)(w/V)  +  C„pS3s"(pd/2V)(w/V) 

+  C„p558*(pd/2V)(w/V)  +  C„,53e^(w/V)sin(Nr) 

+  Cn,8  5a  sin  (j)  +  Cn,5  6b  COS  (j)  ]  2.28 


Again  N  =  number  of  planes  of  symmetry  ( i.e.  number  of  fins )  and 

y  =  tan‘‘  (vg/Wb) 

where 

Vb  =  V  cos  (j)  +  w  sin  (|) 

Wb  =  -vsin(|)  +wcos(|) 


S  2.29 
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2.4  IDENTITIES 


2.4.1.  Angular  Relationships 
Missile  Angles 

sin  a„,  =  (  sin  +  cos  ii/^  sin  Bn, ) 

sin  e„,  =  sinx„/cosvi/„ 

sin  Gn,  =  sin 

sin  =  -  sin  cos  Xm 

tan  Pn,=  -tanij/m  /  cosGn, 

Missile  Angles  to  Earth  Angles 

0  =  0m  -  sin'^  ( sin  /  cos  ij/m ) 

\\i  =  sin''  ( sin  v)/m  /  cos  )  +  5^ 

Direction  Cosines  ( ,  n,. ,  )  to  Fixed  Plane 

0  =  -  sin'*  ng 

\|/  =  sin  [mg/(mg  +Pe  )  ] 

2.4.2  Velocity  Relationships 

Missile  Velocities  to  Missile  Angles 

V  =  -V  sin  \|/m 
w  =  V  cos  vj/m  sin  0^ 

U  =  V  cos  V|/m  cos  0m 

Body  Fixed  to  Fixed  Plane 

U=Ub 

v  =  V|,cos(|)  -  Wb  sin(|) 
w  =  Vb  sin  <|)  +  Wb  cos  (j) 

P  =  Pb 

q  =  qb  cos  <|)  -  Tb  sin  <j) 
r  =  qb  sin  (|)  +  ib  cos  <j) 

Fixed  Plane  to  Body  Fixed 

Ub  =  U 

Vb  =  vcos(|)  +  wsinij) 

Wb=-vsin<t)  +  wcos(|) 

Pb  =  P 

qb  =  q  cos  <))  +  r  sin  (|) 
rb  =  -qsin(|)  +  rcoscj) 
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2.43  Missile  Angle  Definitions 


ANGLE 


SINE 


COSINE 


TANGENT 


il  2  2 

u  +  V  +  w 


U  +  V 

/T~T~T 

a  +  V  +  w 


^1.1 
U  +  V 


/2V2~2 

a  +  V  +  w 


/2T2 

u  +  w 


U  +  V  +  w 


/2~2 
u  +  w 


/2~2 
u  +  w 


/  2  .  2 
u  +  w 


/  2  .  2 
u  +  w 


/  2  .  2 
u  +  w 


U  +  V 


/  2  ^  2 
U  +  V 


/2~2 
V  +  w 


/  2  .  2  ^  2 
u  +  V  +  w 


/  2  .  2  ^  2 
U  +  V  +  w 


/  2  ^  2 

V  +  w 


2.4.4  Trajectory  Angle  Definitions 


ANGLE 


SINE 


COSINE 


TANGENT 


2  2  ‘ 
y  ? 

V 


2  2  % 

(i  +y  ) 


2  2^ 

(i  +  7  ) 


_ 

(i  -H  7  ) 
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2.4.5  Coordinate  System  Transformations 

Fixed  Plane  -  Missile  to  Earth 

X  -  u  cos  0  cos  v|/  -  V  sin  v|/  +  w  sin  0  cos  v|/ 

=  u  cos  0  sin  V)/  +  V  cos  i|/  +  w  sin  0  sin  i}/ 
i  =  -u  sin  0  +  w  cos  0 

Fixed  Plane  -  Earth  to  Missile 

\i=  X  cos  0  cos  y  +  y  cos  0  sin  v|/  -  i  sin  0 
V  =  -X  sin  ij/  +  y  cos  v|/ 
w  =  X  sin  0  cos  \|/  +  y  sin  0  sin  i|/  +  z  cos  0 

Body  Fixed  -  Missile  to  Earth 

X  =  Ub  cos  0COS  v|/  +  Vb(sin  0sin  (|)cos  v|/  -  cos  <j)sin  \\i)  +  Wb(sin  0cos  (t)COS  v|/  +  sin  (t)sin  \\i) 
=  Ub  cos  0  sin  v(/  +  Vb(sin  0sin  (|)sin  v(/  +  cos  (j)cos  vj/)  +  Wb(sin  0cos  <|)sin  v|/  -  sin  (j)cos  v]/) 

z  =  -Ub  sin  0  +  Vb  cos  0  sin  (fi  +  Wb  cos  0  cos  (j) 

Body  Fixed  -  Earth  to  Missile 

Ub  =  X  cos  0  cos  v|/  +  y  cos  0  sin  v|;  -  z  sin  0 

Vb  =  X  (sin  0  sin  (j)  cos  \\f  -  cos  (j)  sin  V|/)-i-y  (sin  0  sin  ([)  sin  \|;  +  cos  <|)  cos  H/)+i  (cos  0  sin  (|)) 

Wb=  X  (sin  0  cos  (j)  cos  v)/  +  sin  (j)  sin  v|/)+j)  (sin  0  cos  (j)  sin  vj/  -  sin  (j)  cos  v|/)+z  (cos  0  cos  ([>) 
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III.  TEST  PREPARATION 


3.1  TEST  OBJECTIVES 

The  first  step  in  solving  any  problem  is  gaining  an  understanding  of  what  the  problem 
is.  This  encompasses  the  two  questions  of :  why  are  the  tests  to  be  accomplished,  and 
what  needs  to  be  determined  ?  The  test  engineer  cannot  plan  and  design  a  test  /  research 
program  without  knowing  the  answers  to  these  basic  questions.  For  example,  if  the 
primary  purpose  of  the  program  is  to  determine  the  allowable  range  of  the  center-of- 
gravity  for  a  particular  configuration  then  the  center-of-pressure  location  must  be 
determined.  As  will  be  shown  in  Section  7.3.1  this  requires  that  both  the  pitching  moment 
and  normal  force  derivatives  be  accurately  measured.  Although  the  pitching  moment 
derivative  is  normally  one  of  the  easiest  and  most  accurately  determined  measurements 
obtained  from  free-flight  trajectories  the  normal  force  derivative  can  be  difficult  to 
precisely  determine.  The  test  engineer  may  therefore  need  to  design  the  model  and  /  or 
test  procedure  in  such  a  manner  as  to  maximize  his  ability  to  determine  this  derivative. 
Since  the  normal  force  derivative  is  primarily  obtained  from  the  swerving  motion  in  the 
y  -  z  directions,  this  could  mean  minimizing  the  model’s  mass  in  order  to  ensure  that  the 
mavimiim  amount  of  swerving  motion  is  obtained.  If  the  test  engineer  has  no  control  over 
the  model’s  mass,  i.e.,  frequently  he  is  handed  the  models  and  told  to  “shoot  these”,  then 
his  only  other  option  to  maximize  the  swerving  motion  is  to  attempt  to  augment  the 
angular  motion  by  some  means.  The  increased  angles  of  attack  increase  the  normal  force 
acting  on  the  test  item  and  hence  increase  the  resulting  swerving  motion.  Methods  of 
augmenting  the  initial  disturbances  will  be  discussed  later  in  this  section. 

Another  example  of  why  the  test  engineer  needs  to  understand  the  purpose  of  the  test 
program  could  be  that  the  drag  data  are  required  for  a  fire  control  system  and  /  or  the 
round  dispersion  is  also  required  for  probability  of  hit  calculations.  For  this  purpose  it 
may  not  be  desirable  to  augment  the  angular  motions  because  this  could  contaminate  the 
primary  data  required  from  the  tests.  The  increased  angles  of  attack  would  require  a 
correction  to  be  made  to  the  measured  drag  coefficients  in  order  to  obtain  the  zero  yaw 
drag  coefficient  needed  for  the  fire  control  algorithm.  These  corrections  can  certainly  be 
made,  see  Section  6.2.1 .  However,  they  can  potentially  add  an  error  source  into  the  Cdq 
results.  Also,  the  round  dispersion  is  a  function  of  the  initial  angles  of  attack,  see  Section 
7.4.1,  and  artificially  increasing  the  angular  motions  will  likewise  increase  the  resulting 
dispersion. 

The  purpose  of  this  discussion  is  to  underscore  the  need  for  the  test  engineer  to 
understand  the  purpose  of  the  tests  and  what  data  /  measurements  are  required  to  obtain 
the  needed  information.  This  understanding  between  the  test  engineer  and  the  user  / 
sponsor  should  be  one  of  the  primary  goals  of  the  initial  test  plan  meeting  or  test  request 
correspondence. 
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3.2  PRETEST  PLANNING 


Once  the  test  engineer  understands  what  data  are  required  by  the  user  it  then  becomes 
necessary  to  consider,  early  in  the  planning  process,  how  this  may  affect  the  facility 
configuration,  the  launcher  to  be  used,  and  the  instrumentation  requirements.  Facility 
considerations  may  include  whether  or  not  the  complete  range  is  to  be  used  or  if  the  flight 
is  to  be  terminated  early.  Many  flights  are  terminated  prior  to  the  end  of  the  instrumented 
seetion  of  the  facility.  Usually  these  flights  are  terminated  because  the  drop  in  the 
trajectory  due  to  gravity  means  that  the  model  will  fall  below  the  instrumentation  window 
( i.e.  below  the  camera’s  field  of  view  or  the  IR  detection  system).  These  flights  are 
therefore  terminated  before  the  model  impacts  any  of  the  instrumentation,  i.e.  cameras 
mounted  in  the  pits  on  the  floor  of  the  facility.  Other  flights  may  be  terminated  early  to 
avoid  hitting  the  facility  walls  caused  by  high  dispersion  or  high  lift  eonfigurations. 

A  portable  projectile  trap  is  used  to  terminate  these  flights  and  installation  requires 
scheduling  a  fork  lift  and  ensuring  that  the  armor  plating  is  sufficient  to  stop  the  test  item. 
There  may  also  be  considerations  assoeiated  with  any  potential  blast  chamber 
modifications.  These  may  encompass  the  size  and  loeation  of  the  stripper  hole  leading 
fi-om  the  blast  chamber  into  the  instrumented  range,  or  any  required  witness  panels  in  the 
blast  chamber,  or  range  proper.  These  preparations  may  require  considerable  lead  times  in 
order  to  procure  or  manufacture  the  needed  components.  Hence,  the  planning  for  such 
preparations  should  be  accomplished  as  early  as  possible. 

The  launcher  to  be  used  for  the  proposed  tests  should  also  be  identified  early  in  the 
planning  process.  The  general  capabilities  of  the  available  launchers  are  shown  in  Fig. 

3.1 .  Once  the  test  conditions  are  determined  the  appropriate  launcher  can  be  seleeted 
using  this  figure.  Where  the  figure  indicates  that  more  than  one  launeher  is  capable  of 
achieving  the  desired  test  conditions,  considerations  such  as  simplicity  of  design,  ease  of 
operation,  or  minimizing  launeh  loads  should  be  involved  in  the  decision  as  to  which 
system  is  to  be  used.  For  example  the  two-stage  light  gas  gxin  generally  has  a  lower  in 
bore  peak  acceleration  than  the  single  stage  powder  guns.  However,  it  is  a  much  more 
complex  and  costly  system  to  operate  and  requires  several  additional  expendables  other 
than  the  normal  model,  sabot  and  propellant  needed  for  the  single  stage  guns.  These 
additional  expendables  include  pistons,  diaphragms,  and  Helium.  Nevertheless  the  trade 
off  between  lower  acceleration  and  complexity  are  the  kinds  of  decisions  that  must  be 
considered  when  selecting  the  system  to  be  used. 

For  some  test  conditions  Fig.  3.1  indicates  that  there  is  only  one  available  system 
capable  of  aehieving  the  desired  launch  conditions  and  this  certainly  simplifies  the 
decision  process.  Also,  the  figure  illustrates  that  there  are  many  launch  conditions  that 
cannot,  at  present,  be  achieved  using  any  of  the  available  systems,  (large  mass  at  high 
velocities).  We  should  always  be  interested  in  investigating  new  and  advanced  systems 
which  have  the  potential  of  expanding  our  capabilities  into  these  areas,  i.e., 
electromagnetic  launchers,  wave  guns,  or  ram  accelerators. 
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Figure  3.1  Launch  Capability  of  Available  Gun  Systems 


The  requirements  for  any  special  instrumentation  systems  should  also  be  considered 
very  early  in  the  planning  process.  This  could  include  Doppler  radar,  interferometry, 
schlieren,  or  any  other  specified  systems  which  require  advanced  preparation  and 
planning.  Prior  thought  into  these  various  aspects  of  the  proposed  tests  will  improve 
significantly  the  chances  of  success  and  save  time  and  resources. 

3.3  MODEL  /  SABOT  DESIGN 

Some  outstanding  tools  have  been  developed  and  provided  to  the  test  engineer  in  order 
to  assist  in  designing  models  and  sabots.  One  of  these  tools  is  the  Projectile  Design  and 
Analysis  System ,  PRODAS  It  is  not  our  purpose  herein  to  discuss  all  of  the 
capabilities  of  this  system  or  to  instruct  the  reader  in  its  use.  However,  it  will  be  used  to 
illustrate  the  general  applications  and  some  of  the  considerations  of  which  the  designer 
must  be  alert  to.  Figure  3.2  shows  a  simple  spin  stabilized  configuration  designed  using 
this  system.  The  total  model  mass  properties  along  with  the  elemental  mass  properties  of 
this  design  are  also  shown  in  this  figure.  The  projectile  shown  was  designed  using  only 
one  material,  aluminum,  with  no  fins  or  internal  cavities  but  these  do  not  represent 
restrictions  within  PRODAS.  This  routine  can  be  used  to  design  very  complex 
configurations  with  many  elements,  using  different  materials,  with  or  without  cavities. 
The  primary  restriction  on  PRODAS  is  that  only  symmetrical  projectiles  with  circular 
cross  sections  can  be  designed.  PRODAS  has  many  other  capabilities  helpful  to  the 
designer  /  test  engineer  and  some  of  these  will  be  mentioned  in  the  following  sections. 
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The  location  of  the  model’s  center-of-gravity  (eg)  is  one  of  the  major  design 
considerations.  Obviously  the  model  must  be  either  spin  stabilized  or  aerodynamically 
stable  in  flight,  see  Section  7.2,  because  unlike  wind  tunnel  tests  the  model  must  possess 
static  stability  to  successfully  fly  through  the  instrumented  facility.  Also,  it  can  be 
advantageous  to  match  the  eg  location  of  the  full  scale  vehicle.  Although  as  mentioned 
previously,  the  pitching  moment  can  be  shifted  to  any  moment  reference  location  using 
the  relations  derived  in  Section  7.3.2  this  shift  can  induce  errors.  Especially,  when  large 
shifts  are  required  or  when  the  normal  force  data  are  ill  defined.  Some  of  the  other 
derivatives,  i.e.  Magnus  and  damping-in-pitch,  cannot  be  easily  shifted  and  in  fact  usually 
are  not.  Therefore  the  data  obtained  about  the  model’s  eg  location  is  assumed  to  be 
representative  of  the  full  scale  vehicle.  Hence,  if  the  model’s  eg  is  designed  to  be 
consistent  with  the  full  scale  vehicle  these  potential  error  sources  are  eliminated. 

Normally  for  most  test  configurations  it  is  good  design  policy  to  minimize  the 
projectile  mass.  This  not  only  improves  the  precision  of  the  normal  force  as  previously 
mentioned  but  many  of  the  other  coefficients  as  well.  For  example,  the  accuracy  of  the 
drag  measurements  are  improved  by  increasing  the  model’s  deceleration  experienced 
during  the  flight  through  the  instrumented  range.  It  should  be  noted  that  like  the  pitching 
moment  derivative  the  drag  measurements  are  normally  very  precise,  less  than  1%  error. 
But  if  the  model  experiences  very  little  deceleration,  i.e.  in  a  reduced  density 
environment,  minimizing  the  model  mass  may  be  advantageous  in  accurately  measuring 
the  drag. 

The  reduced  mass  also  minimizes  the  angular  momentum  of  the  model  thereby 
maximizing  the  aerodynamic  damping  effects  on  the  model’s  motions.  This  in  turn  is 
reflected  in  better  determined  damping  derivatives,  pitch-yaw  and  roll.  Minimizing  the 
model’s  mass  usually  produces  several  other  beneficial  effects.  Firstly,  the  lesser  mass 
reduces  the  acceleration  loads  experienced  by  the  model  during  launch.  This  increases  the 
probability  that  the  model  will  survive  launch  without  any  structural  failures.  It  also 
simplifies  the  sabot  design  which  should  result  in  a  lighter  and  stronger  sabot  since  the 
forces  acting  on  the  sabot  are  likewise  reduced.  The  lighter  masses,  both  model  and  sabot, 
makes  it  easier  to  terminate  their  respective  flights  without  causing  undue  damage  to  the 
blast  chamber  walls  or  the  projectile  trap  at  the  end  of  the  range.  Lastly,  the  reduced 
model  /  sabot  package  mass  requires  less  of  a  propellant  charge  and  hence  less  of  an  over 
pressure  in  the  blzist  chamber  of  the  facility. 

Nevertheless  it  may  not  always  be  wise  to  minimize  model  mass.  When  testing  high 
lift  configurations  a  heavy  model  may  be  desirable  in  order  to  reduce  the  lift  effect  on  the 
model’s  trajectory.  These  configurations  can  exhibit  large  deviations  in  their  flight  paths 
and  culminate  in  the  model  impacting  the  facility  walls  or  even  worse  the  facility 
instrumentation  systems.  One  technique  frequently  used  to  control  the  trajectories  of 
high  lift  configurations  is  to  purposely  induce  a  roll  rate  during  laxmeh  and  /  or  flight. 

This  ensures  that  the  lift  force  is  not  always  pointed  in  the  same  direction  and  tends  to 
cause  the  model  to  fly  in  a  spiral  trajectory.  Obviously  the  goal  is  to  make  the  diameter  of 
the  spiral  less  than  the  diameter  of  the  instrumentation  window  of  the  facility.  Other 
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common  shapes  which  actually  require  heavy  models  are  blunt  entry  configurations,  i.e. 
planetary  probes.  These  configurations  have  such  large  drag  coefficients  and  hence  high 
decelerations  that  the  test  engineer  usually  attempts  to  minimize  the  deceleration  the 
model  experiences  in  the  range. 

Spin  pins  are  frequently  inserted  into  the  model’s  base  in  order  to  measure  the  roll 
orientation,  and  therefore  the  roll  rate  (spin)  as  the  model  traverses  the  instrumented 
range.  See  Section  5.2.3  for  a  discussion  of  these  measurements  and  the  design  criteria 
that  should  be  considered.  When  testing  subscale  models  of  missiles  or  projectiles  that 
operate  at  high  roll  rates  it  is  desirable  to  match  as  best  as  possible  the 
nondimensionalized  spin  parameter,  pd  /  2V  ,  as  mentioned  in  the  Section  1 .4.  This 
simulates  the  rotational  flow  about  the  full  scale  missile  in  flight.  Some  of  the  more 
common  methods  used  to  accomplish  this  match  include  using  rifled  launch  tubes, 
grooved  sabots,  and  canted  fins. 

One  more  consideration  the  model  designer  should  take  into  accoimt  is  the  possibility 
that  aerodynamic  heating  may  cause  the  model  to  ablate  during  flight.  This  ablation 
process  can  cause  shape  changes  to  occur.  These  shape  changes  can  be  especially  severe 
on  the  model’s  nose  or  fin  leading  edges  and  effectively  alters  the  configuration  being 
tested.  This  can  be  a  serious  problem  at  hypersonic  Mach  numbers  ( above  6)  eind  may 
require  certain  heat  resistant  materials  to  be  used  for  these  conditions.  This  ablation  can 
be  predicted  and  materials  selected  to  minimize  this  problem  using  the  curve  presented  in 
Fig.  3.3  in  combination  with  the  data  shown  in  Table  3.1.  Figure  3.3  and  Table  3.1  were 
extracted  from  Ref.  3.2  and  illustrate  that  ablation  is  a  function  of  the  material  properties, 
atmospheric  pressure,  and  flight  distance.  Note  that  ablation  occurs  for  all  conditions 
below  the  curve  in  Fig.  3.3.  Obviously  when  making  predictions  using  this  curve  the 
flight  length  should  be  assumed  to  be  the  length  of  the  range.  By  judiciously  selecting 
materials  the  designer  should  attempt  to  minimize  these  ablation  effects  during  the  test 
and  therefore  eliminate  them  from  the  resulting  aerodynamic  coefficients  and  derivatives. 

At  the  present  time  the  designer  does  not  have  a  user  fi-iendly  software  package  to 
assist  in  designing  sabots,  computing  stress  levels,  and  predicting  potential  sabot  failure 
modes.  Although  there  is  an  on  going  effort  along  these  lines  and  such  a  tool  may  be 
available  in  the  near  future.  PRODAS  can  be  used  to  model  the  sabot  and  compute  the 
mass  properties  and  there  are  several  outstanding  references,  i.e.  Ref  1.3  and  Refs.  3.3 
and  3.4  ,  which  can  also  be  used  to  assist  the  designer.  Unfortunately  some  aspects  of 
sabot  design  are  still  more  of  an  art  than  a  science  but  this  is  rapidly  changing  with  the 
use  of  modem  finite  element  programs.  There  is  a  considerable  body  of  literature  which 
discusses  various  point  designs  also  available^ 
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Figure  3.3  Flight  Conditions  Which  Will  Give  Material 

Melting  Temperatures  at  the  Stagnation  Point 
(  Fig.  26  ,  Ref.  3.3  ) 
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Table  3.1 
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3.4  PRETEST  PREDICTIONS 


It  is  extremely  advisable  to  make  a  concerted  effort  at  predicting  the  various  aspects 
and  potential  results  associated  with  the  test  program.  These  predictions  should  be 
accomplished  prior  to  testing  and  include  such  things  as  the  required  launcher  conditions, 
test  item  aerodynamics,  expected  trajectories,  and  any  other  aspect  pertinent  to  the 
program.  Predictions  are  generally  thought  of  as  computations  of  some  sort  and  as  such 
will  be  discussed  in  the  following  sections.  However,  predictions  can  also  include 
personal  experiences  or  experiences  documented  by  other  individuals  conducting  related 
programs.  For  example  if  good  records  are  available  and  maintained  from  previous  test 
programs  these  can  be  an  excellent  source  of  predictions  for  futvue  tests.  This  is  one  of 
the  primary  reasons  why  all  test  organizations  maintain  test  logs  which  include  such  data 
as  are  listed  below. 

a.  Launcher  conditions  (i.e.) 

bore  diameter 
rifling  rate(twist) 
propellant  type(s)  and  charge(s) 
piston  mass 

diaphragm  rupture  pressure 
pump  tube  gas  and  pressure 

b.  Sabot  description  (i.e.) 

pusher,  puller 
number  of  petals 
mass 
materials 
sketches 

c.  Model  description  (i.e.) 

mass 

moments-of  inertia 
center-of-gravity 
lengths,  diameters 
materials 
dravidngs 

d.  Facility  test  conditions 

temperature 
relative  humidity 
pressure 

e.  Results 

velocity  /  Mach  number 
Reynolds  number 
flight  (length,  time) 
instrumentation  obtained  or  missed 
model  /  sabot  integrity 
general  comments 
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Published  data  reports  and  technical  papers  of  previous  tests  are  also  an  excellent 
source  of  predictions.  Hence  a  literature  search  should  be  conducted  for  archival 
information  of  similar  previous  test  programs.  In  fact  any  new  data  resulting  from  the 
proposed  tests  should  be  compared  with  any  previous  data  of  a  similar  configuration  and 
test  conditions.  Such  comparisons  provide  insights  as  to  the  viability  of  the  test  results 
and  hence  provides  a  measure  of  reliability  to  the  test  methods  and  processes  used.  In 
fact  once  the  tests  are  in  progress  any  and  all  predictions  should  be  compared  with  the 
actual  results  as  soon  as  possible  and  the  differences  analyzed  by  the  test  engineer.  Such 
comparisons  can  yield  valuable  information  which  can  alert  the  test  engineer  to  potential 
problems  or  anomalies  present  during  testing.  These  comparisons  can  also  lead  to 
potential  advances  of  the  existing  predictive  capabilities  which  will  assist  the  designer  in 
future  programs. 

3.4.1  Launcher 

Once  the  test  engineer  has  identified  the  launcher  to  be  used  either  from  Fig.  3.1  or 
other  means  there  exists  numerous  techniques  to  predict  launcher  performance  and  the 
resulting  muzzle  conditions.  Certainly  the  operations  crew  need  to  know  the  required  load 
/  charge  conditions  to  achieve  the  desired  muzzle  velocity.  Likewise  the  model  /  sabot 
designer  needs  to  know  what  in-bore  accelelerations  are  going  to  be  applied  to  the  test 
package.  Depending  on  the  type  of  launcher  selected  different  prediction  methods  are 
used.  For  example  if  a  single  stage  powder  gun  is  used,  PRODAS^ '  is  frequently  used  to 
predict  these  interior  ballistic  characteristics.  A  simple  PROD  AS  prediction  is  shown  in 
Fig.  3.4.  The  model  configuration  shown  in  this  figure  is  the  same  one  shown  in  Fig.  3.2 
when  the  model  design  features  of  PRODAS  were  discussed.  The  interior  ballistic 
calculations  shown  in  Fig.  3.4  represent  an  abbreviated  output  and  much  more 
information  is  available  to  the  user  if  required.  Again  it  is  not  our  purpose  herein  to 
instruct  the  reader  on  how  to  use  this  computer  program  but  only  to  indicate  that  it  can  be 
a  source  of  interior  ballistics  predictions. 

If  the  two  stage  light  gas  gun  is  used  as  the  launcher  of  preference  other  predictive 
tools  are  available.  Some  of  these  two  stage  light  gas  gun  programs  are  complex  and 
requires  considerable  expertise  in  successfully  obtaining  accurate  estimates.  However, 
some  relatively  simple  routines  are  also  available,  for  example,  the  “Light  Gas  Gun 
Program,  LGG”  was  written  specifically  for  the  small  launcher  operated  by  the 
Vehicles  Branch,  Wright  Laboratory,  Eglin  AFB,  FL.  This  code  is  applicable  for  20,  30, 
and  40  mm  launch  tubes  and  only  requires  5  additional  inputs  as  listed  below: 

1 .  amount  of  propellant, 

2.  piston  mass, 

3.  helium  pressure  in  the  pump  tube, 

4.  diaphragm  rupture  pressure, 

5.  model  /  sabot  package  weight. 

All  of  the  other  geometric  aspects  associated  with  the  launcher  are  fixed.  Obviously  the 
estimates  provided  by  this  routine  are  not  as  precise  as  the  more  complex  programs  but 
this  is  more  than  made  up  for  in  ease  and  simplicity  of  operation.  A  summary  plot  of 
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Figure  3.4  Example  of  Interior  Ballistics  Using  PRODAS 


propellant  weight  vs.  muzzle  velocity  is  shown  in  Fig.  3.5  for  convenience.  If  peak 
accelerations  are  required  or  more  specific  information  needed  the  user  is  encouraged  to 
execute  this  simple  routine  or  any  other  routine  he  may  have  access  to. 


PACKAGE  WT,  gms 


Figure  3.5  Propellant  vs  Velocity  Relationship  for  the  Eglin  Light  Gas  Gun 


If  the  test  models  are  large  and  the  launch  velocities  relatively  low  then  a  compressed 
gas  gun  may  be  the  launcher  of  choice  ( see  Fig.  3.1).  Subsonically  the  launcher 
performance  can  be  predicted  using  the  relationship  below: 

Vp=  Ef  [2P,V,(l-K'-^)/(y-l)mj''2  3.1 

where  K  =  (  Ab  L],  +  )  /  Vc 

and  Vp  =  muzzle  velocity  of  the  projectile 

Pc  =  chamber  pressure 
Vc  =  chamber  volume 
Y  =  ratio  of  specific  heats  ( air  =  1 .4 ) 

=  projectile  mass 

A(,  =  launch  tube  cross  sectional  area 

Lb=  launch  tube  length 

Ef=  efficiency  factor  ( assumed  to  be  0.87 ) 
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Equation  3.1  was  obtained  from  Ref.  3.8  and  is  applicable  for  subsonic  launches.  A 
unique  variation  of  this  relationship  was  derived  for  a  double  barreled  compressed  air 
launcher  as  is  shown  in  Ref  3.9.  The  efficiency  factor,  Ef ,  was  experimentally  measured 
as  0.87  and  included  in  Eq.  3.1  during  the  tests  discussed  in  Ref.  3.9. 

If  a  supersonic  launch  is  required  ( M  ^  1.5  )  using  the  compressed  gas  gun  then  Eq. 

3 . 1  should  not  be  used  and  the  compressibility  of  the  gas  must  be  included.  These 
relationships  have  been  developed^  and  a  predictive  routine  ,  GUNPAI,  exists  \vhich 
utilizes  these  relationships.  Example  performance  curves  using  GUNPAI  are  shown  in 
Fig.  3.6  and  illustrate  the  difference  in  performance  when  helium  is  used  in  place  of  air  as 
the  propelling  gas.  Also  shown  in  this  figure  are  experimental  data  for  comparison 
purposes. 

Obviously  any  prediction  should  be  evaluated  using  existing  data  from  previous  tests. 
Especially  when  the  user  has  data  of  similar  in-gun  weights  and  fired  at  nearly  the  same 
velocities.  This  provides  an  excellent  mechanism  for  predicting  launcher  performance  in 
future  tests. 

3.4.2  Aerodynamics  /  Trajectory 

The  test  item’s  expected  aerodynamics  should  also  be  predicted  and  used  to  compare 
with  the  experimentally  determined  coefficients  and  derivatives.  These  predicted 
aerodynamics  can  be  extremely  useful  in  illustrating  the  expected  trends,  ( i.e.,  vs  Mach 
number ,  angle  of  attack,  etc. )  and  can  be  used  as  an  indicator  of  the  experimental 
viability.  The  test  engineer  should  carefully  analyze  any  large  deviations  between 
experiment  and  prediction  to  ensure  that  the  experiment  was  properly  conducted  and  that 
the  results  are  reasonable.  Previous  data  on  similar  configurations  at  similar  conditions 
should  also  be  used  for  this  purpose.  Obviously  the  test  engineer  should  be  concerned  if 
similar  configurations  at  similar  conditions  yield  significantly  different  trends  and 
immediate  evaluation  of  these  differences  is  in  order. 

The  PRODAS  program  discussed  previously  is  capable  of  fulfilling  this  requirement 
for  symmetrical  spin  stabilized  or  statically  stable,  i.e.  finned,  configurations.  Table  3.2 
illustrates  the  PRODAS  aerodynamic  predictions  for  the  example  model  shown  in  Fig. 

3.2.  This  routine  is  based  on  previous  data  and  theoretical  relationships  are  used  to 
extrapolate  between  configurational  and  test  condition  variations.  Hence  if  the  test  item 
configuration  falls  within  the  bounds  of  the  existing  data  reasonable  predictions  are 
normally  obtained. 

Once  the  aerodynamics  have  been  predicted  this  routine  can  also  generate  trajectories 
using  these  aerodynamics  in  conjunction  with  the  model’s  physical  characteristics  and  the 
atmospheric  conditions.  Either  6DOF  or  2DOF  trajectories  can  be  generated  for  most 
launch  conditions.  A  typical  2DOF  trajectory  representing  a  flight  through  the  ARF  is 


37 


FNF  =  Coefficient  of  friction  between  sabot  and  gun 


f  =  Fluid  friction  coefficient 


(PSIG) 


a.  Air  Propellent  (  Fig.  5  ,  Ref.  3.10  ) 


(PSIA) 


Figure  3.6  Comparison  of  GUNPAI  Predictions  with 
Experimental  Results  ,  151mm  I.D., 

3.66  m  liong  Launc^liTube  ,  1  Kg  Package 


Aerodynamic  Predictions  Using  PRODAS 
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presented  in  Fig.  3.7.  This  trajectory  illustrates  the  expected  drop  due  to  gravity  and  the 
velocity  decay  of  the  previously  used  example  model  (see  Fig  3.2  ). 

3.4.3  Initial  Motion  Augmentation 

Frequently  it  is  necessary  to  augment  the  model’s  initial  angular  disturbance  as  it 
exits  the  launch  tube.  This  is  required  because  the  naturally  occurring  disturbance  is  not 
adequate  to  impart  an  angular  motion  large  enough  for  accurate  determination  of  the 
aerodynamic  coefficients  and  derivatives.  Over  the  years  numerous  techniques  have  been 
used  to  increase  the  initial  disturbances  including  modifications  to  the  muzzle  of  the 
launcher,  plates  near  the  muzzle  to  cause  asymmetric  flow  upon  exit,  offsetting  the  model 
in  the  sabot,  etc.  All  of  these  techniques  and  many  others  have  had  some  success; 
however,  quite  often  what  worked  for  one  test  configuration  did  not  work  for  the  next. 
Also  the  amount  of  disturbance  could  not  be  predicted  prior  to  testing.  Many  of  these 
techniques  tended  to  impart  an  undesirable  transverse  velocity  as  well  as  the  desired 
angular  disturbance.  This  transverse  velocity  can  result  in  the  model  impacting  the  walls 
of  the  facility,  or  worse  the  instrumentation  systems  in  the  facility. 

The  problem  outlined  above  still  applies  for  statically  stable  configurations  because  no 
systematic  procedure  has  ever  been  developed  which  is  predictable  and  reliable  for  these 
configurations.  On  the  other  hand  a  reliable  and  predictable  method  has  been  developed 
for  spin  stabilized  projectiles.  This  method  uses  the  model’s  gyroscopic  stability  as  a 
means  of  imparting  the  initial  angular  disturbance  without  a  corresponding  transverse 
velocity.  The  gyroscopic  stability  factor ,  Sg ,  is  discussed  in  some  detail  in  Section  VII 
but  for  our  purposes  here  it  is  sufficient  to  note  that  if  Sg  is  less  than  1 .0  the  model  is 
gyroscopically  unstable  and  the  angular  motion  will  increase.  The  relationship  developed 
in  Section  VII  for  Sg  is  , 

Sg  =  (2I,P')/(IyP  ;r  d^C^J  (seeEq.7.11) 

where  : 

Ix  =  axial  moment  of  inertia 

ly  =  transverse  moment  of  inertia 

P  =  spin  rate 

d  =  model  diameter 

Cma  =  pitching  moment  derivative 

p  =  atmospheric  density 

From  this  equation  it  can  be  seen  that  if  the  atmospheric  density  is  increased  sufficiently 
Sg  can  be  forced  to  a  value  less  than  unity. 

This  represents  the  crux  of  the  technique  and  operationally  is  accomplished  by  flying 
the  model  through  a  relatively  short  tube  ( ~  3  m )  containing  a  high  density  gas.  The 
tube  has  both  ends  sealed  with  a  diaphragm,  the  air  is  pumped  out,  and  a  sufficient 
amount  of  high  density  gas  is  bled  into  the  evacuated  tube.  An  environmentally  safe 
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Subscale  lOSmm  with  Ml  Modified  Fuse 


Figure  3.7  Facility  Trajectory  Predictions  Using  PRODAS 
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Freon  compound  is  used  for  this  piupose.  When  computing  the  required  amount  of  gas, 
using  the  equation  above,  an  estimate  of  must  be  available.  This  can  be  obtained 
either  from  similar  data  or  PRODAS  as  mentioned  previously.  It  should  be  recognized 
that  varies  with  Mach  number  and  the  Mach  number  in  the  tube  will  be  significantly 
different  than  the  Mach  number  in  air  ( for  Freon  about  twice  as  high  )  .  Once  the  model 
pxmctures  the  diaphragm  and  enters  the  tube  the  angular  motion  will  increase  because  it  is 
gyroscopically  unstable.  Upon  exit  from  the  tube  the  augmented  angular  motion  can  be 
measured  and  analyzed  as  discussed  in  the  later  sections  herein. 

3.5  PRETEST  MEASUREMENTS 

It  should  be  obvious  that  whatever  can  be  measured  (i.e.,  model’s  physical 
characteristics,  range  and  test  conditions,  instrumentation  calibration,  etc. )  should  be 
measured  and  recorded  prior  to  any  testing  being  accomplished.  Nevertheless  the  lack  of 
obtaining  such  measurements  causes  more  difficulties  during  the  post  test  analysis  than 
any  other  single  cause.  For  some  complex  full  scale  flight  tests  these  measurements  can 
be  difficult  to  obtain.  For  example  the  atmospheric  conditions,  i.e.  temperature,  pressure, 
wind,  etc. ,  are  all  functions  of  altitude  and  the  model’s  geometry  can  be  variable  as  well, 
i.e.  eg  movements,  fin  deflections,  etc.  Difficult  or  not  these  measurements  are 
absolutely  necessary  in  order  to  completely  analyze  the  test  results.  This  is  just  as  true  for 
the  simpler  case  of  testing  subscale  models  within  an  indoor  facility.  Model 
measurements  and  facility  conditions  are  absolute  requirements  and  not  only  obtained  but 
should  also  be  examined  for  reasonableness  and  consistency  prior  to  testing.  During  this 
pretest  measurement  process  it  is  also  a  good  time  to  obtain  photographs  of  the  models 
and  sabots  for  later  documentation  purposes. 

3.5.1  Model  Physical  Properties 

The  physical  characteristics  of  the  models  include  the  linear  dimensions,  mass,  center- 
of-gravity,  and  moments  of  inertia.  There  are  numerous  methods  and  instruments  by 
which  the  dimensional  measurements  can  be  obtained.  Not  only  should  the  obvious 
measurements  of  length  and  diameter  be  obtained  but  also  any  other  potentially  pertinent 
measurements  such  as  nose  bluntness  or  fin  cant.  Small  variations  in  nose  bluntness  can 
significantly  affect  the  aerodynamic  drag  measurements  and  changes  in  fin  cant  can 
seriously  alter  the  rolling  motion  and  hence  the  overall  flight  dynamics  of  the  test  model. 
The  test  engineer  should  be  alert  to  any  potential  model  variation  which  could  affect  the 
results  and  insure  that  the  required  measurements  are  obtained. 

The  distance  between  the  model’s  nose  and  the  center-of-gravity  (Lgg)  is  usually 
measured  using  a  simple  beam  balance  type  system  as  illustrated  in  the  sketch  on  the 
following  page.  Writing  the  moment  equation  about  the  pivot  point,  A ,  leads  to  the 
relation: 

X,  =  WX2/ma 

from  which  the  models  L^g  distance  can  be  determined.  Note  that  m^  is  the  model  mass 
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and  W  is  the  counter  balancing  weight.  In  application  such  systems  can  become  fairly 
complex^  although,  this  is  the  basic  principle  which  is  normally  applied. 

The  model’s  moments  of  inertia,  transverse  and  axial,  are  usually  obtained  using  a 
torsional  pendulum  system.  Sketches  illustrating  this  principle  are  shown  below. 


Once  the  constant  K  is  determined  by  measuring  the  frequency  of  an  object  with  a  known 
inertia  ( i.e.  a  solid  homogenous  cylinder )  an  unknown  inertia,  for  example  a  test  item, 
can  be  determined  using  Eq.  3.2  and  measuring  its  rotational  frequency  about  the  axis  of 
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the  torsion  bar.  There  are  several  commercially  built  instruments  presently  available  for 
this  purpose^  and  most  use  the  principle  as  outlined  above. 

3.5.2  Atmosphere 

When  accomplishing  free-flight  tests  in  an  indoor  facility  these  measurements  are 
relatively  easy  and  consists  of  temperature,  pressure,  and  relative  humidity.  For  large 
( long )  facilities  it  is  recommended  that  multiple  measurements  along  the  flight  path 
be  recorded  at  about  the  same  time  the  tests  are  conducted.  As  previously  mentioned 
when  testing  on  an  outdoor  facility  or  full  scale  flight  profiles  are  duplicated  these 
atmospheric  measurements  along  with  the  existing  winds  should  be  obtained  throughout 
the  altitude  range  of  interest.  Since  these  are  also  frequently  time  dependent  it  is  desirable 
to  obtain  the  measurements  as  near  the  test  time  as  possible.  A  whole  host  of  commercial 
instruments  are  readily  available  for  these  purposes.  However,  those  which  do  not  require 
human  recordings  and  /  or  interpretation  are  preferred.  Regardless  of  which  instruments 
are  chosen,  these  measurements  are  highly  recommended  because  failure  to  obtain  or 
record  these  parameters  can  make  the  post  test  analysis  extremely  difficult  if  not 
impossible. 

3.5.3  Range  /  Instrumentation  Calibration 

As  pointed  out  in  the  previous  sections  it  is  critical  to  measure  and  record  the  physical 
characteristics  of  the  models,  range  conditions  (atmosphere)  during  the  tests  and  any 
other  parameter  that  could  be  expected  to  affect  the  final  results.  Obviously  the 
instruments  used  to  obtain  these  measurements  should  be  periodically  calibrated  to  ensure 
that  the  measurements  are  reliable.  Most  research  /  test  organizations  have  processes  in 
place  which  accomplish  this  periodic  evaluation  and  these  should  be  faithfully  adhered  to 
because  unreliable  measurements  can  render  the  final  results  suspect  and  hence  useless. 
Some  instruments  cannot  easily  be  sent  to  a  calibration  lab  for  evaluation  and  procedures 
must  be  established  to  calibrate  /  evaluate  these  instruments  /  systems  in  place.  This  can 
be  accomplished  by  periodically  evaluating  a  standard  and  maintaining  a  running  log  of 
the  measurement  which  will  highlight  any  unexpected  variation.  These  periodic 
evaluations  can  be  for  simple  instruments  imbedded  in  a  more  complex  system  or  the 
complex  system  itself  For  example,  the  calibration  of  the  facility  and  the  determination 
of  the  its  overall  accuracy  in  positions,  attitudes  and  time. 

Free-flight  spark  range  calibration  is  generally  defined  as  the  precise  determination  of 
the  spatial  locations  of  the  facility’s  reference  system,  spark  sources,  cameras,  and  other 
critical  components.  In  the  early  days  of  free-flight  spark  ranges  these  determinations 
were  accomplished  by  physically  measuring  all  of  these  critical  components  using 
calibrated  tapes,  optical  scopes,  and  a  myriad  of  special  fixtures  eind  tools,  i.e.  see  the 
discussion  in  Ref  1.2.  More  recently  these  determinations  have  been  done  by 
mathematical  means.  Refs.  3.13  to  3.16.  The  mathematical  technique  developed  in  Ref 
3.13  and  utilized  in  Ref  3.16  essentially  uses  known  points  in  space  to  compute  the 
spatial  coordinates  of  the  reference  beads  and  spark  locations,  see  the  abbreviated 
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discussion  in  Appendix  2.  The  technique  discussed  in  Ref.  3.15  determines  the  biases  at 
each  station  by  maintaining  a  record  of  the  various  fits  to  the  experimental  data  and 
subtracting  these  biases  from  the  raw  measurements  during  analysis.  In  practice  both 
methods  should  be  utilized.  The  more  precise  technique  of  Ref.  3.13  should  be 
periodically  utilized,  once  every  year  or  so,  unless  the  reference  system  is  damaged  by  a 
projectile  impact,  then  it  should  be  accomplished  as  soon  as  the  repair  is  completed.  The 
method  of  Ref.  3.15  should  be  used  on  a  continuing  basis  for  evaluating  the  facility 
accuracy  and  updating  the  station  bias  files.  This  technique  is  very  powerful  and  caution 
must  be  used  to  ensure  that  the  residuals  between  the  experimental  data  and  the  fitted 
trajectories  are  caused  by  random  measurement  noise  and  not  modeling  inadequacies. 

If  the  range  calibration  becomes  suspect  it  can  readily  be  examined  by  launching 
several  standard  spin  stabilized  projectiles  whose  aerodynamics  and  flight  dynamics  are 
known.  If  these  residuals  demonstrate  a  noticeable  trend  (i.e.  a  sine  wave  function)  this 
usually  indicates  an  aerodynamic  modeling  problem.  However,  if  the  residuals  at  an 
individual  station  have  a  tendency  to  always  be  positive  (  or  negative )  then  this  indicates 
a  potential  calibration  problem.  The  decision  as  to  when  a  facility  recalibration  is  required 
and  which  method  to  use  hinges  on  several  factors.  These  include  but  are  not  limited  to: 

1)  the  purpose  of  the  tests  and  whether  or  not  the  required  results  would  be  noticeably 
impacted  by  the  suspected  errors,  2)  the  magnitude  and  quantity  of  the  suspected  errors, 

3)  how  much  time  is  available  for  the  calibration. 
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IV.  TESTING 


4.1  SAFETY  /  OPERATING  INSTRUCTIONS 

Safety  should  always  be  of  primary  concern  and  rigorously  following  established 
procedures  will  minimize  if  not  eliminate  hazardous  situations.  Considerable  effort  has 
been  expended  over  the  years  in  developing  well  thought  out  procedures  and  the 
establishment  of  operating  instructions  ( OI’s )  for  the  test  facilities,  launchers, 
instrumentation  systems,  and  freon  tube.  These  OI’s  not  only  discuss  safety  concerns  but 
also  specifically  provide  the  standard  operating  procedures  for  the  facility  and  its 
systems.  Faithfully  following  these  guidelines  will  result  in  a  safe  and  successful  test 
program.  A  hurried  chaotic  process  is  inherently  unsafe  and  leads  to  errors  which 
culminate  in  test  objections  not  being  accomplished. 

When  situations  arise  which  are  not  adequately  covered  by  the  existing  OI’s  these 
should  be  noted  and  included  in  periodic  reviews  and  amendments  to  the  instructions. 
These  situations  occur  because  new  instrumentation  and  launcher  systems  are  developed 
or  procured  and  new  test  requirements  mandate  operations  which  had  never  before  been 
attempted.  Nevertheless,  these  experiences  should  be  included  in  the  OI’s  as  soon  as 
practical  and  while  memories  are  fresh. 

4.2  PRELIMINARY  EVALUATIONS 

Frequently  it  is  necessary  to  conduct  a  preliminary  evaluation  of  the  models,  sabots  or 
instrumentation  prior  to  the  primary  test  program.  These  preliminary  evaluations  are 
normally  conducted  on  an  outdoor  range,  i.e.,  the  Ballistic  Experimentation  Facility 
( BEF  ).  The  desire  to  accomplish  such  tests  in  an  imconfined  environment  arises  from 
the  inherent  uncertainties  associated  with  such  tests  and  for  protecting  the  facility  and  the 
instrumentation  systems  from  unnecessary  risks.  Occasionally  these  preliminary  tests 
must  be  done  within  the  confines  of  an  indoor  range  because  that  is  where  the  required 
systems  are  and  the  impracticality  of  moving  or  operating  these  systems  in  an  outdoor 
environment.  Nevertheless  these  preliminary  evaluations  should  be  accomplished  in  a 
manner  such  that  the  attendant  risks  are  minimized. 

4.2.1  Models  and  Sabots 

It  is  good  practice  to  launch  new  model  /  sabot  packages  first  on  the  outdoor  range  in 
order  to  evaluate  the  performance  of  the  package.  This  includes  the  separation  of  the 
sabot  petals  from  the  model,  the  amount  of  separation  at  a  distance  equivalent  to  the 
stripper  or  blast  chamber  wall,  the  trajectory  of  the  petals,  and  the  structural  integrity  and 
trajectory  of  the  model  itself.  Even  model  and  sabot  designs  that  have  been  used  on 
previous  occasions  should  be  reevaluated  prior  to  a  follow  on  test  entry.  Sometimes 
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materials  are  substituted  for  in  the  manufacturing  process  or  the  machine  tolerances  are 
different.  These  small  variations  can  and  frequently  do  result  in  drastically  different 
performance  and  cause  imexpected  problems  during  testing. 

Numerous  systems  have  been  used  over  the  years  to  evaluate  model  /  sabot 
performance  during  launch  and  subsequent  flight.  These  include  muzzle  X-rays,  high 
speed  cameras,  radar,  witness  panels  and  others.  As  mentioned  the  primary  purpose  of  the 
evaluation  is  to  determine  model  and  sabot  integrity  during  launch,  sabot  separation,  and 
component  trajectories  after  launch.  X-rays  are  fi’equently  used  at  the  muzzle  to  view  the 
model  /  sabot  package  upon  exit  because  of  the  debris  cloud  ( smoke,  gasses,  etc.) 
usually  accompanying  the  package.  An  example  X-ray  is  shown  in  Fig.  4.1.  This  X-ray 
image  illustrates  an  excellent  launch,  both  model  and  sabot  appear  to  be  in  perfect 
condition  and  the  sabot  is  separating  as  designed.  Such  muzzle  X-rays  do  not  assist  in 
evaluating  how  much  separation  is  achieved  at  a  given  distance,  i.e.  blast  chamber  wall  or 
stripper,  or  what  the  trajectory  of  the  sabot  components  might  be.  One  simple  technique 
of  evaluating  these  concerns  is  the  placement  of  witness  panels  at  desired  distances. 

Figure  4.2  is  a  typical  example  of  the  information  that  can  be  obtained  from  these  witness 
panels.  This  figure  illustrates  that  the  sabot  petals  have  separated  about  14  inches  in  30 
feet  and  would  be  trapped  in  the  blast  chamber  of  the  facility;  thereby ,  eliminating  any 
potential  damage  to  the  instrumentation  systems  due  to  petal  impacts.  However,  other 
components  of  the  sabot  (  3  plugs  )  would  likely  make  it  through  the  stripper  creating  a 
possible  hazard  to  the  downrange  systems. 

Frequently  high  framing  rate  cameras  are  used  to  evaluate  the  performance  of  the 
model  /  sabot  package  during  the  initial  portion  of  the  flight.  Although  these  cameras 
cannot  view  the  package  immediately  upon  exit  from  the  muzzle,  due  to  the  debris  cloud, 
the  model  and  sabot  petals  can  be  viewed  a  few  feet  from  the  muzzle  and  this  separation 
evaluated.  This  technique  is  especially  useful  for  relatively  large  packages  (  several 
inches)  and  low  velocities  ( less  than  2000  ft  /  sec ).  Image  resolution  and  motion  blur  can 
be  a  problem  when  a  large  field  of  view  is  being  observed,  especially  for  high  velocity 
launches. 

4.2.2  Test  Conditions 

Another  purpose  of  these  preliminary  checkout  tests  is  to  determine  the  launch 
conditions  required  to  achieve  the  desired  test  velocity.  Although  the  corresponding 
launcher  conditions  are  normally  predicted  as  discussed  in  Section  3.4.1  it  is  desirable  to 
confirm  these  predictions  using  the  actual  model  /  sabot  package.  This  desirability  arises 
from  several  variables  which  can  significantly  affect  the  expected  muzzle  velocity  as 
determined  from  predictions  or  even  from  similar  previous  launches.  For  example 
variations  in  sabot  tolerances  can  drastically  alter  the  muzzle  velocity  because  the  start 
pressure  can  change  significantly.  The  start  pressure  is  the  pressure  on  the  base  of  the 
model  or  sabot  when  the  package  initially  begins  to  move  down  the  launch  tube.  If  the 
package  is  a  “loose  fit”  it  can  start  to  accelerate  too  early  in  the  laimch  cycle  and  exit  the 
muzzle  before  all  of  the  energy  in  the  propellant  is  expended.  Quite  often  imbumed 
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propellant  is  found  in  the  blast  chamber  when  this  occurs  and  the  resulting  muzzle 
velocity  is  lower  than  expected.  In  practice  designers  usually  attempt  to  tailor  the  start 
pressure  in  order  to  achieve  maximum  performance.  Other  variations  such  as  launch  tube 
wear,  propellant  conditions(  damp  or  dry ),  diaphragm  tolerances,  and  a  host  of  other 
causes  can  affect  the  expected  muzzle  velocity. 

This  velocity  determination  has  been  accomplished  on  the  outdoor  range  in  the  past 
using  several  different  techniques,  i.e.  circuit  break  /  make  paper,  multiple  pulse  X-rays, 
high  speed  framing  cameras,  digital  CCD  camera  systems,  and  Doppler  radar.  All  of 
these  techniques  have  their  advantages  and  depend  on  the  resources  available,  i.e.  funds, 
equipment,  and  man  power.  For  example  the  circuit  break  /  make  paper  technique  ( see 
Ref  4.2  )  is  a  basic,  simple,  inexpensive  system  but  is  manpower  intensive  and  is 
seriously  affected  by  weather  conditions  ( wind,  rain,  etc. )  on  the  outdoor  range. 

Whereas  the  Doppler  radar  unit  is  a  more  sophisticated  system  and  initially  more 
expensive.  This  system  is  also  more  versatile,  see  sample  data  in  Fig.  4.3,  and  has  a 
multitude  of  uses,  i.e.  muzzle  velocities,  in  bore  velocities,  spin  determination  and  long 
range  elevated  trajectory  measurements'* 


Time ,  sec 

Figure  4.3  Doppler  Velocity  Tracks  vs  Time  ( 25  mm  Projectile ) 


The  primary  advantage  of  the  photographic  techniques  ( X-ray  and  cameras  )  in 
obtaining  initial  velocities  is  that  the  designer  /  test  engineer  can  also  simultaneously 
evaluate  model  /  sabot  integrity  during  launch.  In  practice  even  when  the  radar  unit  is 
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being  utilized  it  is  usually  desirable  to  operate  one  or  more  of  these  photographic  systems 
for  this  evaluation.  Operating  both  types  of  systems  when  possible  provides  an 
independent  check  on  velocity  as  well  as  a  backup  in  case  one  system  malfunctions.  It 
can  also  permit  the  determination  of  the  initial  disturbances  transmitted  to  the  model 
during  launch.  This  is  often  an  important  reason  for  conducting  the  preliminary  outdoor 
tests  in  the  first  place.  The  model  /  sabot  designer  needs  to  insure  that  at  least  a  certain 
amount  of  angular  motion  will  be  present  for  the  trajectory  analysis  and  coefficient 
extraction  process  to  be  successful.  If  adequate  initial  disturbances  are  not  imparted  to  the 
model  then  augmentation  may  be  required  as  discussed  in  Section  3.4.3. 

Finally  useful  aerodynamic  data  can  also  be  obtained  during  these  outdoor  tests  if 
adequate  preparation  and  thought  are  expended.  In  fact,  organizations  that  do  not  have 
readily  available  access  to  an  indoor  ffee-flight  facility  can  accomplish  fairly 
sophisticated  analysis  on  an  outdoor  range  using  the  techniques  described  herein  and  in 

Refs.  4.4  and  4.5. 

4.2.3  Instrumentation 

Occasionally  it  is  necessary  to  do  preliminary  testing  in  order  to  checkout,  adjust,  set 
delays,  etc.,  for  special  instrumentation  systems  required  for  a  particular  tests.  Normally 
this  can  not  be  accomplished  on  the  outdoor  range  but  must  be  done  in  the  indoor  facility 
because  of  lighting,  protection  from  the  weather,  or  the  impracticality  of  moving  the 
systems  to  an  outdoor  environment.  For  these  occasions  the  general  practice  is  to  fire 
known,  well  behaved  items,  i.e.  20  or  30  mm  target  practice  (  TP)  ammunition,  through 
the  systems  for  instrumentation  evaluation.  Once  the  instrumentation  is  deemed  ready 
then  the  correct  launcher  is  installed  and  the  actual  test  items  laxmched.  Using  the  well 
behaved  TP  ammunition  not  only  minimizes  the  risk  of  damage  to  the  instrumentation 
systems  but  saves  the  test  items  for  launching  when  the  chances  of  success  are  the 
highest. 

The  process  discussed  above  is  also  followed  after  periods  of  downtime  in  order  to 
ensure  that  the  basic  shadowgraph  system  is  in  a  state  of  readiness.  After  periods  of 
extended  downtime,  due  to  calibration  or  holidays,  camera  shutters  tend  to  stick,  spark 
sources  can  collect  moisture  and  become  erratic,  etc.  A  few  checkout  flights  using  well 
behaved  and  inexpensive  projectiles  can  be  useful  in  isolating  these  problems  and  getting 
the  facility  back  in  an  operational  status. 

4.3  TESTS 

The  actual  testing  can  be  the  shortest  and  most  straight  forward  phase  of  the  overall 
program  as  compared  to  the  pretest  and  preparation  or  the  post  test  analysis.  However, 
this  can  be  the  phase  which  receives  the  most  visibility .  Often  user  representatives, 
management  personnel,  and  contractor  employees  are  present  and  their  impressions  of 
what  they  observe  during  this  time  forms  the  basis  of  their  lasting  memories  of  the 
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program.  A  deliberate  professionally  executed  program  which  follows  established 
procedures  will  provide  a  positive  image  and  result  in  a  safe  and  successful  test  program. 

4.3.1  Test  Plan 

As  a  result  of  the  pretest  planning  discussed  in  the  previous  section  and  the 
preliminary  tests  discussed  above  the  test  engineer  is  in  a  position  to  provide  the 
operations  crew  with  a  test  plan.  This  plan  can  be  informal,  but  should  include  the 
information  required  by  the  crew  to  effectively  and  efficiently  conduct  the  test.  For 
example  it  should  include: 

1 .  proposed  test  matrix,  i.e. 

a.  expected  velocities,  number  of  shots 

b.  launcher  information,  propellant,  charge  conditions,  etc. 

2.  facility  configuration,  i.e. 

a.  portable  catcher  requirement 

b.  freon  tube  requirement,  charge  weight 

c.  sabot  stripper 

3.  instrumentation  requirements,  i.e. 

a.  direct  shadowgraph 

b.  interferometry 

c.  radar 

4.  any  other  information  that  is  pertinent  to  the  test  program. 


4.3.2  Documentation  and  Data  Collection 

The  test  engineer  should  ensure  that  the  activities  are  adequately  documented  during 
the  test  process.  Although  much  of  this  documentation  has  been  computerized  over  the 
years  personalized  notes  and  observations  can  also  be  extremely  valuable  later  when  the 
trajectories  are  being  analyzed.  Also  the  information  being  recorded  and  saved  in  the 
computerized  shot  log  must  be  checked  for  accuracy  and  completeness.  An  example  of  a 
daily  computer  log  is  shown  as  Fig.  4.4.  As  can  be  seen,  shot  date  and  time,  facility 
atmospheric  conditions,  model  measurements,  launcher  description,  propellant  type,  and 
velocity  achieved  are  all  recorded.  Some  of  this  information,  atmospheric  conditions  and 
model  measurements,  are  passed  to  the  image  processing  routine  and  subsequently  used 
in  the  trajectory  calculations.  Immediately  after  the  shot  the  hall  and  pit  time  information 
is  recorded  and  also  passed  to  the  film  reader  routine.  A  sample  print  out  of  the  time 
information  is  shown  as  Table  4.1  and  will  be  discussed  in  more  detail  in  Section  V.  The 
information  in  Fig.  4.4  and  Table  4.1  and  any  notes  the  test  engineer  may  have  recorded 
along  with  all  of  the  photographs,  i.e.  shadowgrams,  direct  photographs,  holograms,  etc., 
obtained  represents  the  primary  shot  documentation  during  and  immediately  after  the  test. 
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♦•♦int#*###-*************#****#***  ARF  SHOT  LOS  ******************************** 
SHOTAS  95073104  DATE;  7/31/95  TIME:  9:4(:i 

TITLE:  WAF/ISL 

REMARKS; 

*##**»#******##**•****»#♦****■*•#***•»•****♦**•*•**■*•****##******♦■*•♦***♦**♦**♦***■*■«■♦**** 
TEMPS  (deg  F) :  72.00  72.00  71.00  71.00  71.00  70.00  69.00  71.00 

AIR  PRESSURE;  1017.27  mb  TEMP;  21.60  deg  C  REL  HUMIDITY;  53.00  7. 

DENSITY;  1.201873  kg/m-'S  VISC;  1.822xl0''-5  N*s/ni''2 
**•*■*****♦*♦***•***■*•*♦*******♦**•»**♦*******»***#*♦**********■*•■*•********♦*********** 
MODEL  NUM;  4 

MODEL  MASS;  86.000  gm  SABOT  MASS;  490.500  gm  PACKAGE  MASS;  576.500  gm 

CG  FROM  NOSE;  59.120  mm  MODEL  LENGTH;  114.355  mm  MODEL  DIA:  18.954  mm 

Ix  ;  45.9998  gm*cm''2  ly;  879.4477  gm*cm''2  Iz:  879.4477  gm»cm-'2 

Ixy;  .  gffl*cm'^2 

MODEL  DESCRIPTION:  ALUMINIUM  W/STEEL  NOSE  THICK  WAF 

SABOT  DESCRIPTION:  4  PIECE  LEXAN  W/  2  PIECE  POLY/STEEL  PUSHER 

*#***##*.**##*♦*****#**#♦*****#**♦***#*♦♦*♦**#**♦#*##***♦***************#*******•* 

GUN  TYPE:  2. SIN  *  POWDER  MASS:  400.00  gm  POWDER  TYPE;  GAU-3 

2AS  PRES;  .  psi  PISTON  WT:  .lb  BURST  PR;  .  psi 

EXPECTED  V:  2200.00  ft/s  FREON;  0.00  gm  TWIST:  .  cal /rev 

+♦  RESULTS  **#»♦♦**■*•«****♦♦*******♦*■**»*♦*♦*♦**********♦♦•*■**«•♦*•*•**♦*****•*•*♦**« 
V  MUZZLE:  2063.62  -ft/s  V  MID;  1575.19  ft/s  V  IMPACT;  1193.05  ft/« 

mCH:  1.39  CD  MID:  1.253 

K  (DIA):  0.6003x10^-6  RE  (LTH) :  3.6217x  10-'6 

Figure  4.4  Example  of  Daily  Shot  Log  Entry 

All  of  the  photographic  results  should  be  examined  carefully  by  the  test  engineer  to 
ensure  that  no  damage  occurred  to  the  test  item  during  the  laimching  process.  Generally 
speaking  models  that  have  been  damaged  in  some  way  should  be  deleted  from  the 
analysis  process.  Retaining  results  obtained  from  damaged  models  could  contaminate  the 
final  conclusions  associated  with  the  actual  configuration  of  interest,  ( the  undamaged 
models).  Some  test  items  are  purposely  altered  during  the  laimching  process  and  the 
resulting  data  from  such  configurations  are  meaningful.  For  example,  spin  stabilized 
projectiles  are  designed  to  engage  the  rifling  of  the  launch  tube.  Hence  the  rifling  in  the 
barrel  cuts  grooves  into  the  projectile’s  rotating  bands  and  frequently  scores  the 
cylindrical  section  of  the  projectile.  Since  this  ‘damage’  is  typical  of  the  operational 
ammunition  the  data  resulting  from  such  trajectories  are  of  importance.  However,  such 
variations  from  model  to  model  may  not  be  totally  consistent  and  the  final  aerodynamic 
coefficients  and  derivatives  may  exhibit  more  apparent  scatter  than  would  have  otherwise 
occurred  for  a  saboted  model. 

Another  problem  the  test  engineer  should  be  alert  to  is  possible  reflected  shock 
waves  off  the  walls  of  the  facility.  These  can  be  evident  for  models  which  decelerate 
through  Mach  1  during  their  flight.  At  this  condition  the  near  normal  bow  shock  can 
reflect  off  the  walls  of  the  facility  and  impinge  on  the  model.  This  reflected  shock  can 
then  add  an  impulse  to  the  model’s  motion.  Such  trajectories  should  then  be  treated  as 
two  separate  flights.  The  first  limited  to  that  portion  of  the  flight  prior  to  the  impingement 
of  the  reflected  shock  and  the  second  after  the  impingement  of  the  reflected  shock.  As 
will  be  discussed  in  Section  V  the  shadowgraph  images  are  digitized  and  the 
experimentally  measured  trajectories  computed  using  the  CADRA'*  ’  system. 
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Table  4.1 


Sample  Hall  and  Pit  Time  Printout 


3HGT  NUMBER:  91 

AEROBALLISTIC  RESEARCH  FACILITY^  , 

5  650  TIME!  12:17  DATE:  6/  5/91  VMUZ:  7U49.76  VMID:  4042.51  VIMP 

MGDEL  DESCRIPTION;  HART 

MODEL  13  850SRS 

3331  4Kdg.  400  psi  HEM.  LLG 

STATION 

DELAY 

HALL 

PIT 

F 

RAW 

F 

SMOOTHED 

DRAG 

NUMBER 

TIME 

TIME 

TIME 

1 

VELOCITY 

VELOCITY 

COEF 

1 

0.01 

247294.6 

247293.5 

* 

0.  00 

* 

6707. 16 

4.888 

2 

0.01 

248130.3 

248190.7 

* 

5770.65 

* 

6670.49 

4.915 

T 

0.01 

248965.3 

248967.0 

* 

6206.18 

* 

6633.82 

4.942 

4 

0.01 

250609.7 

250608. 0 

* 

6075. 35 

* 

6560.62 

4.997 

5 

0.01 

253944.7 

253948.6 

* 

5991.98 

* 

6413.93 

5.112 

6 

0.01 

255616.6 

255631.2 

* 

5962. 14 

* 

6340,59 

5.  171 

7 

0.01 

257274.0 

257273.5 

6012.66 

* 

6267.83 

5.231 

3 

0.21 

258103.8 

258110.7 

* 

6375.60 

* 

6228.73 

5.264 

9 

0.01 

258946.2 

258947.5 

* 

5769.90 

* 

6193.31 

5.294 

iO 

0.01 

262288.6 

262289.7 

* 

6007.86 

6046.03 

5.423 

11 

0.01 

265674.5 

265690.0 

* 

5967.98 

* 

5897.51 

5.559 

12 

0.01 

267336.5 

267335.4 

* 

5944.24 

* 

5825.41 

5.628 

13 

0.01 

269054.8 

269031.8 

* 

5903.85 

* 

5751.48 

5.700 

14 

0.01 

270741.1 

270741.8 

5841.68 

* 

5678.72 

5.773 

15 

0.01 

271575.4 

271605.7 

* 

5888.63 

* 

5642.05 

5.811 

16 

0.01 

275017.8 

275017.8 

5858.90 

* 

5494.77 

5.967 

17 

0.01 

276755.6 

276756.6 

5850.61 

* 

5420. IS 

6.049 

18 

0.01 

277639.0 

277637.6 

* 

5622.28 

* 

5333. 80 

6.090 

19 

0.01 

232186.5 

282185.8 

5459.71 

* 

5201.69 

6 . 303 

20 

0.01 

283969.3 

283966.5 

* 

5567.56 

5128.93 

6.392 

21 

0.01 

285388.2 

285906.2 

* 

5141.82 

* 

5056. 17 

6.484 

oo 

0.01 

287733.7 

287823.2 

* 

5197. 12 

* 

4983.41 

6.579 

23 

0.01 

289710.5 

289707.1 

* 

5344.59 

* 

4908.82 

6.679 

24 

0.01 

291711.5 

291711.3 

4908.60 

* 

4836.72 

6.778 

25 

0.01 

293765.9 

293764.8 

* 

4946.55 

* 

4762.20 

6.884 

26 

0.01 

294795.6 

294795.4 

5009. 03 

* 

4724.36 

6.940 

27 

0.01 

297984.2 

297980.2 

* 

4653.77 

4615.59 

7.103 

28 

0.01 

301312.2 

301314.2 

* 

4503.12 

4505.57 

7.277 

29 

0.01 

304699.0 

304700.4 
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Note:  Stations  at  the  far  end  of  the  range  will  take  up  to  3.3  micro  sec's 
longer  to  receive  the  common  start  pulse.  This  time  delay  has  been 
measured  at  each  station  and  is  compensated  for  in  the  data  reduction. 
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When  tests  are  accomplished  in/on  a  facility  which  does  not  use  a  daily  shot  log  of 
some  type  the  test  engineer  is  primarily  responsible  for  documenting  the  process  and 
maintaining  an  adequate  record  of  the  test.  This  can  be  the  situation  for  some  outdoor  test 
environments;  but,  even  here  a  standardized  record  should  be  maintained  covering  the 
equivalent  information  that  is  shown  in  Fig.  4.4.  The  importance  of  documentation  can 
not  be  over  stated  because  a  lack  of  an  adequate  record  can  result  in  useless  shots 
requiring  repeat  firings  wasting  resources  and  manpower. 
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V.  TRAJECTORY  DETERMINATION 


5.1  IMAGE  ANALYSIS 

As  the  model  travels  downrange  it  passes  the  photographic  stations  and  orthogonal 
images  are  obtained  at  each  station  as  was  illustrated  in  Fig.  1.2.  Each  camera  records  the 
shadow  of  the  model  on  the  reflective  screen  as  cast  by  the  light  from  the  associated 
spark  units.  Each  image  also  includes  the  facility’s  reference  wires  and  beads  as  shown  in 
the  typical  shadowgraph  of  Fig.  5.1.  At  the  present  time  film  cameras  are  used  as  the 
primary  recording  medium  and  the  time  of  flight,  shot  and  station  identification  are  also 
included  on  the  film  format  as  shown  in  Fig.  5. 1 .  Recently  the  microdensitometer 
previously  used  to  read  the  film,  Ref  5. 1 ,  has  been  replaced  with  the  film  scanner  and 
image  processing  system  shown  in  Fig.  5.2.  Using  this  system  the  film  images  are 
scanned  and  the  points  of  interest  on  each  image  are  numerically  coded  and  the  range 
coordinates  computed.  Once  the  film  images  are  digitized  using  the  Ektron  model  1412 
digital  imaging  camera  the  digital  image  is  then  processed  so  that  only  the  area  of  interest 
(projectile  and  reference  beads)  need  be  saved.  The  image  processing  and  storage  are  then 
accomplished  using  a  Macintosh  Quadra  800  computer  and  peripherals  with  the  CADRA 
interactive  software  system.  The  film  scanning  and  CADRA  system  is  described  in  detail 
in  Ref.  5.2.  Once  the  images  are  scanned  and  saved  they  may  be  retrieved  as  required  for 
the  actual  film  reading. 

It  is  planned  to  replace  the  film  cameras  with  solid  state  CCD  cameras  in  the  near 
future  ’  .  This  plan  calls  for  cameras  with  Ik  x  Ik  pixel  arrays  in  the  up-range  end  of 
the  ARF  and  2k  x  2k  pixel  arrays  in  the  downrange  portion  of  the  facility.  It  is  expected 
that  some  of  these  cameras  will  be  installed  during  1997  but  that  it  will  require  several 
years  to  completely  outfit  the  facility.  Once  this  upgrade  has  been  accomplished  the 
digitized  images  will  be  transferred  directly  to  the  CADRA  system  for  analysis  and  the 
film  processing  and  film  scanning  will  be  bypassed.  At  present  the  film  processing 
generally  requires  waiting  overnight  to  see  the  film  from  the  previous  day’s  testing. 
Whereas  the  film  scanning  requires  about  2  hours  for  a  full  range  flight.  This  upgrade  will 
therefore  enable  the  test  engineer  to  have  access  to  the  experimentally  measured 
trajectories  for  analysis  shortly  after  the  launch,  certainly  on  the  same  day. 

Although  the  pixel  arrays  do  not  provide  equal  image  resolution  as  film,  the  technique 
used  in  processing  the  electronically  stored  images  utilizes  the  complete  image  in  contrast 
to  the  previous  film  method  where  only  specific  points  on  the  model’s  shadow  were 
determined,  i.e.,  nose,  comers  of  the  base,  and  spin  pin.  Both  techniques  require  locating 
the  reference  beads  since  these  represent  the  facility’s  axis  system.  As  a  result  of  utilizing 
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the  complete  image  the  precision  to  which  the  model’s  positions  and  attitudes  are 
determined  has  actually  increased  even  though  the  image  resolution  has  decreased.  This 
technique  is  illustrated  using  the  following  idealized  sketch  of  a  typical  shadowgram. 


Figure  5.3  Idealized  Sketch  of  a  Typical  Shadowgram 
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The  equations  for  the  centroid  of  the  shadow  are  given  by: 


N 

IxfndA  i;  xfnAA 

_  0 

— —  N  ' 

Jf„dA  p„AA 


and  N 

J  y  fn  d  A  S  y  fn  AA 

Cy  ~  “  - N 

JfndA  ZfnAA 


> 


5.1 


Where: 

AA  =  area  of  one  pixel 

fn  =  fraction  of  each  pixel  area  filled  by  the  model  shadow 
N  =  number  of  pixels  in  the  shadow 


Once  the  centroid  of  the  shadow  has  been  determined  the  angle  X  can  be  computed  by 
finding  the  axis  for  which  the  second  moment  of  inertia  of  the  shadow  area  is  a  minimum. 


Or, 


tan2X,  =  -2  Ixys  /  ( Ixs  -  lys ) 


5.2 


where, 


^ys 


J  x'^  dA 


N 

Zx'^  AA 


^XS 


=  jy'^dA 


N 

=  ly 

1 


a 


AA 


N 

Ixys  =  J  y'  x'  dA  =  Z  y'  x'  AA 


All  that  remains  is  to  calculate  the  x,y  coordinates  of  the  nose  and  the  center  of  the 
base.  Both  of  these  points  lie  on  the  line  defined  by  the  centroid  and  the  angle  X  and  can 
be  located  by  scanning  along  this  line  vmtil  the  gray  scale  levels  change  rapidly  indicating 
the  limits  of  the  projectile  shadow.  A  more  complete  description  of  this  technique  can  be 
found  in  Ref.  5.3. 
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5.2  POSITION  AND  ATTITUDE  MEASUREMENTS 


In  order  to  calculate  the  model’s  position  and  attitude  in  space  it  is  necessary  to 
convert  the  planar  measurements  obtained  from  the  images  and  discussed  in  the  previous 
section  to  the  three-dimensional  geometry  of  the  facility.  This  process  will  be  described 
for  a  typical  shadowgraph  station  as  illustrated  in  Fig.  5.4.  In  this  figure  it  should  be 
recognized  that  the  reference  bead  coordinates  ( Xgi ,  Ybi  ,  Zbi  )  and  the  spark  gap 
locations  ( Xj ,  Yj ,  Zj )  where  i  =  H,  P  referring  to  the  hall  and  pit  planes  respectively  are 
known  from  calibrating  the  range  as  was  discussed  in  Section  3.5.3  and  Appendix  2. 


-  -  -  ^ 

Hall  Spark 


(Xh.Yh.Zh) 


Figure  5.4  Projectile  Attitude  /  Positioning  Scheme 


The  three  dimensional  range  coordinates  of  the  model’s  shadow  can  be  related  to  the 
image  measurements  of  the  previous  section  using  the  relationships  shown  in  Table  5.1. 
Note  that  the  hall  screens  are  approximately  0.25  inches  behind  the  beads  on  the  reference 
system  wires.  Since  the  reference  wires  define  the  zero  Y  coordinate  the  range  Y 
coordinate  for  all  shadows  on  the  hall  screen  are  -0.25  inch.  Since  the  overhead  screen  is 
about  0.25  inch  above  the  beads  the  Z  coordinates  of  all  the  overhead  shadows  are  Zj  + 
0.25. 
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Table  5.1 

Relationship  Between  Film  /  Image  and  Range  Coordinates 


Point 

Hall  Image 

Pit  Image 

Film  /  Image 
Coordinarte 

Range 

Coordinate 

nim  /  Image 
Coordinate 

Range 

Coordinate 

Reference  Bead 

^BH 

^B 

^P 

^BH 

^B 

"'^BP 

^BP 

Nose 

^NS 

"nh 

X 

NS 

X 

NP 

'^NS 

^NH 

^NS 

^NP 

Y^H  =  -0.25 

Znp+0.25 

Center  of  Base 

\s 

X 

TS 

X 

TP 

^TS 

^TH 

Y 

TS 

'^TP 

Y  =  -0.25 

TH 

Z.pp+0.25 

Spin  Pin 

^PS 

^PH 

PS 

^PP 

"^PS 

^PH 

Y 

PS 

^PP 

VpH  =-0-25 

Zpp+0.25 

Note  that  the  hall  screens  are  approximately  0.25  inchs  behind  the  beads  on  the  reference 
system  wires.  Since  the  reference  wires  define  the  zero  y  coordinate  the  range  y 
coordinate  for  all  shadows  on  the  hall  screens  are  -0.25  inch.  Also  note  that  the  overhead 
screens  are  approximately  0.25  inch  above  the  overhead  reference  wire. 

5.2.1  Direction  Cosines 

The  mathematical  process  of  determining  the  direction  cosines  (me,  ng,  Pe )  of  the 
model’s  principal  axis  becomes  a  problem  of  determining  the  intersection  of  two  planes. 
These  two  planes  are  defined  as  shown  in  Fig.  5.4  or,  plane  1  is  specified  by  the  hall 
spark  location  and  the  centerline  of  the  model’s  shadow  on  the  wall  reflective  screen  ( H 
screen  in  Fig.  5.4).  Plane  2  is  specified  by  the  pit  spark  location  and  the  centerline  of  the 
model’s  shadow  on  the  ceiling  reflective  screen  ( P  screen  in  Fig  5.4 ).  It  is  assumed  that 
the  centerline  of  the  model’s  shadow  represents  the  projection  of  the  model’s  principal 
axis  on  the  reflective  screens.  For  asymmetric  configurations  or  even  symmetric  shapes 
that  may  project  a  shadow  which  is  asymmetric,  i.e.  a  missile  with  three  fins,  a  template 
matching  technique  has  been  developed  and  is  discussed  in  Ref  5.2.  Essentially  that 
technique  computes  how  the  shadow  of  the  model  should  appear  on  the  reflective  screens 
at  various  pitch,  yaw  and  roll  angles  until  the  computed  shadow  matches  the  images 
obtained  from  the  shadowgrams.  However,  for  the  illustrative  discussion  herein  we  will 
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assume  a  symmetrical  shadow  projection.  It  should  also  be  noted  that  there  are  other 
methods  of  determining  the  direction  cosines  of  the  model’s  principal  axis  but  the  method 
of  intersecting  planes  utilizes  the  model’s  complete  shadow  in  identifying  the  angle  X  as 
defined  in  Eq.  5.2. 


The  general  equation  for  a  plane  can  be  written  as, 


PX  +  QY  +  RZ  =  1 


5.3 


where  P,  Q,  R  are  the  direction  numbers  of  a  line  normal  to  the  plane.  If  we  apply  Eq.5.3 
to  the  plane  associated  with  the  horizontal  plane  ( H  screen )  then  using  the  nomenclature 
in  Fig.  5.4  we  can  write  the  three  relations  below: 

Ph  ^nh  +  Qh  Ynh  +  Rh  ^nh  “  1 
P XjH  +  Qh  Yth  +  Rh  ^th  ”1  5.4 

Ph  ^H  +  Qh  Yh  +  Rh  Zh  =  1 


Equations  5.4  may  be  solved  for  P^,  Qh,  and  Rh  and  written  in  matrix  form  as: 


-1 


■Ph" 

Xnh 

Ynh 

Znh 

~r 

Qh 

= 

X-PH 

Yth 

Zth 

1 

5.5 

LrJ 

LXh 

Yh 

Zh  _ 

i_ 

If  we  apply  Eq.  5.3  to  the  vertical  plane  ( P  screen )  then  we  can  similarly  write: 


-1 


~Pp" 

^NP  Y^P  ^NP 

"1 

"I" 

Qp 

= 

X-pp  Y-pp  Zjp 

1 

5.6 

LRpJ 

.Xp  Yp  Zp  _ 

_  1_ 

and  solve  for  Pp,  Qp,  and  Rp. 

Once  the  two  planes  are  defined  we  need  to  find  the  line  that  is  common  to  both 
planes.  This  line  represents  the  principal  axis  of  the  model.  To  accomplish  this  we  will 
define  two  points  that  lie  on  this  line  as  (  X,,  Yj,  Zj  )  and  ( X2,  Y2,  Z2  ).  Then  using  Eq. 
5.4  we  can  further  write: 

QhYi  +  RhZi  =  1  -  PhXi 

Qp  Y]  +  Rp  Z]  =  1  -  PpX|  5.7 

QhY2  +  RnZa  =  1  -  PhX2 

QpY2  +  RpZ2  =1 -PpX2  5.8 

If  we  solve  the  two  5.7  equations  for  Yj ,  set  equal  and  then  solve  for  Zj  we  can  obtain 
after  some  manipulation, 

Z,=  [(1-PpX,)Qh-(1-PhX,)Qp]/(QhRp-QpRh)  5.9 
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Now  if  we  solve  the  two  5.7  equations  for  Z,  ,  set  equal  and  solve  for  Yj  we  can  similarly 
arrive  at, 

Y,  =  [( 1  -  PhX,  )Rp  -  ( 1  -  PpX,  )Rh  ]  /  (  QhRp  -  QpRh  )  5.10 

Duplicating  these  derivations  for  Eqs.  5.8  we  get, 

Z2=  [(1-PpX2)Qh-(1-PhX2)Qp]/(QhRp-QpRh)  5.11 
and 

Y2  =  [( 1  -  PhX2  )Rp  -  ( 1  -  PpX2  )Rh  ]  /  (  QhRp  -  QpRh  )  5.12 

Substituting  Eqs.  5.9  through  5.12  into  the  definitions  of  the  direction  numbers,  m„, 
Un,  and  pn  ( Eqs.  5.13  ), 

m„  =  Y2-Y, 
nn=  Z2-Z1 

p„  =  X2-X, 

one  arrives  at,  after  some  manipulation,  the  relations  for  the  direction  numbers  shown 
below. 

m„  =  -  pn  [(  RpPh  -  RhPp  )  /  (  QhRp  -  QpRh  )] 

Hn  =  -  Pn  [(  QnPp  ■  QpPn  )  /  ( QnRp "  QpRh  )] 


If  we  selectively  choose  the  two  points  lying  on  this  line  of  intersection  such  that 
X2  -  Xi  =  1,  or  Pn  =1 ,  then  the  direction  cosines  of  the  models  principal  axis  is, 

m.  =  mn  /  D,  ) 

ne=n„/De  5.14 

Pe=  1  /De  j 

where 

/  1  I  2  ,  2  xl/2 

DeS(l+mn  +nn  ) 


Equations  5.14  are  used  to  compute  the  direction  cosines  of  the  model’s  principal 
axis  at  each  shadowgraph  station  after  the  model  has  traversed  the  instrumented  range.  It 
should  be  recognized  that  by  selectively  choosing  the  two  points  such  that  X2  -  Xj  =  1  we 
have  assumed  that  the  model’s  angle  of  attack  will  never  be  greater  than  90  degrees  or 
that  the  model  can  not  tumble.  This  assumption  has  never  been  a  problem  in  the  past. 
However,  if  in  the  future  a  model  is  launched  ‘backwards’  in  order  to  obtain  very  high 
angle  of  attack  data  then  this  assumption  might  have  to  be  removed  from  the  derivation. 
The  geometric  principles  used  in  this  and  the  following  derivations  have  been  extracted 
from  various  text  books,  see  for  example  Ref  5.5. 
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5.2.2  Position 


Now  that  we  have  the  projectile’s  orientation  defined  we  need  to  determine  the  range 
coordinates  of  the  projectile.  For  discussion  purposes  we  will  select  the  nose  tip  of  the 
projectile  as  the  point  that  we  will  determine.  However,  it  should  be  imderstood  that 
the  method  described  can  be  used  to  determine  the  range  coordinates  of  any  point  on  the 
projectile  including  the  spin  pin  or  the  centroid  of  the  shadow  as  defined  in  Section  5.1. 
The  range  coordinates  of  the  projectile’s  nose  ( X^,  Y^,  )  are  defined  as  the  nearest 

point  of  intersection  of  the  two  lines  connecting  the  hall  spark  coordinates  (  Xh,  Yh,  Zh  ) 
and  the  shadow  of  the  nose  on  the  hall  screen  (  X^hs  Ynh»  Z^h  )  and  the  pit  spark  location 
(Xp,  Yp,  Zp  )and  the  shadow  of  the  nose  on  the  overhead  screen  (  X^p,  Y^p,  Z^p), 
see  Fig.  5.4.  Recognize  that  there  is  a  point  on  each  of  these  lines,  ( Xj  ,Yi ,  Zj )  where 
i  =  1,2  respectively,  such  that  the  distance  between  these  two  lines  is  a  minimum. 


Using  the  general  equation  for  a  straight  line,  the  following  relationships  can  be 
written: 


and 


x.-x.  Y,-Yh_  Z^-Z, 

y  ^  =  V  -  V  “  7.7 

NH  H  'nH  H  H 


Xp-Xp  ^  Y,-Yp  Z,-Zp 

^NP-Xp-  Y^p-  Z^p-Z, 


The  distance  between  points  1  and  2  is  defined  as. 


5.15 


5.16 


d  =  [(X2 

-X,f  +  (Yj-Y,)"  +  (Z2-Z,)y“ 

5.17 

If  we  let. 

Xnh  ■  Xh  =  Aj 

>  Ynh-Yh  =  Bi  ,  Znh-Zh  =  Ci 

and. 

Xnp  -  Xp  =  A2 

,  Ynp  -  Yp  =  B2  ,  Znp  -  Zp  =€2 

Then  Eqs.  5.15  and  5.16  can  be  rewritten  as  shown  below: 


X.-Xh  Y,  -Yh  ^  Z^-Zh 


B 


5.18 


and 


Xp-Xp  _  Y,-Yp 


B. 


5.19 


Solving  Eqs.  5.18  for  Yj  yields, 

(X,-Xh)B,  =  (Y,-Yh)A,  =i>  Y.  =  (B,/A.)Xi-D,/A. 


5.20 
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where, 


D]  =Xh  B]  -  Yh  Aj 


Similarly  solving  for  Y2  ,  Eq.  5.19,  one  arrives  at: 

Y2  =  ( B2  /  A2)  X2  -  D2/  A2  ;  where,  D2  =  Xp  Bj  -  Yp  A2 

Subtracting  Eq.  5.20  from  5.21  we  then  get, 

Y2  -  Y,  =  ( B2/  A2)  X2  -  ( Bj/  A,)  Xi  +  K  ;  where,  K  =  Dj/  A,  -  D2/  A2 


5.21 


5.22 


Following  the  same  logic  as  outlined  above,  relations  for  Zj,  Z2  and  Z2  -  Zj 

can  also  be  derived  from  Eqs.  5.18  and  5.19.  The  resulting  equations  are  listed  below: 

Z,  =  (C,/  Ai)  Xi  -  E,/  A,  ;  where,  Ei  =  Xh  C,  -  Zh  A,  5.23 

Z2  =  (C2/ A2)  X2  -  E2/ A2  ;  where,  E2  =  XpC2-ZpA2  5.24 

Z2  -  Z,  =  ( C2/  A2)  X2  -  ( C,/  Ai)  X,  +  L  ;  where,  L  s  Ej/  Aj  -  E2/  A2  5.25 

Substituting  Eqs.  5.22  and  5.25  into  Eq.  5.17  the  distance  d  can  be  expressed  as 
d  =  {(X2  -  X,f  +  [(B2/A2)X2  -  (B,/A,)X,  +  K]'  +  [(C2/A2)X2  -  (C,/A,)X,  +  5.26 

Now  if  we  expand  the  terms  in  Eq.  5.26,  collect  like  terms,  take  the  partial  derivatives 
of  d  with  respect  to  both  Xi  and  X2,  set  these  partial  derivatives  equal  to  zero,  we  can 
arrive  at  the  two  relations  shown  below: 

HX, -1X2-0  =  0  5.27 

MX2- JXi  +  S  =  0  5.28 

where, 

H=l+(B,/Aif +  (C,/Aif 


J  =  1  +  (  B2B,)/(A2A,)  +  (C2C,)/(A2Ai) 
G  =  K(B,/A,)  +  L(Ci/A,) 

M  =  1  +  (B2/A2f +  (C2/A2)^ 

S  =  K(B2/A2)  +  L(C2/A2) 


Solving  Eqs.  5.27  and  5.28  for  X2  and  X,  we  arrive  at: 

X2  =  (  SH  -  JG )  /  ( -  MH ) 
Xi  =  ( JX2  +  G  )  /  H 


5.29 

5.30 


Once  we  have  determined  X2  and  Xi  we  can  substitute  these  back  into  Eqs.  5.20,  5.21  and 
Eqs.  5.23,  5.24  to  compute  Y,,  Y2  and  Z,,  Z2  respectively.  These  are  the  coordinates  of 
the  two  points  on  lines  1  and  2  which  represent  the  closest  point  of  intersection.  We  now 
assume  that  the  range  coordinates  of  the  nose  are  halfway  between  the  two  points.  Or 

Xn  =  (X2-X,)/2  ) 

Yn  =  (Y2-Y,)/2  \  5.31 

Zn  =  (Z2-Z.)/2  J 

We  can  now  compute  the  range  coordinates  of  the  model’s  nose,  or  any  other  point  on  the 
body  of  the  projectile,  in  terms  of  the  measured  shadow  coordinates  and  the  calibrated 
spark  gap  and  reference  bead  coordinates. 
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5.2.3  Roll  Orientation 


When  measuring  the  roll  orientation  of  a  projectile  in  range  coordinates  it  is  usually 
necessary  to  modify  the  model  in  some  manner.  Normally  this  modification  consists  of 
inserting  a  reference  point,  i.e.  spin  pin,  into  the  base  of  the  projectile.  This  spin  pin 
provides  a  radial  reference  location  from  which  the  model’s  roll  orientation  can  be 
determined  as  will  be  demonstrated  in  this  section.  In  order  to  minimize  the  roll 
measurement  error  the  spin  pin  should  be  positioned  radially  as  far  from  the  center  of  the 
base  (or  model’s  centerline)  as  practical.  Trailing  edge  fin  tabs  are  sometimes  used  in 
place  of  spin  pins  on  finned  models  because  of  this  requirement.  The  error  in  the  roll 
orientation  measurements  as  a  function  of  the  distance  (L)  between  the  model’s  centerline 
and  the  spin  pin  is  illustrated  in  Fig.  5.5. 


Figure  5.5  Accuracy  of  Roll  Orientation  Measurements 


The  curve  shown  in  Fig.  5.5  represents  a  worst  case  scenario  since  it  assumes  a 
measurement  error  ( )  of  0.05  cm  in  both  the  location  of  the  spin  pin  and  the  center  of 
the  base.  The  centerline  of  the  model  is  computed  using  the  complete  shadow  image  as 
discussed  in  Section  5.2;  hence,  the  precision  to  which  the  range  coordinates  of  the  center 
of  the  base  is  determined  is  significantly  better  than  that  of  locating  a  point.  Nevertheless, 
the  curve  does  demonstrate  the  importance  of  positioning  the  spin  pin  as  far  from  the 
model’s  centerline  as  practical. 

It  should  be  noted  that,  in  the  past,  frequently  two  spin  pins  were  used  in  determining 
the  roll  orientation.  Essentially  the  advantage  of  using  two  spin  pins  was  to  double  L  in 
the  figure  above.  The  disadvantage  of  using  two  spin  pins  was  in  identifying  which  pin 
was  which.  In  order  to  expedite  that  identification,  spin  pins  of  differing  diameters  were 
employed.  However,  because  of  flow  turbulence  in  the  model’s  wake,  trailing  shock 
waves,  etc.  it  was  often  difficult  to  correctly  differentiate  between  the  spin  pins  and  large 
errors  in  the  computed  roll  orientations  resulted.  Since  the  incorporation  of  CADRA^  ^  as 
the  image  processing  system  only  one  spin  pin  has  been  utilized  and  the  technique  of 


determining  the  roll  orientation  using  only  one  spin  pin  will  be  described  herein.  The 
process  would  be  the  same  if  two  spin  pins  were  used  except  the  range  coordinates  of  the 
center  of  the  base  would  be  replaced  with  the  coordinates  of  the  second  pin. 

The  range  coordinates  of  the  spin  pin  and  the  center  of  the  base  are  assumed  to  have 
been  determined  using  the  methods  described  in  the  previous  sections.  The  coordinates 
of  these  two  points  are  therefore: 

(Xp,Yp,Zp)  =  Spin  Pin 
( Xf ,  Yf ,  Zj)  =  Center  of  Base(Tail) 

The  direction  cosines  of  the  vector  connecting  the  spin  pin  with  the  center  of  the  base  are 
defined  as  shown  below: 

mp  =  (Yp  -  Yt)  /  Dp  =  cos  p 
Op  =  ( Zp  -  Zj)  /  Dp  =  cos  Y 
Pp  =  ( Xp  -  Xt)  /  Dp  =  cos  p 

Where  the  distance  between  the  pin  and  the  centerline  of  the  body  is, 

Dp  =  V  ( Xp  -  Xt)^  +  (  Yp  -  Yt)^  +  (  Zp  -  Zjf 
and  the  angles  P,  p,  and  y  are  defined  on  the  sketch  below. 
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If  we  assume  that  the  length  of  the  D  vector  defined  by  the  centerline  of  the  body  and 
the  tip  of  the  spin  pin  is  unity,  then  the  coordinates  of  the  end  point  are  Dx  ,Dy  X>z-  Then 
the  projection  of  the  vector  endpoints  on  the  Y-Z  plane  is  at  X=0,  Y=Dy,  and  Z=Dz. 
Where,  Dy  =  cos  p  =  mp  and  =  cos  y  =  np .  The  angle  ^  shown  in  the  sketch  is  the 
angle  of  the  D  vector  projected  on  the  Y-Z  plane.  Then, 

tan  (|)  =  Dy  /  Dz  =  mp/  np 
or 

(j)  =  tan*’  (mp  /  Op )  5.33 

Equation  Eq.  5.33  is  the  relationship  used  to  determine  the  model’s  roll  orientation  at 
each  of  the  shadowgraph  stations. 


5.2.4  Center-of-Gravity 


In  Section  5.2.2  the  model’s  nose  was  chosen  to  illustrate  the  method  used  to 
determine  the  range  coordinates  of  specific  points  on  the  body.  However,  what  the  data 
reduction  routines  require  are  the  range  coordinates  of  the  model’s  center-of -gravity. 
This  is  the  point  about  which  the  equations  of  motion  have  been  derived  and  therefore  is 
the  point  for  which  the  X,Y,  Z  position  measurements  are  required.  One  of  the  pretest 
measurements  as  discussed  in  Section  3.5.1  was  the  distance,  L(.g,  between  the  nose  and 
the  center-of  gravity  of  the  model.  This  measurement  in  conjunction  with  the  range 
coordinates  of  the  nose  as  obtained  using  Eqs.  5.3 1  and  the  direction  cosines  from  Eqs. 
5. 14  is  used  to  compute  the  range  coordinates  of  the  model’s  center-of-gravity.  These 
relationships  are  shown  below:  s 


^cg  ~  “  Teg  (  Pe  ) 

Yjg  —  Yj,|  -  Egg  ( mg) 
Zgg  ~  Zfvj  -  Egg  ( Og ) 


i  5.34 


The  range  coordinates  of  the  model’s  center-of  gravity,  Xgg ,  Ygg  ,  Zgg ,  are  also 
determined  at  each  shadowgraph  station  and  provided  to  the  data  reduction  routines  for 
analysis  of  the  aerodynamic  coefficients  and  derivatives. 


5.3  TIME  OF  FLIGHT  MEASUREMENTS 


The  present  chronograph  system  provides  the  required  times  of  flight  as  the  projectile 
traverses  the  instrumented  range.  Event  times  corresponding  to  each  spark  source 
discharge,  hall  and  pit,  are  obtained  to  a  resolution  of  ±  0.01  jxsec  by  electronic 
chronographs  at  each  of  the  50  shadowgraph  stations.  The  chronographs  at  each  station 
operate  under  the  control  of  the  sequencer  and  in  conjunction  with  the  IR  detection 
system  and  the  spark  gap  assembly.  The  timing  system  in  all  chronographs  is 
synchronized  to  a  master  10  Mhz  clock  located  in  the  control  console.  The  sequencer 
provides  the  pulses  listed  below  to  the  chronographs  at  the  following  times  in  the  firing 
sequence: 

(1)  -3.5  seconds,  solenoid  operated  camera  shutters  open 
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(2)  -2.0  seconds,  clear  (resets  all  chronograph  counters  to  zero) 

(3)  -1.5  seconds,  start  (counters  synchronized  to  master  10  Mhz  clock) 

(4)  -1 .0  seconds,  fire  (fire  pulse  sent  to  launcher) 

(5)  +3.5  seconds,  solenoid  operated  camera  shutters  closed 

(6)  +10.5  seconds,  strobe  (accumulated  times  transferred  and  strobed  on  film) 

When  the  projectile  is  detected  in  the  vicinity  of  a  shadowgraph  station  a  15-volt 
pulse  initiates  the  discharge  of  the  capacitor.  A  light  sensing  diode  located  inside  each  of 
the  spark  gap  cavities  senses  the  discharge  and  subsequently  stops  the  chronographs  for 
both  the  hall  and  pit  counters.  This  results  in  the  total  elapsed  time  from  when  the 
common  start  pulse  was  sent  to  the  arrival  of  the  projectile  at  the  shadowgraph  station.  A 
sample  printout  of  the  accumulated  times  for  both  the  hall  and  pit  cameras  were  shown 
previously,  see  Table  4.1.  The  hall  and  pit  times  are  averaged  to  obtain  the  time  of  flight 
which  is  provided  to  the  data  reduction  routines  along  with  the  position-attitude  data  as 
discussed  in  the  preceding  section. 

Some  of  the  other  data  that  were  listed  in  Table  4.1  should  be  mentioned.  For 
example  the  raw  velocity  profile  is  computed  using  the  averaged  hall  and  pit  times  and 
the  nominal  distances  between  stations  as  listed  in  Appendix  1 .  This  raw  velocity  profile 
is  also  ‘fitted’  using  a  polynomial  least  squares  curve  fit  and  the  resulting  smoothed 
velocity  profile  also  shown.  These  velocity  profiles  are  provided  in  order  to  assist  the 
facility’s  operational  crew  and  test  engineer  in  an  immediate  evaluation  of  that  particular 
flight  and  are  not  used  in  the  final  data  analysis  routines.  Also  included  in  Table  4.1  is  a 
flag  (FI)  indicating  whether  or  not  the  hall  and  pit  times  are  within  1  psec  of  each  other. 
Variations  in  these  times  greater  than  1  psec  should  be  examined  by  the  operations  crew 
and  the  appropriate  adjustments  made.  The  column  of  drag  coefficients  is  determined 
from  the  smoothed  velocity  profile  and  is  presented  as  a  quick  look  evaluation 
mechanism  and  not  as  a  final  determination  of  the  drag  coefficient. 

5.4  TRAJECTORY  OUTPUT 

The  final  trajectory  data  determined  using  the  techniques  previously  discussed  are 
provided  to  the  aerodynamic  data  analysis  routines  as  computed  using  CADRA  .  A 
sample  raw  data  file  is  shown  as  Table  5.2.  The  time  of  flight  data  were  measured  directly 
as  discussed  in  the  previous  section  where  the  time  at  the  first  station  was  subtracted  from 
each  of  the  following  times  yielding  a  time  profile  with  zero  time  at  the  first  station.  The 
x,y,z  position  columns  were  obtained  using  Eqs.  5.34  and  the  roll  orientation,  <|),  was 
computed  from  Eq.  5.33.  The  fixed  plane  angles  are  derived  from  the  direction  cosines, 

( m^,  n^,  Pe ),  as  shown  in  Eqs.  5.14  and  the  fixed  plane  identities  as  listed  in  Section 
2.4.1,  see  Eqs.  2.8.  Or, 

2  2  1/2 

sin0  =  -ne  and  sin\|/=  me/(me  +Pe  )  5.35 

This  now  completes  the  methodology  of  determining  the  experimentally  measured 
trajectories  which  are  provided  to  the  analysis  routines  for  extraction  of  the  aerodynamic 
parameters. 


71 


Table  5.2 

Raw  Trajectory  Data  From  CADRA 


950704  ISL/WAF 

Range  data  file  saved  in  S9 50704. xyz 

Average  Reynolds  Number  =  2050000;  Average  Mach  Number  =  1.48 


station 

time 

V(m/sec)  x{m) 

y{m) 

z(m) 

-sin(\|;) 

sin(0) 

1 

0.0000000 

618.26 

2.12497 

1.41345 

0.73105 

0.10361  -0.15549 

226.66 

2 

0.0025025 

615.98 

3.66820 

1.41971 

0.73672 

-0.07315  -0.18116 

181.90 

3 

0.0050034 

613.71 

5.20546 

1.42374 

0.73858 

-0.19982  -0.07842 

147.92 

4 

0.0099451 

609.29 

8.22347 

1.42528 

0.74674 

-0.05133 

0.18431 

30.82 

5 

-0.2514225 

0.00 

14.28941 

1.43166 

0.77003 

0.00278  -0.18781 

151.52 

6 

0.0251753 

596.05 

17.39730 

1.43791 

0.77688 

-0.17858 

0.02295 

342.49 

7 

0.0302340 

591.78 

20.39437 

1.43961 

0.78392 

0.05431 

0.17573 

168.58 

8 

0.0328044 

589.63 

21.91187 

1.43994 

0.79007 

0.16627 

0.08201 

92.74 

9 

0.0354859 

587.41 

23.48780 

1.44302 

0.79432 

0.16190  -0.07419 

341.13 

10 

0.0459116 

578.93 

29.55916 

1.45214 

0.81134 

-0.11680 

0.13258 

345.73 

11 

0.0566330 

570.45 

35.70624 

1.46220 

0.82884 

0.06102  -0.15324 

286.38 

12 

0.0619041 

566.37 

38.69989 

1.46673 

0.83496 

-0.16408  -0.01770 

113.81 

13 

0.0673124 

562.25 

41.74600 

1.47039 

0.83913 

0.01794 

0.17214 

285.94 

14 

0.0727530 

558.17 

44.79242 

1.47272 

0.84591 

0.15227  -0.04873 

88.69 

15 

0.0754616 

556.15 

46.30227 

1.47567 

0.85318 

0.06332  -0.13776 

357.28 

16 

0.0865332 

548.07 

52.40311 

1.48547 

0.86514 

0.02475 

0.15914 

333.02 

17 

0.0921448 

544.07 

55.46419 

1.48800 

0.87129 

0.15005  -0.04899 

145.67 

18 

0.0949342 

542.10 

56.98063 

1.49013 

0.87708 

0.03830  -0.10338 

73.39 

19 

0.1091995 

532.25 

64.62859 

1.49987 

0.88627 

0.15391 

0.05611 

296.69 

20 

0.1150409 

528.32 

67.72624 

1.50745 

0.89362 

0.01094  -0.16405 

145.54 

21 

0.1207103 

524.56 

70.71324 

1.51155 

0.89643 

-0.12654 

0.03632 

339.91 

22 

0.1264541 

520.80 

73.71564 

1.50739 

0.90060 

0.08501 

0.12370 

166.66 

23 

0.1323776 

516.98 

76.79142 

1.52084 

0.90080 

0.08317  -0.10487 

7.21 

24 

0.1382248 

513.27 

79.80318 

1.52611 

0.90757 

-0.11228  -0.04451 

210.44 

25 

0.1442620 

509.49 

82.89021 

1.52940 

0.90714 

0.01446 

0.14294 

47.99 

26 

0.1472299 

507.65 

84.40366 

1.52956 

0.90960 

0.12048 

0.09482 

333.54 

27 

0.1562597 

502.14 

88.96892 

1.53878 

0.91372 

-0.08002  -0.09976 

144.69 

28 

0.1654281 

496.67 

93.55601 

1.54535 

0.91592 

0.06102 

0.12929 

295.23 

29 

0.1746733 

491.27 

98.13047 

1.55271 

0.91964 

-0.00966  -0.12632 

113.88 

30 

0.1840232 

485.93 

102.71127 

1.55765 

0.91537 

-0.00296 

0.12628 

285.85 

31 

0.1935652 

480.59 

107.32884 

1.56425 

0.92208 

0.03112  -0.12796 

113.04 

32 

-0.0922876 

0.00 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

180.00 

33 

0.2225252 

465.10 

121.05094 

1.58288 

0.91885 

-0.07198 

0.12088 

339.86 

34 

0.2322087 

460.14 

125.54861 

1.59135 

0.91409 

0.11390  -0.09023 

159.34 

35 

0.2422016 

455.13 

130.12985 

1.59589 

0.91235 

-0.11643 

0.07853 

358.38 

36 

0.2524416 

450.11 

134.77364 

1.60659 

0.90783 

0.13468  -0.08008 

160.02 

37 

0.2625836 

445.24 

139.32494 

1.61151 

0.90246 

-0.13006 

0.07177 

348.46 

38 

0.2726660 

440.51 

143.79599 

1.61415 

0.89701 

0.15049  -0.05835 

185.87 

39 

0.2832808 

435.64 

148.44504 

1.62340 

0.88753 

-0.14166 

0.06864 

342.70 

40 

0.2938210 

430.90 

153.01378 

1.62996 

0.88074 

0.16352  -0.07236 

166.57 

41 

0.3045064 

426.20 

157.59068 

1.63824 

0.87054 

-0.13382 

0.08178 

354.85 

42 

0.3153282 

421.55 

162.17097 

1.64027 

0.86018 

0.15472  -0.09490 

152.74 

43 

0.3263592 

416.91 

166.78352 

1.64999 

0.84774 

-0.12224 

0.12451 

322.49 

44 

0.3374380 

412.35 

171.35765 

1.65641 

0.83529 

0.12631  -0.12955 

144.50 

45 

0.3561066 

404.89 

178.94235 

1.66385 

0.81138 

0.18168  -0.01659 

206.80 

46 

0.3675855 

400.43 

183.53200 

1.67269 

0.79490 

-0.16948 

0.05072 

347.20 

47 

-0.4100693 

o:oo 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

180.00 

48 

0.3906884 

391.76 

192.59471 

1.68348 

0.75577 

-0.13874 

0.12300 

347.60 

49 

-0.0101338 

0.00 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

180.00 

50 

0.4147032 

383.13 

201.77632 

1.69648 

0.71092 

-0.06531 

0.17812 

318.64 
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VI.  TRAJECTORY  ANALYSIS 


Trajectory  analysis  as  referred  to  in  ffee-flight  testing  is  the  process  by  which  the 
aerodynamic  coefficients  and  derivatives  are  obtained  from  the  experimentally  measured 
trajectory.  This  process  consists  of  fitting  the  theoretical  equations  of  motion  to  the 
experimental  position-attitude-time  measurements  assuming  that  the  aerodynamic 
coefficients  and  derivatives  are  unknowns  in  the  process.  There  are  various  techniques 
that  have  been  used  to  accomplish  this  process.  The  two  techniques  to  be  described 
herein  are  generally  referred  to  as  ‘Linear  Theory’  and  ‘Numerical  Integration’. 

6.1  LINEAR  THEORY 

When  Charters’  *  first  began  to  extract  aerodynamic  coefficients  and  derivatives  from 
free-flight  spark  range  data  at  BRL  in  the  1940’s  computer  technology  was  in  its  infancy. 
This  lack  of  computer  power  necessitated  that  the  equations  of  motion  be  derived  in  a 
manner  such  that  closed  form  solutions  could  be  written  from  the  differential  equations. 
These  closed  form  solutions  were  then  used  in  the  fitting  process  to  determine  the 
associated  aerodynamics.  In  fact  in  those  early  days  the  fitting  process  usually  consisted 
of  analog  computer  techniques  or  frequently  graphical  methods  of  extracting  the 
aerodynamics.  In  writing  these  closed  form  solutions  it  was  necessary  to  make  various 
assumptions  which  will  be  highlighted  in  the  following  sections.  The  two  primary 
assumptions  inherent  in  these  derivations  are  symmetry  and  linear  aerodynamics.  Linear 
aerodynamics  means  that  derivatives  of  the  aerodynamic  moment  and  force  coefficients 
with  respect  to  the  angle-of-attack  are  constants,  i.e.,  vs  a  can  be  represented  as  a 
straight  line.  This  assumption  of  linearity  also  means  that  small  angles  are  usually 
assumed. 

The  various  linear  theory  techniques  discussed  herein  were  developed  by  numerous 
individuals.  However,  the  two  people  primarily  associated  with  these  techniques  are  Dr. 
Charles  Murphy®  *’  ®  ^  and  Dr.  John  D.  Nicolaides®  Many  of  the  derivations  and  much 
of  the  discussion  in  the  following  sections  were  extracted  from  these  and  other  references, 
i.e.,  also  see  Refs.  6.5  and  6.6.  These  linear  theory  techniques  provide  insights  into  the 
free-flight  dynamics  that  cannot  be  achieved  using  the  complex  numerical  integration 
methods  that  will  be  discussed  later  in  Section  6.3.  Linear  theory  is  also  typically  used  to 
provide  the  initial  estimates  of  the  aerodynamic  coefficients  and  derivatives  required  by 
the  numerical  integration  methods.  For  many  test  configurations  and  conditions  where  the 
inherent  assumptions,  required  in  deriving  the  linear  theory  equations  of  motion,  are  not 
severely  violated  these  techniques  provide  very  accurate  estimates  of  the  resulting 
aerodynamics. 

6.1.1  Drag 

The  total  aerodynamic  drag  coefficient,  C^,  is  obtained  from  an  analysis  of  the 
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distanced  traveled,  x,  verses  time  measurements.  The  total  drag  coefficient  is  defined  as 
Cd  =  D  /  (1/2)  p  A  and  applying  Newton’s  second  law,  D  =  -  m^  a ,  the  following 
relationship  can  be  written; 

CD=-(2ma/pAV^)a 

Since  the  deceleration,  a,  can  be  expressed  as  V  ( dV  /  dx )  the  above  relationship 
becomes, 

CD  =  -(2ma/p  AV)dV/dx  6.1 

If  the  fimctional  dependence  between  t  and  x  is  assumed  to  be  a  cubic  polynomial,  i.e. 

t  =  Cq  +  CiX  +  +  C3X^ 


then 

(dt  /dx)  =  ( 1/V )  -  Cl  +  2c2X  +  3c3X^ 

6.2 

and 

(dV/dx)  =  -  (  2c2  +  6C3X )  /  (ci  +  2c2X  +  3c3X^)^ 

6.3 

Substituting  Eqs.  6.2  and  6.3  into  Eq.  6.1  produces  the  following  equation: 

Cd  =  (2  nia  /  p  A)  [(2c2  +  6C3X)  /  (cj  +  2c2X  +  3C3X  )] 

6.4 

Equation  6.4  is  the  common  linear  theory  relationship  used  in  computing  the  total  drag 
coefficient.  The  process  is  to  least  squares  fit  the  measured  time  and  distance  traveled 
data  using  the  cubic  polynomial  as  shown.  Once  the  coefficients  (Cj’s)  of  the  polynomial 
have  been  determined  they  are  substituted  into  Eq.  6.4  and  the  Cp  evaluated 
at  the  midrange  x  location. 


Occasionally  the  analyst  is  provided  velocity  data  in  lieu  of  time  vs  distance.  This  can 
come  from  velocity  screens  where  the  At  and  Ax  measurements  are  used  to  compute 
average  velocities  between  screens.  The  total  drag  coefficient  can  be  determined  from 
these  measurements  using  the  familiar  relation: 


V  =  Voe‘’^  where,  K  =  CDpA/2ma 


6.5 


ln(V/Vo)  =  -Kx  6.6 

Graphically  this  can  be  described  as  plotting  V  /  Vq  on  a  log  scale  vs  x  (semi-log  paper) 
and  fitting  the  measurements  with  a  straight  line  as  shown  on  the  following  sketch. 


The  slope  of  the  straight  line  is  K,  see  Eq.  6.6,  thus  Cp  =  K  ( 2  mg  /  p  A ).  The  term 
“ballistic  coefficient”,  p,  is  frequently  used  and  this  is  generally  defined  as  the  projectile 
weight  divided  by  Cp  A,  or 

p  =  g  ma  /  Cd  A 

Equation  6.5  is  frequently  used  in  estimating  the  velocity  of  an  object  as  a  function 
of  distance  when  the  drag  coefficient  is  assumed  to  be  known.  This  relation  provides 
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accurate  velocity  estimates  even  when  relatively  large  distances  are  involved  or  large 
decelerations  are  present.  For  those  cases  where  the  velocity  decay  is  small,  less  than  10 
percent,  Eq.  6.1  is  frequently  used  directly.  Historically  Eq.  6.1  was  used  because  of  the 
simplicity  of  plotting  velocity  vs  distance,  drawing  a  straight  line  through  the  data, 
computing  the  slope  of  the  line  and  calculating  Cq  assuming  that  dV/dx  =  AV/Ax.  This 
could  be  accomplished  graphically  with  ease  and  precision  even  prior  to  the  advent  of 
computer  technology,  see  sketch  below. 


All  of  the  techniques  mentioned  above  result  in  the  estimate  of  the  total  drag 
coefficient  of  the  free-flight  object.  If  that  object  is  flying  at  a  non-zero  angle-of-attack, 
normally  an  oscillatory  angular  motion  is  present  and  possibly  trimmed,  then  the 
measured  Cp  is  affected  (usually  increased)  by  this  non-zero  angle-of-attack  flight 
condition.  This  change  in  caused  by  the  object’s  angle-of-attack  can  be  estimated  and 
the  total  drag  coefficient  adjusted  yielding  the  zero  angle-of-attack  drag  coefficient,  Cdq- 
This  process  will  be  described  in  Section  6.2.1 . 

6.1.2  Roll 

The  analysis  of  the  missile’s  rolling  motion  is  accomplished  using  the  measured  roll 
orientation,  (j),  and  x  data.  Since  the  roll  moment  is  defined  as  C]  =  1  /  (1/2)  p  A  d  and 
the  rotational  equivalent  of  Newton’s  second  law  is  1  =  Ix  it  follows  that: 

Ix^  =(l/2)pV^  AdC,  6.7 

If  we  assume  that  C|  is  comprised  of  a  driving  moment,  due  to  fin  cant,  and  a  roll 
damping  moment,  caused  by  skin  friction  and  the  paddle  wheel  effect  of  the  fins,  then  C| 
can  be  expanded  as, 

C,-5C,6  +  (^  d/2V)C,p  6.8 

Substituting  Eq.  6.8  into  6.7  yields: 

Ix^  =  (1/2)  p  A  d  [5  C|8  +  (^  d  /  2V  )  C,p]  6.9 

Since  we  have  both  ^  and  time  measurements  a  closed  form  solution  to  Eq.  6.9  could 
be  written  and  in  conjunction  with  a  curve  fitting  process  the  unknown  aerodynamic 
coefficients  C15  and  Cip  determined.  This  closed  form  solution  would  require  assuming 
that  the  coefficients  of  the  differential  equation  are  constant.  As  a  result  of  this 
assumption  the  velocity  of  the  missile  would  also  be  assumed  to  be  constant,  or  Cp  =  0. 
This  very  restrictive  assumption  can  be  relaxed  by  changing  the  independent  variable  in 
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Eq.  6.9  from  time  to  distance.  This  will  permit  the  velocity  to  vary  and  hence  Cp  to  be 
non-zero.  Eventually  however,  we  will  be  required  to  assume  that  Cq  is  a  constant. 
Nevertheless,  it  is  believed  that  this  is  a  prudent  change  of  variables  and  is  accomplished 
as  shown  below.  Since 

<!>  =  d<|)/dt  =  (d(|)/dx)(dx/dt)  =  V 
then  (f)  can  be  expressed  as, 

^  =  (!>"  +  F'  ^ 

Substituting  the  above  relationships  for  (p  and  p  along  with  V  as  obtained  from  Eq.  6.1 
into  Eq.  6.9  and  collecting  like  terms,  one  can  arrive  at  the  following  differential 


equation: 

(Zi"  +  Kp^'-Kg  =  0 

6.10 

where 

Kp  =  -(p  A/  2m  Cd  +  (d^  I,)  C,p] 

and 

K5  =  (pAd/2y6C,8 

Assuming  that  Kp  and  Kg  in  Eq.  6.10  are  constants,  or  Co,  Cjp,  and  Cjg  are  constants, 
the  two  roots  of  the  characteristic  equation  are  zero  and  -  Kp.  The  particular  solution  is 
(Kg  /  Kp)  X  and  after  some  manipulation  the  solution  of  Eq.  6.10  can  be  written  as  shown 
below: 


(l)  =  (t)o  +  (K5/Kp)x  +  C(e-‘"’''-l)  6.11 

The  unknowns  in  Eq.  6.1 1  which  require  determination  are  (|)o.  Kg,  Kp  and  the  constant  C. 
Because  Kp  appears  in  a  nonlinear  fashion,  the  process  of  fitting  the  above  equation  to  the 
measured  and  x  data  requires  using  a  differential  corrections  technique.  This  is  the  same 
process  used  in  fitting  the  oscillatory  pitch  and  yaw  motions,  see  Appendix  5.  Historically 
this  complexity  was  eliminated  by  approximating  the  e'*^’’^  term  with  the  first  four  terms 
of  the  infinite  series  expansion,  i.e., 

e-Kpx  =  1  .  Kp  X  +  (1/2)  Kp^  x^  -  (1/6)  Kp^  x^ 

Substituting  this  expansion  into  Eq.  6.1 1  and  collecting  like  terms,  one  arrives  at  the 
equation  shown  below: 

(|)  =  (|)0  +  (|)i  X  +  <t»2  x^  +  (j)3  x^  6.12 

where,  <j),  =  (Kg  /  Kp)  -  C  Kp 

(|)2  =  (l/2)CKp2 

(1.3  =  -(l/6)CKp' 
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The  unknown  (t)i’s  in  Eq.  6.12  can  now  be  determined  using  a  simple  least  squares 
fitting  process.  Then  C,p  and  5C,5  are  determined  using  the  <|)i ,  (jia ,  and  relations  (Eqs. 

6. 1 2)  in  conjunction  with  the  Kp  and  K5  definitions  (Eqs.  6. 1 0)  arriving  at; 

C,p  =  ( Ix/ d^  ma)[  (6  (t)3  ma/ (1)2  P  A)  -  Cd ] 
and  6C,5  =  (2Ix/pAd)[2(|)2-(3(|)3(|),/(l)2)  6.13 

Both  Eqs.  6.10  and  6.13  demonstrate  that  prior  to  computing  the  roll  damping  derivative, 
C|p,  the  drag  coefficient  must  be  known. 


6.1.3  Pitch  and  Yaw 

The  differential  equation  governing  the  angular  oscillatory  motion  of  a  slightly 
asymmetric  missile  as  a  function  of  distance  traveled  has  been  rigorously  derived  in  Refs. 
6.1  to  6.4  and  is  not  repeated  herein.  However,  this  differential  equation  is  shown  below: 

(p  +  ia)”  +  (H  -  iP)  (P  +  ia)'  -  (M  +  iPT)  (p  +  ia)  =  iA^  e*^”  6.14 


where  H  =  (pA/2ma)  [  -  2Cd  -  k,'^(  C^q  +  Cn,&)]  6. 1 5 

M  =  (pA/2mad)kt-'C„,„  6.16 

T  =  (pA/2ma)  [  CNa  -  Cd  +  ka'"  C„p„  ]  6.17 

P  =  (Ix/Iy)p  6.18 


At=  (pA/2mad)  {k,-'(  Cn,o  +  iCno)  +  Pd[(  Ix  /  ly )  -1](Cyo+  iCzo)}  6.19 


Also  ka  and  kt  are  the  nondimensionized  axial  and  transverse  radiuses  of  gyration 
respectively,  or 

ka  =  ( Ix  /  ma  d^  )'^  and  k,  =  ( ly  /  ma  d^ 

The  closed  form  solution  to  Eq.  6.15  assuming  that  H,  M,  T,  P,  A^  and  the  roll  rate, 
p,  are  all  constant  can  be  immediately  written  as  shown  below: 

P  +  ia  =  Kp  exp[(pp  +  i(t)p')x]  +  K^  exp[(pN  +  i<|)N')x]  +  Kj  exp(ipx)  6.20 

Equation  6.20  is  frequently  referred  to  as  the  ‘tricyclic  equation  of  motion’  because  it  can 
be  graphically  represented  by  three  rotating  vectors  in  the  complex  P  +  ia  plane,  as 
illustrated  in  the  sketch  below. 
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This  sketch  demonstrates  that  the  trim  vector,  Kj,  has  a  constant  amplitude  and  rotates  at 
the  spin  rate,  p.  The  nutational  vector,  K^,  rotates  with  a  rate  of  (j)^'  and  changes 
amplitude  at  a  rate  of  Likewise,  the  precessional  vector  rotates  at  a  rate  of  (j)p'  and 
changes  amplitude  at  a  rate  of  pp.  Note  that  if  the  p’s  are  negative  both  vectors  decrease 
in  amplitude  and  the  motion  damps.  By  definition  (p^H-  icl)^')  and  (pp+i(|)p')  are  the  roots 
of  Eq.  6,14.  Also,  the  assumption  that  H,  M,  and  T  are  constant  necessitates  that  the 
aerodynamic  coefficients  are  linear  with  angle-of-attack  and  hence  the  name  ‘Linear 
Theory’. 

This  ‘Linear  Theory’  technique  of  analyzing  the  angular  motion  of  a  missile  in  ffee- 
flight  requires  that  the  measured  P  +  ia  data  be  fitted  as  a  fimction  of  distance  traveled,  x, 
and  determine  K^,  Kp,  Kj,  p^,  pp,  and  <|)p'  in  Eq.  6.20  assuming  that  they  are  the 
unknowns  in  the  process.  This  fitting  process  is  derived  and  discussed  in  Appendix  3. 
Therefore,  in  this  section  we  will  proceed  as  if  the  Kj’s,  pj’s,  and  <t>j’s  have  already  been 
determined . 

We  can  relate  the  determined  pj’s  and  ())j’s  to  the  aerod5mamic  coefficients  as  defined 
in  the  constants,  Eqs.  6.15  to  6.18,  of  the  differential  equation,  Eq.  6.14,  by  analyzing  the 
roots  of  the  characteristic  equation  using  the  operator  method.  Or, 

+  (H  -  iP)  D  -  (M  +  iPT)  =  0 


then  the  roots  are , 

Di,2  =  (1/2)  [  -  (H  -iP)  ±  V(H  -  iP)^  +  4(M  +  iPT)]  6.21 

By  adding  the  two  roots  we  can  obtain. 


D,+D2  =  -(H-iP) 


6.22 


80 


and  by  multiplying  the  two  roots, 


6.23 


D1D2  =  -  (M  +  iPT) 

By  using  the  definition  of  the  roots  we  can  also  write, 

Di  +  D2  =  ( |J.N  +  l^p)  +  i(<l>N'  +  <t>p') 

and  D1D2  =  Pn  |4.p  -  <|)n^  (t>p'  +  i(<l>N'  Pp  *t*p  M-n) 

If  we  now  equate  the  real  and  imaginary  parts  of  Eq.  6.22  with  the  real  and  imaginary 
parts  of  Eq.  6.24,  after  some  manipulation  and  using  Eqs.  6.16  through  6.19  we  can  arrive 


at  the  following  equations: 

-  (2ma  /  pA)(  Pn  +  Pp)  =  Cwa  -  2Cd  -  ( C^q  +  6.26 

and  P  =  (Iy/Ix)(V  +  't)p')  6.27 

Likewise  if  we  equate  the  real  and  imaginary  parts  of  Eq.  6.23  and  6.25,  and  using  Eqs. 
6.15  through  6.18,  after  manipulation  we  obtain: 

*l*p^  *  M-n  Pp)  ( 2Iy  /  pAd )  =  6.28 

and  ((jiiq'  Pp  +  (jip'  Pn)  /  (<!>n'  +  ‘t>pO(  pA  /  2ma)  =  Cp  -  Cno  -  kg  ^  Cnp„  6.29 


Equations  6.26  through  6.29  illustrate  that  the  aerodynamic  coefficients  and 
derivatives  are  a  function  of  the  ‘fitted’  pj’s  and  <|)i’s  of  Eq.  6.20.  The  drag  coefficient 
would  be  provided  using  the  techniques  of  Section  6.1.1  and  the  normal  force  derivative, 
Cnoi,  using  the  method  discussed  in  the  next  section.  Eqs.  6.26, 6.28,and  6.29  would  then 
be  used  to  determine  Cn,q+  ,  C^a ,  and  C„p„ .  It  is  of  interest  to  note  that  the  spin  rate, 
p,  of  Eq.  6.27  is  a  function  of  the  nutational  and  precessional  frequencies.  For  this  reason 
the  spin  rate  is  not  treated  as  an  unknown  in  Eq.  6.20.  It  should  also  be  noted  that 
normally  the  product  of  the  damping  rates,  p^  pp ,  in  Eq.  6.28  is  assumed  negligible  and 
the  pitching  moment,  C^a  >  is  written  only  as  a  function  of  the  frequencies.  This  was  a 
particularly  useful  simplification  in  the  days  when  graphical  techmques  were  used. 
Although  this  product  is  very  small  and  ignoring  it  results  in  a  negligible  error  to  C^a 
there  is  really  no  reason  why  this  assumption  still  persists. 

Most  linear  theory  routines  in  use  today  have  been  modified  from  that  shown  as  Eq. 
6.20  by  allowing  the  nutational  and  precessional  frequencies  to  be  linear  functions  of 
distance.  These  are  assumed  to  be, 

<|)p'  =  (|)'po  +  <|)p"  X  and  (|)n'  =  (j)'No  +  (|>n”  x  6.30 
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If  we  substitute  Eqs.  6.30  into  Eq.  6.27  the  spin  rate  also  becomes  a  linear  function  of  x. 
Or, 


P  =  (  ly  /  Ix)  [f  NO  +  f  PO  +  X]  6.3 1 

The  aerodynamic  coefficients,  and  ,  as  defined  by  Eqs.  6.28  and  6.29  then 
become  functions  of  x  and  are  normally  evaluated  at  the  midrange  ( Xnu- )  location  using 
Eqs.  6.30.  This  variable  frequency  expansion  to  the  classic  ‘Linear  Theory’  technique  is  a 
powerful  tool  and  is  often  utilized  by  the  analyst  as  a  first  order  approximation  in 
accounting  for  nonlinear  and  variable  roll  effects  in  the  experimentally  measured  motion 
profiles.  Notice  that  the  variable  frequencies  do  not  affect  the  computation  of  Cn,q+C„^; 
from  Eq.  6.26.  However,  if  the  wavelength  of  the  motion  is  not  well  matched  during  the 
fitting  process  this  will  usually  result  in  significant  errors  in  the  fitted  damping  rates,  Pj’s, 
and  hence  cause  significant  errors  in  the  computed  C^q+C^i  values.  Therefore  the  use  of 
this  variable  frequency  expansion  helps  to  ensure  that  the  errors  in  the  computed 
aerodynamics  are  minimized. 

It  should  be  recognized  that  the  inclusion  of  these  variable  frequencies  is  strictly  an 
intuitive  approach  and  does  not  represent  an  actual  mathematical  solution  to  the 
differential  equation.  In  fact,  the  variable  spin  rate  and  the  variable  frequencies  violate  the 
assumption  that  H,  M,  P,  T,  and  p  are  all  constant  when  writing  the  closed  form  solution 
to  the  differential  equation. 

Some  linear  theory  routines  have  also  been  modified  to  permit  the  trim  vector  to 
change  magnitude  by  including  a  damping  rate,  Px ,  on  the  trim  term  of  Eq.  6.20.  This 
extension  of  the  classic  ‘Linear  Theory’  can  also  be  useful,  especially  in  analyzing  free- 
flight  trajectories  for  missiles  which  may  be  approaching  a  spin-yaw  resonate  condition, 
p->(t)N'.  When  this  condition  occurs  the  aerodynamic  trim  increases  in  amplitude  as 
shown  in  the  sketch  below. 
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This  increase  in  amplitude  is  described  by  the  amplification  factor,  AF,  and  can  be 
shown  to  be  defined  by  the  relationship  below: 

AF  « -( m,d  /  Iy(t)N')  /  [  Cn„  -  (m/  /  ly)  (€„,,+  C^)]  6.32 

For  most  missile  configurations  this  amplification  factor  turns  out  to  be  a  value  between 
8  and  14.  Because  of  this  amplification  effect  the  inclusion  of  a  damping  rate  on  the  trim 
term  of  Eq.  6.20  can  sometimes  assist  in  fitting  the  tricyclic  equation  of  motion  to  the 
experimentally  measured  trajectories.  If  the  model  actually  passes  through  resonance 
where  the  trim  initially  grows  and  then  decreases  a  simple  damping  rate  added  to  the  trim 
term  will  not  appreciably  aid  the  fitting  process.  Other  problems  associated  with  this 
resonance  condition  can  also  arise  when  fitting  the  angular  motion.  For  example  by 
definition  this  condition  means  that  both  the  nutational  and  trim  vectors  are  rotating  in  the 
same  direction  at  about  the  same  rate.  Mathematically  this  motion  can  be  described  by  a 
single  vector  rotating  at  near  resonance  and  the  fitting  routine  can  and  frequently  does 
diverge  because  of  the  absence  of  a  unique  solution.  Nevertheless,  various  analysts  have 
taken  considerable  liberties  with  the  classic  ‘Linear  Theory’  over  the  years  and  one 
should  be  alert  for  these  modifications. 


6.1.4  Swerve 

In  contrast  to  the  rigorous  derivation  for  the  angular  motion  in  Refs.  6.1  to  6.4  here 
we  will  do  an  intuitive  derivation  of  the  swerving  motion  starting  with  Newton’s  second 
law,  a  =  F/ma  ( see  Ref.  6.1  for  a  complete  derivation).  We  will  use  the  complex 
nomenclature  similar  to  Eq.  6. 1 5,  or  a  =  (y  +  iz),  then 

(y  +  iz)  =  ( 1  /2)  p  A  Cp  /  ma  6.33 

The  force  coefficient ,  Cp,  can  be  expanded  as  a  function  of  the  lift  force  coefficient.  Cl, 
and  a  trim  force  which  rotates  with  the  spin  rate,  or 

Cp  =  Cl  +  (Cyo  iCzo)  ®  ^  6.34 

The  Cl  coefficient  is  further  expanded  as  CLa(  P  +  ia)  and  substituting  this  and  Eq.  6.34 
into  Eq.  6.33  we  obtain: 

(y  +  iz)  =  (pA/  2ma)  V'  [Cl^CP  +  ia)  +  (Cyo  +  iCzo)  6.35 

If  we  now  use  the  change  of  variables  technique  as  was  done  for  the  roll  equation, 
Eq.  6.9,  Eq.  6.35  can  be  written  as, 

(y  +iz)"  =  ( pA  /  2ma)  [CLa(P  +  ia)  +  (Cyo  +  iCzo)  6.36 
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Now  if  we  substitute  Eq.  6.20  for  (  p  +  ia),  collect  like  terms,  integrate  twice  by  parts, 
and  evaluate  from  0  to  x,  after  some  manipulation  we  can  arrive  at  the  equation  shown 
below: 

(y  +  iz)  =  (  pA  /  2ma){CLa  [  Kn  exp<|)N  x)  +  (t)p‘^(  Kp  exp(|)p  x)  -  (Km  +  Kp  (j)p'‘)  x 

-  (Kn +  Kp <j)p'^)]  +  (  Kt  CLa  +  Cyo  +iCzo)(p'^)[  1-  cos  px  +  i(px  -  sin  px)]} 

+  (y  +  iz)'o  X  +  (y  +  iz)o  6.37 

Where:  (j)p  =  pp  +  i(|)p'  ,  (j)^  =  Pn  ^nd, 

(y  +  iz)'o  =  initial  velocities  in  the  y  and  z  directions 

(y  +  iz)o  =  initial  displacements  in  the  y  and  z  directions 

In  arriving  at  the  above  equation  it  was  also  necessary  to  assume  that  ())p  ,<j)N  and  p  were 
not  equal  to  zero. 

Equation  6.37  is  fit  to  the  measured  position  data,  y  +  iz,  assuming  that  Ci^  ,  Cyo  , 
Czo  and  the  initial  velocities  and  displacements  are  unknown.  Note  that  the  Kj’s,  (|)i’s,  and 
p  are  assumed  to  be  known  and  provided  by  the  analysis  discussed  in  the  preceding 
section.  It  should  be  noted  that  for  spin  stabilized  projectiles  it  is  normally  assumed  that 
Kp,  Cyo,  and  Czq  are  zero  and  Eq.  6.37  is  somewhat  simplified. 


6.2  QUASI-NONLINEAR  ANALYSIS 

The  linear  aerodynamics  assumption  required  in  writing  closed  form  solutions  to  the 
differential  equations  have  placed  certain  restrictions  on  these  analysis  techniques.  These 
restrictions  have  been  recognized  by  many  researchers  and  various  quasi-nonlinear 
methods  have  been  proposed  by  several  individuals,  see  Refs.  6.7  to  6.10.  All  of  these 
quasi-methods  attempt  to  take  the  aerodynamic  coefficients  and  derivatives  obtained 
using  the  linear  theory  techniques  and  infer  from  those  estimates  what  the  nonlinear 
effects  are.  Murphy,  Refs.  6.1,  6.7,  and  6.10,  has  developed  methods  which  assume 
certain  forms  of  the  nonlinear  aerodynamics  and  these  will  be  described  in  the  following 
sections.  These  methods  are  easily  applied  by  the  analyst  and  have  been  shown  to 
produce  very  good  results.  It  is  not  our  purpose  herein  to  derive  all  of  these  nonlinear 
methods;  however,  we  will  attempt  to  provide  the  reader  with  the  general  processes 
involved  in  arriving  at  the  quasi-nonlinear  techniques. 

6.2.1  Drag 

The  determination  of  the  total  drag  coefficient,  Co,  was  discussed  in  Section  6.1.1.  As 
was  noted  in  that  section,  Cq  is  a  function  of  the  angle-of-attack  history  experienced  by 
the  model  as  it  traverses  the  instrumented  range.  Therefore,  the  total  drag  coefficient  is  an 
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average  value  corresponding  to  an  effective  angle  of  attack  experienced  during  the 
model’s  flight.  It  is  assumed  that  Cp  can  be  represented  as  a  quadratic  function  of  this 
effective  angle-of-attack,  or 


C0  -  C0O  +  C02  S 

6.38 

X2 

Then  from  Ref  6.7, 

=  [1/ ( X2  -  X,)]  J^^dx 

6.39 

Xl 

where,  5  =  a  +  P 

Cdo  =  zero  angle-of-attack  drag 
Cd2  =  quadratic  slope  term 


The  effective  angle-of-attack  squared  term,  S^,  is  normally  obtained  from  a  fit  of  the 
angular  motion  by  integrating  Eq.  6.21  for  6  and  accomplishing  the  summation  indicated 
by  Eq.  6.39.  It  should  be  recognized  that  the  same  portion  of  the  flight  used  to  determine 

Cd  must  also  be  used  to  accomplish  the  angular  motion  fit  and  hence  the  S  calculation. 

The  process  of  determining  the  zero  angle-of-attack  drag  coefficient,  C0O,  and  the 
quadratic  drag  term,  Cd2>  then  consists  of  plotting  the  determined  Cq  values  as  a  function 

of  the  respective  S  as  illustrated  in  the  sketch  below. 


As  can  be  inferred  from  the  above  sketch  this  process  requires  several  flights  of  the 
identical  configuration  at  each  test  condition.  Normally,  once  the  €^2  term  is  determined 
from  the  group  of  flights  at  each  test  condition,  the  C^o  term  is  computed  for  each  flight 
by  using  the  transpose  of  Eq.  6.38,  or 

t^DO  ~  ^0  -  C02  S  6.40 

Obviously  the  higher  the  angles-of-attack  experienced  by  the  projectile  during  flight  the 
more  critical  an  accurate  determination  of  C02  becomes.  If  C02  cannot  be  determined 
because  of  too  few  data  points  or  relatively  small  angles-of-attack  it  can  be  estimated 
using  Eq.  6.41  below  (see  derivation  in  Section  7.1). 

Cd2*Cx2  +  Cn„  6.41 
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As  shown  in  the  derivation  of  Section  7.1  this  relation  assumes  small  angle  theory  and 
that  both  Cd  and  Cx  can  be  represented  as  a  quadratic  expansion.  The  technique  outlined 
above  for  determining  C^o  and  C^yi  from  the  total  drag  measurements  has  proven  over  the 
years  to  be  quite  adequate  for  most  configurations,  see  Ref  6.11. 

6.2.2  Pitching  Moment  and  Normal  Force 

Again  what  is  required  is  to  determine  an  effective  angle-of-attack  for  the  pitching 
moment  derivative,  C^a  ^and  the  normal  force  derivative,  Cn„  ,  in  a  similar  manner  as 
was  accomplished  in  the  previous  section  for  the  drag  coefficient.  These  effective  angles- 
of-attack  will  correspond  to  the  polynomial  expansions  and  hence  are  a  frmction  of  the 
assumed  nonlinearities.  This  process  will  be  illustrated  herein  by  deriving  the  effective 
angle-of-attack  for  an  assumed  cubic  pitching  moment.  The  process  is  identical  for  a  fifth 
order  expansion  as  shown  in  Ref  6.7  and  6.10.  The  cubic  expansion  is  chosen  here  for 
the  pitching  moment  because  it  represents  the  expansion  which  is  most  frequently  used  to 
analyze  free-flight  data.  Even  though  nonlinearities  are  frequently  found  in  free-flight 
spark  range  data,  the  angles-of-attack  are  normally  below  20  degrees  and  higher  order 
nonlinearities  are  not  prevalent.  Nevertheless  analysts  should  recognize  that  higher  order 
expansions  may  be  required  and  utilize  the  corresponding  effective  angle-of-attack 
parameters  as  derived  in  Refs.  6.7  and  6.10. 

If  we  assume  that  the  pitching  moment  can  be  represented  as  a  cubic  expansion  of  the 
angle-of-attack,  then 

Cm  =  Cmo+  C„c,o  a  +  c„,„3'  6.42 

and  the  derivative  with  respect  to  a,  becomes, 

Cma  ~  Cmao  Cma3  ^  where,  C„jq3  =  3  6.43 

The  technique  requires  a  determination  of  the  effective  angle-of-attack  squared,  6^^, 
which  will  correspond  to  in  Eq.  6.43.  This  determination  is  shown  below. 

If  we  rewrite  Eq.  6.14  assuming  zero  damping  and  zero  trim,  we  see  that 
(p  +  ia)"  +  (-  iP)  (P  +  ia)'  -  (M)  (p  +  ia)  =  0 
With  the  assumed  cubic  expansion  of  ,  M  in  Eq.  6.44  becomes 

M  =  Mo  +  M2  6e^ 

where 

Mo  =  (  pA  /  2mad)  k^^ 
and 

M2  =  (pA/2m^d)kt'^C„„3 

The  solution  to  Eq.  6.44  can  be  written  similarly  to  Eq.  6.20,  or 

P  +  ia  =  Kp  exp(  i(j)p'  x)  +  exp(  i<f»N'x)  6.46 

The  derivatives  with  respect  to  x  can  therefore  be  written  as  shown  below: 

(P  +  ia  )'  =  Kp  i(t)p'  exp(i(|)p'  x)  +  i<|)N'exp(i(j)N'x)  6.47 


6.44 

6.45 
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and 


(p  +  ia  )''  =  -  Kp  ((|)p' exp(  i(t)p'  x)  -  exp(i<|)N'x)  6.48 


In  Ref.  6.7,  Murphy  defines  6^  as  =  (p  +  ia )  (p  +  ia )  and  accomplishing  the  indicated 
multiplication  using  Eq.  6.46  one  obtains, 

5^  =  Kn^  +  Kp^  +  KnKp  [exp  (jipOx  +  exp  i((t)p'-  ^^')x  ]  6.49 

If  Eqs.  6.46  through  6.49  are  substituted  into  Eq.  6.44  and  like  terms  collected  one  can 


arrive  at  the  lengthy  relationship  shown  below: 

Kn  (exp  i(t)N')[-  -  Mo  -  M2  (  Kn^  +  2Kp^)] 

+  Kp  (exp  i(t)p')[  -  +  P(t)p'  -  Mo  -  M2  (  2Kn^  +  Kp^)] 

-  M2  {  KN%[exp  i(2(t)N  -  (|)p)]  +  KnKp'  [exp  =  0  6.50 

The  third  term  containing  the  mixed  frequencies  in  Eq.  6.50  is  neglected;  and,  since  Kjg 
and  Kp  cannot  be  zero  the  quantities  within  the  brackets  must  be  zero.  Or 

-  (|)n'^  +  P^n'  =  Mo  +  M2  (  Kn^  +  2Kp^)  6.51 

-  (|)p'^  +  P  (t)p'  =  Mo  +  M2(  2Kn^  +  Kp^)  6.52 

If  we  define  +  2Kp^  and  =  2Kn^  +  Kp^  and  solve  both  Eqs.  6.5 1  and 

6.52  for  P  and  set  equal  to  each  other  we  can  show  that , 

=  Mo  +  M2  [(<t)p'  5e2^  -  <t>N'  Sei^)  /  ((jlp'  -  (jlN')]  6.53 

Since  Eq.  6.53  is  equivalent  to  Eq.  6.45, 6^  must  be  defined  as: 

Se'-(<t)p'Se2'-‘|)N'5el')/(<t)p'-‘|)N')  6.54 

2  2 

Finally,  if  the  definitions  of  5ei  and  5^2  are  substituted  into  Eq.  6.54  the  effective  angle- 
of-attack  squared  can  be  written  as  shown  below: 

5e'  =  Kn'  +  Kp'  +  [((^n'  Kn'  -  (t)p'  Kp')  /  ((Iin'  -  6.55 


What  Eq.  6.55  represents  is  an  equivalent  angle-of-attack  squared  parameter 
consistent  with  a  cubic  pitching  moment  expansion.  In  practice  how  this  would  be 
utilized  is  similar  to  the  drag  analysis  of  the  previous  section.  For  example  the  determined 
Ctna  values  as  obtained  by  fitting  Eq.  6.20  and  using  Eq.  6.28  are  plotted  as  a  fimction  of 
5e^  as  illustrated  in  the  sketch  below. 


The  midrange  values  for  Kn  ,Kp  ,  (|)n',  and  (j)p'(  if  the  variable  frequencies  are  utilized, 
see  Eq.  6.30)  are  used  in  computing  5^'.  By  using  the  midrange  values  we  are  accounting 
for  the  damping  in  the  system  thus  relaxing  somewhat  the  assumption  of  zero  damping 
when  we  derived  the  relationship  for  8/.  Again  it  is  obvious  that  several  flights  of  the 
same  configuration  at  similar  test  conditions  must  be  available  to  accurately  determine 
the  and  Cn,a3  terms  using  this  technique.  If,  once  the  vs  5^  plot  has  been 
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accomplished,  the  plotted  data  cannot  be  represented  by  a  linear  function  then  a  higher 
order  nonlinearity  must  be  assumed.  Murphy  heis  derived  these  various  higher  order 
effective  amplitudes  in  various  references,  see  for  example  Ref.  6.7. 


A  similar  process  is  again  applied  to  the  measured  normal  force  derivatives,  Cn„, 
data.  Although  not  derived  herein,  the  effective  amplitude  parameter  assuming  a  cubic 
polynomial  expansion  that  results  from  this  process  is  shown  below. 

6e/  =  Kn'  +  Kp'  +  [Kn^  Kp'((t)N'"  +  (|)p'"  )  /  (Kp^  (1)n''  +  Kn'  ^/)]  6.56 

Again  plotting  Cn^  vs  5es^  yields  both  CNao  and  Cnos  ,see  sketch  below. 


These  quasi-nonlinear  techniques  for  Cq,  (previous  section),  C^a  and  are 
relatively  straight  forward  and  easy  to  apply  once  the  aerodynamic  coefficients  have  been 
determined  and  the  effective  amplitude  parameters  calculated.  This  process  has  been 
shown  to  provide  good  results  with  the  primary  disadvantage  being  the  number  of  flights 
required  to  adequately  determine  both  the  zero  angle-of-attack  coefficients  and  their 
respective  nonlinear  terms.  Also,  when  large  nonlinearities  are  present  the  very  nature  of 
‘Linear  Theory’  means  that  the  fits  to  the  experimentally  measured  trajectories  will  not  be 
within  the  measuring  capability  of  the  facility.  These  fits  will  normally  result  in  residuals 
that  display  a  definite  sinusoidal  variation.  Hence,  these  poor  quality  fits  to  the  measured 
motions  are  not  conducive  to  establishing  a  positive  feeling  for  the  accuracy  of  the  final 
results.  Therefore,  even  though  the  final  results  may  be  quite  good  it  can  be  difficult  to 
appreciate  the  quality  of  results  in  light  of  the  poor  fits  to  the  measured  motions. 

6.2.3  Damping  and  Magnus 

As  defined  previously,  the  damping  derivatives  being  extracted  from  free-flight 
trajectory  measurements  are  in  actuality  the  combination  of  (€„,<,  +  Cn,^).  Analysts  have 
become  lax  in  their  identification  of  these  derivatives  and  quite  frequently  refer  to  the 
combination  as  Cn,q  only.  This  simplification  is  often  used  in  writing  the  equations  of 
motion  and  in  tabulating  and  presenting  results.  This  abbreviated  nomenclature  was 
used  in  Section  2.3  and  will  be  used  again  in  this  section  but  with  the  understanding  that 
it  refers  to  the  combined  damping-in-pitch  derivatives. 
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Murphy,  Ref.  6.1,  again  derives  an  effective  angle-of-attack  associated  with  the 
determined  C^q  values  similar  to  the  process  previously  discussed.  This  expansion  has 
the  form  shown  below: 

Cmq  =  ^mqo  +  C„q2  A  6.57 

where 

A  =  -(  ly  /  Ix)  [(k'  +  W)  /  -  W)]  (  Kn'  -  W) 

Although  A  has  the  units  of  degrees^  (or  rad^)  it  can  take  on  negative  values  depending  on 
the  magnitudes  of  the  nutational  and  precessional  vectors.  Because  negative  values  of  A 
can  be  disconcerting  to  the  analyst  an  alternate  approach  is  frequently  used  to  evaluate 
potential  nonlinear  effects  associated  with  Cn,q.  This  alternate  approach  is  also  used  for 
the  Magnus  moment  derivative,  C„p„ ,  as  well.  The  effective  amplitudes  squared  for  the 
nutational  and  precessional  damping  rates  turns  out  to  be  5^1  and  5e2  ,  as  previously 
defined.  Or, 

6e2'sKN'  +  2Kp' 

5e,'  =  2KN'  +  Kp' 

The  process  then  becomes  plotting  the  nutational  damping  rate  vs  d^2  and  the 
precessional  damping  rate  vs  8gi^  as  shown  in  the  sketch  below. 


The  two  zero  angle-of-attack  damping  rates,  ppo  and  are  then  substituted  into  Eq. 
6.26  to  calculate  (Cn,q+Cn,^i)o  and  into  Eq.  6.29  to  calculate  Cnp„o  •  This  process  eliminates 
the  concern  associated  with  the  sign  of  A  as  previously  mentioned;  but,  does  not  easily 
lend  itself  to  determining  the  higher  order  nonlinear  coefficients. 


It  may  appear  to  the  reader  that  these  quasi-nonlinear  techniques  are  somewhat 
desperate  attempts  at  eliminating  the  severe  assumptions  of  linear  aerodynamics  in  the 
derivations  of  the  equations  of  motion.  Although  there  may  be  some  truth  in  this,  these 
quasi-nonlinear  techniques  have  proven  to  be  extremely  useful  over  the  years  and  have 
provided  very  good  estimates  for  both  the  zero  angle-of-attack  coefficients  as  well  as  the 
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nonlinear  terms.  A  good  example  of  the  application  of  these  techniques  are  provided  in 
Ref  6.12.  As  mentioned  previously  the  primary  disadvantage  has  been  the  large  number 
of  flights  required  to  accurately  determine  the  nonlinearities. 

6.3  NUMERICAL  INTEGRATION 

Linear  Theory  and  the  associated  quasi-nonlinear  techniques  discussed  in  the  previous 
sections  were  developed  in  the  manner  discussed  primarily  because  of  a  lack  of  computer 
power  during  the  period  of  1940  through  the  60’s.  Because  of  this  lack  of  computer 
power  closed  form  solutions  of  the  differential  equations  were  required  in  order  to  extract 
the  aerodynamic  coefficients  and  derivatives.  This  changed  considerably  in  1970  when 
Chapman  and  Kirk  published  a  method  of  data  analysis  they  called  parametric 
differentiation*^  Since  the  differential  equations  were  used  directly  by  numerically 
integrating  the  equations  and  comparing  the  resulting  motions  with  the  experimental 
measurements,  many  of  the  assumptions  required  by  the  linear  theory  techniques  were  no 
longer  necessary.  This  numerical  integration  (NI)  method  has  developed  to  the  point 
where  typically  all  of  the  differential  equations  describing  the  motions  in  the  x,y,z 
directions  along  with  the  angular  motions, 0  ,v}/,  and  (j)  are  simultaneously  fitted  for 
multiple  flights.  These  differential  equations,  both  fixed  plane  and  body  fixed,  have  been 
presented  in  Section  2.2.2  and  will  not  be  repeated  here. 

This  NI  technique  also  permits  the  aerodynamic  coefficients  and  derivatives  to  be 
functions  of  the  angle-of-attack,  Mach  number,  and  roll  orientation  or  roll  rate,  see  for 
example  the  aerodynamic  expansions  shown  in  Section  2.3.  In  fact,  the  primary 
restriction  on  this  technique  is  the  analyst’s  ability  and  ingenuity  in  mathematically 
modeling  the  aerodynamic  variations.  Since  the  assumptions  of  linear  aerodynamics  with 
angle-of-attack  are  no  longer  required  this  technique  is  frequently  referred  to  as  nonlinear 
analysis.  This  technique  has  been  described  in  great  detail  in  various  references,  see  Refs. 
2.2  to  2.5  and  6.13;  therefore,  only  a  simple  explanatory  example  will  be  described 
herein. 

6.3.1  Drag  Example 

We  will  rewrite  Eq.  6. 1  as  shown  below: 

jc=-(pA/2ma)(x)^CD  6.58 

Now  if  we  expand  the  total  drag  coefficient  to  be  a  function  of  angle-of-attack  and 
velocity  as, 

Cd  “  Cqo  +  Cd2  £  Cdv  (^  ■  Vref)  6.59 

we  obtain, 

x=  K( i: )  [  Cpo  C[)2  £  Cqv  ~  Vref)  ]  6.60 

where  K  =  -  (  p  A  / 

8^  =  sin^  a  ,( see  Eq.  2. 1 4 ) 

Vref  ~  reference  velocity 
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The  drag  expansion  shown  in  Eq.  6.59  permits  the  total  drag  to  be  a  quadratic  function  of 
the  angle-of-attack  with  a  linear  variation  of  velocity  about  a  reference  value. 

The  problem  now  is  to  determine  the  unknown  Cdq?  Cd2>  Cpy  coefficients  of  the 
expansion  by  matching  the  theoretical  trajectory  as  defined  by  Eq.  6.60  to  the 
experimentally  measured  time,  t,  and  distance,  x,  data.  It  is  also  assumed  that  the  angle- 
of-attack  history  is  available.  If  we  define  the  partial  derivatives,  P; ,  as  follows: 

?i=  X  /<?Cdo  ,  P2  =  X  /  ,  P3  =  X  /  ,  and 

P4S  (^X/  <^Xo  ,P5=  <?x/ 6.61 

We  will  also  assume  that  the  order  of  differentiation  can  be  reversed,  or 
IdC;  (d^x  /  dt^)  =  d^  /  dt^(^  x  /^Cj  )=  d  x  IdC;  =  P , 

and  similarly, 

d  /^Cj  (dx /dt)  =  d /dt  (<^  X /(^Cj )  =^x/^5’Ci  =  P\ 

Using  the  above  relationships  in  conjunction  with  Eq.  6.60  we  can  write: 

P^  =  d  X  /^?Cdo  =  K;c[2Cdo  +  2Cd2  8^  +  Cdv(3x  -  2Vref)]  (^^/<^Cdo)  + 

After  expanding,  collecting  like  terms,  and  simplifying  we  arrive  at: 

P 1  =  K(io)^  [( 1  3  Cdvo  P  1)  +  2  P  i(Cdo  ^02  ^  ~  Cdv  Vref)o  ]  6.62 

Where  the  subscripted,  0 ,  terms  imply  that  the  quantity  was  evaluated  with  the  given 
coefficients,  either  initial  guesses  or  corrected  values,  after  each  iteration.  Likewise  the 
remaining  partials  can  be  derived  and  the  resulting  relations  are  shown  below: 

P2  “  K( X 0)^  [(s^  +  3  Cdvo  P 2)  2  P 2(Cdo  +  Cd2  s  *  Cqv  Vref)o  1  6.63 

P3  =  K( i 0)^  {[(  X  -  Vr£f)o  3  Cqvo  P3]  +  2  P 3(Cdo  +  C02  ^  -  Cqv  Vref)o  }  6.64 
P4  =  K(xo)'  [( 3  Covo  ^4)  +  2  P4(Cdo  +  Cd2  s'  -  Cdv  Vref)o  ]  6.65 

P  5  =  K(  x  0)^  [(  3  Cdvo  P  5)  +  2  P  5(Cdo  +  Cd2  s  -  C^v  Vref)o  ]  6.66 

Equation  6.60  is  numerically  integrated  utilizing  the  initial  guesses  for  the  unknown 
coefficients  of  the  drag  expansion  and  the  estimated  initial  conditions  (x  q  and  x  q). 
Equations  6.62  through  6.66  are  also  numerically  integrated  such  that  the  partial 
derivatives  with  respect  to  each  of  imknown  coefficients  are  evaluated.  The  method  of 
differential  corrections  is  then  used  to  obtain  the  corrections  to  each  of  the  initial  guesses. 
This  method  consists  of  expanding  the  calculated  value  of  Xj  <,3]  about  a  given  set  of 
coefficients  in  a  Taylor  series,  or 

N 

Xi  cal  =  Xi  caio  X  ^ ^  ‘  terms  6.67 

j  =  i 
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Note  that  the  higher  order  terms  are  neglected,  N  is  the  total  number  of  unknowns,  and 
that  the  subscript,  0,  again  indicates  that  the  term  is  evaluated  using  the  given  set  of 
values,  either  initial  guesses  or  updated  values  after  each  iteration. 

The  sum  of  the  squares  of  the  residuals  can  be  written  as, 

NP 

RSQ=X  6.68 

i  =  l 

The  subscript  i  denotes  the  i-th  measurement  and  NP  the  total  number  of  measurements 
If  we  substitute  Eq.  6.67  into  6.68,  taking  the  derivative  of  RSQ  with  respect  to  each  of 
the  unknown  coefficients(  C|j )  and  setting  equal  to  zero,  we  can  arrive  at  the  following 


equation: 

C  =  (Ajj)-'B» 

where. 

NP 

Ajk=  E  <?X|/^CjL  {<?x,/.?Ci„) 

i  =  1 

and. 

NP 

Bk=  X  (Xi-Xi,calo)(^Xi/^7Ck) 
i=l 

The  C  matrix  in  Eq.  6.69  represents  the  corrections  to  be  added  to  each  of  the  unknowns 
to  be  determined,  including  the  initial  conditions  x  g  and  x  g.  Once  these  corrections  are 
added  to  the  initial  estimates  then  these  new  estimates  are  substituted  back  into  Eq.  6.60 
and  6.62  through  6.66  and  the  process  repeated  until  convergence  is  achieved. 
Convergence  is  defined  as  when  the  change  in  the  sum  of  the  residuals  squared  is  less 
than  a  small  predetermined  value.  This  small  predetermined  value  is  set  by  the  analyst 
depending  on  the  precision  deemed  necessary. 

Several  aerodynamic  modeling  investigations  have  been  accomplished  over  the  years 
using  the  very  routine  described  in  this  section.  These  investigations  have  included  both 
drag  and  roll  expansions  using  various  polynomials  and  splines,  see  Refs.  6.14  and  6.15. 


6.3.2  Maximum  Likelihood  Method 

In  the  simple  drag  example  shown  in  the  previous  section  the  data  correlation 
technique  was  basic  least  squares.  When  simultaneously  fitting  multiple  sets  of  equations 
of  motion  such  as  those  listed  in  Section  II  to  experimental  data  containing  different 
levels  of  measurement  noise  basic  least  squares  cannot  be  used.  There  does  exist  various 
correlation  techniques  which  can  handle  this  situation,  i.e,  weighted  least  squares,  WLS, 
and  the  maximum  likelihood  method,  MLM,  along  with  others.  An  error  matrix 
containing  the  known  measurement  noise  is  incorporated  into  the  least  squares  theory  to 
account  for  the  different  levels  of  measurement  noise.  Whereas,  the  MLM  can  determine 
the  error  matrix  along  with  the  unknown  parameters  in  the  equations  of  motion.  The 
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correlation  technique  we  have  chosen  for  fitting  the  fixed  plane  and  body  fixed  equations 
of  motion  shown  in  Section  II  is  a  hybrid  of  both  WLS  and  MLM.  This  techmque  is 
referred  to  as  a  hybrid  because  although  the  likelihood  function  is  maximized  the  error 
matrix  is  assumed  known.  This  simplification  to  the  standard  MLM  approach  is  justified 
because  after  analyzing  numerous  ffee-flight  trajectories  the  facility’s  measurement  error 
levels  are  known  to  a  high  degree.  This  modified  MLM  is  an  iterative  procedure  which 
adjusts  the  unknown  aerodynamic  coefficients  until  the  theoretical  equations  of  motion 
match  the  experimental  data. 

The  mathematical  treatment  is  a  sensitivity  analysis  where  the  partial  differential 
equations  form  a  basis  of  quasi-linearization.  These  are  derived  by  differentiating  the 
equations  of  motion  with  respect  to  each  of  the  unknown  coefficients  and  derivatives  as 
was  done  in  the  preceding  simple  drag  example.  Again  the  sensitivity  equations  are 
integrated  in  parallel  with  the  equations  of  motion  to  yield  the  sensitivity  coefficients 
(partial  derivatives)  which  reflect  the  sensitivity  of  the  computed  solution  with  respect  to 
each  of  the  unknown  parameters.  The  differential  corrections  procedure  in  conjunction 
with  the  Taylor  series  expansion  is  again  used  to  form  the  relationship  between  residuals, 
corrections  to  the  unknowns,  and  partial  derivatives.  The  likelihood  function  is  then 
defined  as,  n 

L  =  [1  /  ( I S  r^]  exp[(  -1/2)  X  R  ]  6.70 

where, 

N  =  number  of  data  points 
S  =  covariance  matrix  of  measurement  noise 
R  =  residuals  after  corrections  have  been  made 

The  application  of  this  likelihood  function  in  the  trajectory  matching  process 
eliminates  the  inherent  assumption  in  least  squares  that  the  magnitude  of  the 
measurement  noise  must  be  consistent  between  dynamic  parameters,  irrespective  of  imits. 
This  is  achieved  through  the  covariance  matrix  which  contains  the  information 
concerning  the  measurement  precision.  The  likelihood  function  also  has  the  property  of 
asymptotically  approaching  the  solution  as  compared  to  least  squares  which  approaches 
the  solution  parabotically.  A  good  discussion  of  data  correlation  techniques  can  be  foimd 
in  Ref.  6.16  and  many  texts. 

6.4  DATA  ANALYSIS  SOFTWARE  SYSTEM 

The  software  system  developed  for  analyzing  ffee-flight  trajectory  data  and  extracting 
the  aerodynamic  coefficients  and  derivatives  is  the  Aeroballistic  Research  Facility  Data 
Analysis  System  (  ARFDAS).  This  system  is  illustrated  in  the  Fig.  6.1  below  and 
described  in  detail  in  Refs.  2.2  to  2.4. 
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Figure  6.1  Aeroballistic  Research  Facility  Data  Analysis  System  ( ARFDAS ) 

This  software  system  uses  the  experimental  trajectories,  CADRA  output  -  see  Section  5.4, 
along  with  the  model’s  physical  characteristics  and  the  atmospheric  conditions  in  the 
facility  and  accomplishes  a  preliminary  screening  of  the  raw  trajectory  data.  The  system 
then  does  linear  theory  fits  consistent  with  the  process  as  described  in  Section  6.1  and 
provides  these  results  to  the  appropriate  NI  routine  for  a  nonlinear  analysis.  Once  each  of 
the  experimental  trajectories  have  been  analyzed  separately  this  routine  can  then  be  used 
to  fit  multiple  trajectories. 

Obviously  only  similar  configurations  at  nearly  the  same  test  conditions  should  be 
simultaneously  analyzed.  The  techniques  as  described  herein  can  be  a  complex  and 
lengthy  process  requiring  considerable  computer  power.  This  is  especially  true  when 
multiple  sets  of  equations  are  utilized  such  as  the  fixed  plane  or  body  fixed  equations  of 
motion  as  outlined  in  Section  2.2.2.  Obviously,  not  all  of  the  associated  coefficient 
expansion  terms,  as  shown  in  Section  2.3,  are  determined  from  each  experimentally 
measured  trajectory.  The  analyst  requires  considerable  experience  in  determining  which 
of  the  potential  aerodynamic  expansions  are  likely  to  affect  the  motions  of  the 
configuration  being  tested.  Usually  this  requires  fitting  the  trajectories  several  times  with 
different  combinations  of  the  unknowns  being  examined.  Although  at  first  glance  this 
may  appear  to  be  an  impossible  task,  the  selection  of  the  correct  combination  of 
unknowns  becomes  quite  manageable  and  logical  with  experience. 
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Figure  6.3  Example  of  Aerodynamic  Results 
(Blunt  Planetary  Probe) 


Examples  of  typical  NI  curve  fits  to  experimental  angular  motions  are  shown  in  Fig. 
6.2.  These  fits  were  accomplished  using  ARFDAS  and  demonstrate  both  body  fixed  and 
fixed  plane  analyses.  Many  such  fits  for  the  angular  as  well  as  roll  and  swerve  motions 
are  presented  in  various  papers  and  data  reports,  see  Appendix  6.  Examples  of  some 
aerodynamic  results  from  a  planetary  probe  configuration  are  also  shown  in  Fig.  6.3. 
These  results  were  obtained  fi’om  analyses  of  single  trajectories  as  well  as  multiple  flights 
and  are  shown  here  only  for  illustrative  purposes. 
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VII.  AERODYNAMIC  ANALYSIS 


Once  the  aerodynamic  coefficients  and  derivatives  have  been  determined  from  the 
experimental  trajectories  additional  analysis  is  normally  required  in  order  to  determine 
the  missile’s  center-of-pressure,  gyroscopic  stability  factor,  dynamic  stability  and  other 
aerodynamic  parameters.  Also  in  order  to  compare  these  results  with  data  obtained  from 
previous  tests,  the  moment  data  may  require  shifting  to  a  different  reference  location.  The 
purpose  of  this  section  is  to  develop  some  of  the  relationships  used  for  these  purposes. 

7.1  AERODYNAMIC  DRAG 


The  relationships  between  the  aerodynamic  drag,  axial  force,  lift  and  normal  force  are 
illustrated  in  the  sketch  below. 


If  we  assume  that  both  Cq  and  Cx  can  be  expanded  as  a  quadratic  polynomial  in  a,  then 
Eq.  7.1  becomes, 

Cdo  +  Cd2  =  ( Cxo  +  Cx2  cc^)  cos  a  +CNaasina  7.3 

Now  if  we  apply  small  angle  theory,  cos  awl  and  sin  a  w  a ,  and  recognizing  that  Cxo  is 
equal  to  Cdq  then  Eq.  7.3  can  be  revmtten  and  simplified  to  obtain, 

Cd2  ~  Cx2  +  7.4 

This  relationship  is  occasionally  used  to  estimate  Cd2  when  insufficient  data  are  available 
to  determine  it  from  other  sources.  The  Cx2  term  is  normally  small  ( »!  )  compared  to 
Cnoc  which  can  take  on  values  of  2  to  3/rad  for  spin  stabilized  projectiles,  cones,  and 
other  non  finned  missile  configurations  to  values  as  high  as  8  -12  for  finned 
configurations.  Another  common  relationship  can  be  obtained  from  Eq.  7.2.  Again 
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assuming  small  angle  theory  and  taking  the  derivative  with  respect  to  a,  one  obtains  the 
relations  shown: 


or 


^  Na  *  ^  La 
CxO”Cdo  * 


+  Cxo 
CNa-C 


La 


7.5 

7.6 


7.2  STABILITY 
7.2.1  Static 

The  static  stability  criteria  for  any  system  is  defined  as  the  tendency  to  restore  itself 
to  its  original  state  once  a  disturbance  has  been  input  to  the  system.  The  classic  examples 
are  the  ball  on  a  hill  and  the  ball  in  the  valley,  see  below. 


a.  Statically  Unstable  b.  Statically  Stable 

Here  if  the  ball  in  sketch  ‘a’  above  is  disturbed  from  its  point  of  rest  the  ball  continues  to 
deviate  from  its  initial  position;  hence,  the  system  is  statically  unstable.  However,  when 
the  ball  in  sketch  ‘b’  is  disturbed  it  will  tend  to  return  to  its  initial  position  and  after 
oscillating  back  and  forth  will  again  come  to  rest.  This  system  is  defined  as  statically 
stable.  This  analogy  when  applied  to  an  object  in  flight  implies  that  the  forces  and 
moments  acting  on  that  body  tend  to  return  the  angle-of-attack  to  the  trim  angle  once  the 
object  has  been  disturbed.  There  are  two  distinct  techniques  by  which  static  stability  can 
be  provided  to  a  free-flight  object.  These  are  typically  referred  to  as  aerodynamic  and 
gyroscopic  stability. 

7.2.1. 1  Aerodynamic 

Aerodynamic  static  stability  is  illustrated  in  the  sketch  below. 
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This  sketch  indicates  that  at  an  angle-of-attack  the  resulting  normal  force  creates  a 
moment  about  the  missile’s  center-of-gravity  which  tends  to  return  the  angle-of-attack 
back  to  zero.  Therefore  the  missile’s  static  stability  can  be  decreased  or  increased  by 
shifting  the  center-of-gravity  (forward  shift  provides  more  stability)  or  by  changing  the 
magnitude  and/or  the  location  of  the  resulting  normal  force.  This  alteration  of  the  normal 
force  is  exactly  what  the  addition  of  fins  provide  to  the  static  stability  of  the  missile.  Note 
that  if  the  fins  were  forward  of  the  center-of-gravity  (i.e.  canards)  they  would  provide  a 
destabilizing  moment  to  the  missile. 

The  derivative  of  the  pitching  moment  with  respect  to  a,  » is  often  referred  to  as 
the  static  stability  derivative  because  the  sign  of  the  derivative  defines  whether  or  not  the 
missile  possesses  aerodynamic  static  stability.  If  the  sign  of  is  negative  the  missile 

is  statically  stable.  This  is  true  for  symmetrical  as  well  as  asymmetric  configurations; 
however,  the  normal  convention  used  by  aircraft  designers  for  the  yawing  moment 
derivative,  Cpp ,  is  reversed.  Therefore  if  the  sign  of  Cpp  is  positive  the  aircraft  possesses 
static  stability  in  the  yaw  plane.  Note  that  this  arbitrary  sign  change  is  not  reflected  in  the 
aerodynamic  expansions  for  the  body  fixed  equations  of  motion  as  shown  in  Section 
2.3.1.  and  both  C^a  and  Cnp  ,  as  used  herein,  require  negative  values  for  static  stability  to 
exist  in  the  pitch  and  yaw  planes. 

7.2.1.2  Gyroscopic 

Missiles  or  projectiles  which  are  inherently  aerodynamically  unstable,  i.e.  the  normal 
force  creates  a  destabilizing  moment,  can  be  stabilized  with  spin.  For  this  condition  the 
spin  rate  must  be  high  enough  to  develop  a  gyroscopic  moment  which  overcomes  the 
aerodynamic  instability  and  the  missile  is  said  to  be  gyroscopically  stable.  This  is 
typically  the  case  for  most  gun  launched  projectiles  (handguns,  rifles,  cannons,  etc.) 
where  the  rifling  of  the  barrel  provides  the  required  spin  to  the  projectile. 

In  describing  this  condition  ballisticians  have  developed  a  gyroscopic  stability  factor, 
Sg  ,  which  is  obtained  from  the  roots  of  the  linear  theory  equation  of  motion,  Eq.  6.21 . 
The  boundary  of  the  stability  limits  is  defined  as  the  condition  where  both  and  pp  are 
zero.  Therefore  from  Eqs.  6.15  and  6.26  it  can  be  seen  that, 

Pn  +  Pp  =  -  H  =  0 


Likewise  from  Eqs.  6.17  and  6.29, 

(  ‘I^pVn  /  (*f*p'  ‘I’nO  -  -  t  =  0 

Substituting  zero  for  both  H  and  T  into  Eq.  6.21  yields  Eq.  7.7  below: 

(|)p'  =  (1/2)  P  ±  (1/2)  ( -  4M)^^^  7.7 

If  we  now  substitute  Eq.  6.18  into  Eq.  7.7  above  we  can  obtain, 

(j)N',  (t)p'  =  (I,  p  /  2  ly)  +  (1/2)  [  ( I,  p  /  ly)  2  -  4M  7.8 
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We  will  now  define  the  gyroscopic  stability  factor  as, 


7.9 


Using  this  definition  Eq.  7.8  can  be  rewritten  as, 

'l'N'.'t'p'  =  (IxP/2I,)[l  ±  (l-l/sy'^l  7.10 

Eq.  7.10  demonstrates  that  if  the  frequencies  are  to  be  real  numbers,  Sg ,  must  always 
be  greater  than  1 .  Hence  if  a  missile  is  to  be  gyroscopically  stable  Sg  >  1 .  The  form  of  Sg 
as  shown  in  Eq.  7.9  is  not  a  very  useful  relationship  for  computing;  therefore,  if  we 
substitute  Eq.  6.16  into  Eq.  7.9  we  can  arrive  at  the  equation  shown  below: 

S,  =  (2l/p")/(I,pjid’C„„)  7.11 

Note  that  kj'  in  Eq.  6. 16  can  be  written  as, 

k,'^  =  mad^/Iy  7.12 

Equation  7.1 1  is  the  relationship  typically  shown  or  used  when  discussing  the  level  of  the 
gyroscopic  stability  of  a  spin  stabilized  configuration. 


7.2.2  Dynamic 

Dynamic  stability  is  defined  as  the  condition  where  a  system  is  perturbed  and  the 
ensuing  oscillatory  motion  has  a  tendency  to  either  decrease(damp)  or  increase(undamp). 
Note  that  this  definition  assumes  that  static  stability  is  present,  otherwise  the  oscillatory 
motion  would  not  occur.  This  static  stability  criteria  can  be  provided  by  either  the 
aerodynamic  or  gyroscopic  moments  as  previously  defined.  Some  analysts  and  test 
engineers  feel  that  dynamic  stability  is  a  somewhat  complicated  system  parameter  to 
describe  and  define.  This  section  will  attempt  to  provide  a  simple  derivation  of  the 
dynamic  stability  criteria  and  hopefully  put  at  ease  those  who  feel  uncomfortable  about 
this  very  important  aerodynamic  effect. 

Again  we  will  begin  by  using  some  previously  derived  relations,  see  the  previous 
section  and/or  Section  6.1.3. 


Pn  +  fip  “  •  H 

7.13 

(  ‘I^pVn  Pp  “  ■  PT 

7.14 

(<|)p'+(|)n')  =  P 

7.15 

Solving  Eqs.  7.13  and  7.14  for  pp  we  arrive  at, 

PN  =  (PT-(t,N'H)/((j,N'-(t)p') 
and  Pp  =  ( PT  -  (t)p'  H  )  /  ((|)p'  -  (|)i,,') 

Using  Eq.  7.15  the  two  equations  above  can  be  written  as  shown  below: 

Pn  =  -(<1)n'H-PT)/(2V-P)  7.16 

and  Pp  =  -  ( (|)p' H  -  PT)  /  (2«|)p' -  P)  7.17 
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Equation  7.17  can  now  be  expanded  by  adding  and  subtracting  a  like  term,  or 

Pp  =  -  [( (|)p'  H  -  PT)  /  (2<|)p'  -  P)]  +  [(1/2)  HP  /  (2W  -  P)]  -  [(1/2)  HP  /  (2^  -  P)] 


After  simplifying  we  arrive  at, 

Pp  =  -(l/2)[H-P(2T-H)/(2(t)p'-P)] 

or,  Pp  =  -  (l/2)[  H  -  P(2T  -  H  )  /  ((|)p'  -  7.18 

Likewise  for  Eq.  7.16,  or 

Pn  =  -  (l/2)[  H  +  P(2T  -  H  )  /  (V  -  (t.N')]  7.19 


Equations  7.18  and  7.19  shows  that  p^  +  pp  =  -  H ,  which  agrees  with  Eq.  7.13.  Also  that 
for  zero  spin,  P=0, 

Pn-Pp  =  -(1/2)H 


As  the  spin  rate  increases 


Pp  =  -(1/2)(H-  A) 
Pn  =  -(1/2)(H  +  A) 


where 


A  =  P(2T-H)/((|)p'-<^n') 


7.20 


Now  if  we  substitute  Eq.  7.10  for  (t)p'  and  (1)^'  into  Eq.  7.20  we  can  show  that, 

A  =  -  (2T  -  H  )  /  ( 1  - 1  /  Sg  )‘'^ 
or  A  =  -  H(  2T/H  - 1 )  /  ( 1  - 1  /  Sg  )^'^ 


The  above  derivation  illustrates  that  nutational  and  processional  damping  rates,  p^ 
and  Pp  ,  can  be  represented  as  functions  of  spin  (P)  as  shown  in  the  sketch  below. 


+ 


103 


If  the  missile  is  always  dynamically  stable  then  both  ^nd  |a,p  <  0  and  from  the  sketch  it 
is  obvious  that  A/2  <  H/2  ,  or  A  <  H.  Therefore; 

A  =  -H(2T/H-1)/(1  -1/Sg)'^^  <  H 
or 

(2T/H-l)/(l-l/Sg)‘^^  >1  7.21 

Now  if  we  define  the  dynamic  stability  parameter  ,S(j ,  as  shown  below: 

Sds2T/H  7.22 

Substituting  this  definition  back  into  Eq.  7.21  yields, 

(Sd-lf  >  1-1/Sg 

or  Sd  >  1  ±  (1  - 1  /  Sg  7.23 


Equation  7.23  effectively  maps  the  dynamic  stability  bounds  of  any  missile 
configuration  as  illustrated  in  Fig.  7. 1 


The  area  above  the  Sd  axis  represents  all  gyroscopically  stabilized  projectiles,  i.e.  Sg  >  1 ; 
likewise,  the  area  below  the  Sd  axis  represents  all  other  missile  configurations.  Once  the 
missile’s  aerodynamics  are  known,  either  from  testing  or  predictions,  then  Sg  can  be 
computed  using  Eq.  7. 1 1  and  Sd  obtained  fi’om  the  relationship  shown  below. 
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Sd  =  2(  Cn„  -  Cd  +  k;'  C„p J  /  [  CNa  -  2Cd  -  +  C^s)]  7.24 

Equation  7.24  was  obtained  by  substituting  Eqs.  6.15  and  6.17  into  Eq.  7.22. 

Plotting  1/Sg  vs  Sd  as  indicated  by  Fig.  7. 1  provides  the  analyst  or  missile  designer 
with  an  estimate  of  the  configuration’s  stability  characteristics.  This  figure  also  illustrates 
where  a  precessional  or  nutational  instability  (  p.p  or  >  0 )  would  fall  on  the  stability 
criteria  chart.  Obviously  if  the  aerodynamics  are  nonlinear  with  angle-of-attack  this 
evaluation  may  be  required  over  a  range  of  amplitudes.  It  would  not  be  unusual  to  find 
that  a  missile  is  dynamically  stable  for  some  angles-of-attack  and  unstable  at  others. 
Frequently  these  conditions  lead  to  angular  motion  patterns  that  develop  into  limit  cycles. 
Also,  occasionally  analysts  use  terms  such  as  ‘Magnus  instability’  and  this  would  simply 
be  an  instability  caused  by  a  Magnus  moment  derivative,  Cnp^  ,  of  such  a  magnitude  as  to 
cause  Sd  to  be  outside  the  stability  bounds  as  defined  by  Eq.  7.24.  Since  Cnpa  is  an  out  of 
plane  moment  ( perpendicular  to  the  plane  of  the  angle-of-attack)  this  moment  is  quite 
often  the  cause  of  such  dynamic  instabilities.  Some  configurations  may  exhibit  out  of 
plane  moments  which  are  not  functions  of  spin,  i.e.  wrap  around  fin  configurations^ 

For  these  configurations  the  Cnp^  term  in  Eq.  7.24  is  replaced  with  a  term  and  this 
out  of  plane  moment  can  also  be  a  source  of  dynamic  problems,  see  Ref.  7. 1 . 

7.3  OTHER  USEFUL  RELATIONS 
7.3.1  Center-of-Pressure 

The  definition  of  the  center-of-pressure  is  the  point  on  the  body  where  the  resultant 
aerodynamic  forces  can  be  assumed  to  act.  Therefore  if  we  repeat  the  sketch  shown  in 
Section  7.2. 1.1  the  normal  force,  N,  would  be  acting  through  the  center-of-pressure,  see 
the  sketch  below. 


The  center-of  pressure  can  therefore  also  be  defined  as  the  point  where  the  resultant 
aerodynamic  moment  is  zero.  In  order  to  compute  this  point  from  the  determined  and 
Cn,„  data  we  start  by  writing  the  moment  transfer  equation,  or 
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niA  =  niB  +  F  X 


7.25 


This  relation  says  that  any  moment  about  point  A  is  equal  to  the  moment  about  any  other 
point  B  plus  the  force  times  distance  between  points  A  and  B.  Therefore  if  we  let  point  A 
be  the  center-of-gravity(cg) ,  point  B  the  center-of  -pressure(cp),  N  the  normal  force,  and 
X  the  distance  between  the  eg  and  cp,  then  Eq.  7.25  can  be  rewritten  as  shown: 

m^g  =  mep  +  N  (x  /  L)  L 

Expanding  x  as  x^g  -  x^p ,  substituting  the  definitions  for  the  force  and  moment 
coefficients,  and  note  that  mep  =  0,  the  above  equation  can  be  written  as: 

x,p/L  =  (x,g/L)-(C„/CN)(d/L)  7.26 

Equation  7.26  is  the  relationship  used  to  compute  the  center-of-pressure  location  from 
force  and  moment  coefficient  data  as  obtained  in  a  wind  tunnel;  where,  the  (x^g  /  L)  term 
is  the  moment  reference  location  for  that  wind  tunnel  data.  If  the  force  and  moment 
derivatives  had  been  obtained  instead  of  the  coefficients,  i.e.  ffee-fiight  data,  then  Eq. 

7.26  could  be  similarly  written  as: 

x,p  /  L  =  (x,g/  L  )  -  ( /  CnJC  d  /  L)  7.27 


It  should  be  noted  that  the  distance  between  the  center-of-gravity  and  the  center-of- 
pressure,  x^g  -  X(,p ,  is  referred  to  as  the  missile’s  static  margin. 

7.3.2  Moment  Reference  Shift 

In  contrast  to  the  aerodynamic  forces,  the  moments  are  a  function  of  the  moment 
reference  position.  For  fi’ee-flight  data  this  is  normally  the  model’s  eg  location,  at  least 
during  the  trajectory  analysis  process.  However,  for  wind  tunnel  data  the  moment 
reference  location  may  be  any  position  which  is  convenient  for  the  model  /  balance 
installation.  Therefore  similar  configurations  may  be  tested  at  various  times  and  facilities 
using  different  moment  reference  locations.  In  order  to  compare  the  results  from  these 
different  tests  the  moment  data  must  be  shifted  to  a  common  reference  location.  Typically 
this  common  reference  location  might  be  the  actual  eg  location  for  the  full  scale  vehicle. 

In  order  to  accomplish  this  shift  we  will  again  start  with  Eq.  7.25.  Substituting  the 
force  and  moment  definitions  into  this  equation  the  following  relationship  can  be  written: 

C„.B  =  C„A-CN[(XA/L)-(XB/L)](L/d) 

Again  is  the  moment  about  the  eg;  however,  now  C^b  is  the  moment  about  an 
arbitrary  common  reference  location,  then 

C^mREF  “ Cmcg - Cn  [( Xjg / L) - ( Xrbp /  L)](L/d)  7.28 

The  corresponding  relationship  for  the  derivatives  is  therefore, 

CmaREF  ~  Cn^^gg  -  [(  Xgg  /  L)  -  (  Xref  /  L)](L  /  d)  7.29 
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Equations  7.28  and  7.29  are  the  relationships  used  to  shift  the  moment  data  to  any  other 
reference  location  and  thereby  compare  results  obtained  from  different  facilities  or  at 
different  times. 


7.3.3  Normal  Force  Relation 

As  was  shown  in  Section  6.1.4  the  normal  force  derivative  can  be  obtained  from  an 
analysis  of  the  swerving  motion.  However,  for  cases  where  the  model  /  projectile  does  not 
exhibit  any  significant  swerving  motion,  i.e.  a  heavy  model  and/or  one  which  has  very 
small  angular  motions,  the  accuracy  of  the  determined  values  may  be  suspect.  For 
these  cases  another  method  of  determining  Cnq  exists  but  which  requires  the  moment 
data  to  be  available  for  two  different  reference  locations. 

Rewriting  Eq.  7.29  for  two  arbitrary  reference  locations,  A  and  B,  one  can  obtain 
Eq.  7.30  as  shown  below: 

CN„(L/d)=  (C,„A-C„aB)/[(XA/L)-(XB/L)]  7.30 

This  equation  can  be  illustrated  graphically  as  shown  in  the  sketch  below. 


This  sketch  shows  that  the  normal  force  derivative  can  be  determined  from  the  slope  of 
Cmot  vs  x/L.  In  application  this  method  consists  of  launching  models  with  different  eg 
locations  at  the  same  Mach  munbers,  determining  their  respective  pitching  moment 
derivatives,  and  computing  C^a  from  the  indicated  slope  using  Eq.  7.30.  Since  the 
pitching  moment  derivative  can  be  determined  very  accurately  even  from  relatively  small 
angular  motions  this  method  can  provide  good  estimates  for  the  normal  force  when  the 
results  obtained  by  fitting  the  swerving  motion  are  questionable. 
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7.4  DISPERSION 


Frequently  there  is  confusion  over  the  definitions  of  dispersion  and  accuracy  and  as  a 
result  these  terms  can  be  misused.  Simply  put  dispersion  is  the  measure  of  the 
distribution  of  impacts  whereas  accuracy  is  a  measure  of  the  difference  between  the 
center  of  the  dispersion  distribution  and  the  aim  point.  These  definitions  are  illustrated  in 
the  sketches  shown  in  Fig.  7.2. 
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O  0 
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+  O 
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o 

+  ■^Aim  Point  ( typ ) 

-  Impact  Point  ( typ ) 

a.  High  Dispersion,  Inaccurate  b.  High  Dispersion, Accurate 


c.  Low  Dispersion,  Inaccurate  d.  lqw  Dispersion,  accurate 

Figure  7.2  Illustration  of  Dispersion  and  Accuracy 

Dispersion  can  be  a  function  of  a  multitude  of  parameters  associated  with  the  launcher, 
projectile,  and  platform.  Hence  optimizing  the  dispersion  can  be  a  very  complex  task. 
However,  accuracy  can  normally  be  achieved  by  adjusting  the  sighting  system.  For 
example  a  rifle  possessing  low  dispersion  can  be  ‘sighted  in’  with  a  minimum  number  of 
shots,  sometimes  only  one,  by  simply  adjusting  the  sights  to  the  impact  points. 
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Since  dispersion  can  be  a  critical  performance  parameter  for  any  weapon  system  or 
munition,  i.e.  number  of  projectiles  on  target  or  probability  of  hit  for  a  given  munition,  it 
is  normally  specified  as  an  angular  measurement  (milliradians)  and  experimentally 
measured  to  ensure  it  is  within  specifications.  Caution  must  be  used  in  accomplishing 
these  measurements  because  the  apparent  dispersion  can  be  a  function  of  the  distance 
over  which  the  measurements  are  obtained.  This  is  especially  true  if  the  dispersion 
measurements  are  obtained  over  a  relatively  short  flight  distance  such  as  a  1000  calibers 
or  less,  as  will  be  shown  in  the  following  section.  Generally  the  longer  distance  over 
which  the  measurements  are  taken  the  better,  i.e.  several  thousand  calibers  should  be  the 
goal. 

7.4.1  Predictions 

It  is  difficult  at  best  to  compute  the  dispersion  of  a  complete  weapon  system. 

However,  the  dispersion  associated  with  the  aerodynamic  forces  acting  on  a  ffee-flight 
object  can  be  reasonably  estimated.  If  we  start  with  Eq.  6.36  and  neglect  trim  the 
swerving  motion  is  only  a  function  of  the  lift  force  derivative  and  the  angular  motion,  or: 

( y  +  iz )"  =  (  pA  /  2  mg)  [  CtaC  P  +  ia )] 

This  relation  is  an  oscillatory  function  with  range;  however,  if  the  body  is  dynamically 
stable(i.e.  the  motion  damps)  this  relation  tends  to  degenerate  into  a  straight  line  at 
extended  ranges.  Since  the  first  maximum  angle-of-attack  will  also  be  the  largest  angle-of 
attack  experienced  during  the  flight  this  straight  line  will  not  be  parallel  to  the  bore  sight 
line.  This  deviation  from  the  bore  sight  line  is  called  ‘aerodynamic  jump’,  (A.J.),  and  is 
defined  mathematically  as  shown  below: 

X  X 

A.J.  =  (y  +  iz)  /  X  =  lim  (l/x)(pA/2ma)  Clk  J  J  (P  +  ia)dxdx  7.3 1 

x->oo  0  0 

From  Eq.  6. 14,  assuming  Aj  =  0,  we  can  write 

(P  +  ia)  =  [(P  +  ia)"  +  (  H  -  iP  )  (p  +  ia)']  /  (  M  +  iPT ) 

Now  if  this  equation  is  substituted  into  Eq.  7.3 1  and  integrated  by  parts  twice  we  can 
obtain, 

A.J.  =  lim(pA/2ma)[  CLa/(  M  +  iPT  )]{(l/x)(P  +  ia)'  -  (x/x)(p  +  ia)o' 

X-»oo 

+  [(l/x)(p  +  ia)  -  (x/x)(P  +  ia)o]  ( H  -  iP )} 

Taking  the  indicated  limit  we  obtain, 

A.J.  =  [(pA/2m3)  Cl«/(  M  +  iPT )]{-  (P  +  ia)o'  -  (P  +  ia)o  ( H  -  iP  )} 


where 

(P  +  ia)o'  »  qo  ,  initial  angular  cross  rate 
(P  +  ia)o  »  6o  ,  initial  angular  displacement 
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If  we  neglect  the  density  squared  terms,  i.e.  «  0,  and  assume  that  PT  is  small  compared 
to  M  ( M  is  in  the  Cl  plane  where  PT  is  the  side  moment),  then 

A.J.  =  (CLa/  M){-  (pA/2mJ  qo  +  PSq  } 

Substituting  Eq.  6.16  for  M,  Eq.  6.18  for  P,  and  Eq.  7.12  for  kt'^  we  finally  arrive  at  the 
relation  shown  below: 

A. J.  =  (d  CJ  C,„){-  (ly  qo  /  d'mj  +  i(I,  p5o/dVj  }  7.32 

Equation  7.32  can  be  used  to  estimate  the  dispersion  caused  by  the  aerodynamic  lift 
force  derivative  at  extended  ranges.  Note  that  this  relationship  does  not  include  the 
oscillatory  portion  of  the  swerve  caused  by  the  transient  angular  motion.  Also  recognize 
that  this  equation  ignores  the  dispersion  due  to  trim.  For  bodies  with  near  zero  spin 
containing  aerodynamic  or  mass  asymmetries  the  resulting  aerodynamic  trim  can  cause 
large  dispersions  since  this  trim  force  is  always  oriented  in  the  same  direction.  However, 
this  trim  dispersion  can  be  negated  by  purposely  causing  the  model  to  spin  thereby 
changing  the  orientation  of  the  trim  vector. 

If  the  spin  rate  becomes  very  large,  i.e.  spin  stabilized  projectiles,  Eq.  7.32  can  be 
used  to  compute  the  aerodynamic  jump  caused  by  the  in-bore  clearances  of  the  launch 
tube.  For  this  condition  the  initial  cross  rate(qo)can  be  assumed  to  be  equal  to  ip6o  where 
now  5o  is  the  angle  in  the  launch  tube  due  to  gun  clearances.  Therefore  Eq.  7.32  can  be 
written  as, 

A. J.  =  -  i  Clcc  6o  p(Iy  -  Ix)  /  (C„„  d  mj  7.33 

The  angle  due  to  gun  clearances,  6o,  can  be  easily  computed  as  shown  in  the  sketch 
below. 


d 


^ 

■-WB 

Or,  6o  =  (db-d)/LwB  7.34 

where:  d^  =  diameter  of  the  barrel 

d  =  diameter  of  the  projectile 
L^b  =  length  of  the  projectile’s  wheel  base 
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The  wheel  base  length  can  sometimes  be  difficult  to  estimate  since  most  projectiles 
do  not  have  easily  defined  points  of  contact.  For  example  the  model’s  base  may  be  boat- 
tailed  with  a  rotating  band  present  and  the  nose  section  may  consist  of  an  ogive  shape. 
Nevertheless  L^b  can  be  estimated  and  Eq.  7.33  in  conjunction  with  7.34  has  been 
successfully  used  to  compute  the  dispersion  of  numerous  configurations  of  various  sizes 
and  shapes  of  ammvmition.  The  reader  should  refer  to  Ref.  6. 1  for  a  much  better  treatment 
of  aerodynamic  jump. 

7.4.2  Measurements 

When  experimentally  measuring  dispersion  several  factors  should  be  remembered 
and  considered.  Some  of  these  are  outlined  below. 

a.  The  measurements  should  be  obtained  at  a  far  enough  distance  such  that  the  transient 
angle-of-attack  effects  have  dampened  out.  As  mentioned  previously  this  should 
normally  be  several  thousand  calibers. 

b.  The  measurements  are  obtained  with  respect  to  the  center  of  the  impact  pattern,  not 
the  aim  point. 

c.  The  laimcher  and  mounting  system  should  not  be  altered  in  any  way  during  the  tests. 
If  it  is  absolutely  necessary  to  move  the  gun  or  mounts  the  resulting  dispersion  data 
must  be  corrected  in  some  manner,  i.e.  the  data  could  be  corrected  using  the  before 
and  after  aim  point.  However,  it  is  preferable  to  start  a  new  statistical  base  based  on 
the  new  impact  pattern. 

d.  If  the  muzzle  velocities  vary  significantly  during  the  tests,  the  data  must  be  corrected 
for  gravity  drop. 

e.  Each  projectile  configuration  will  have  its  own  dispersion  pattern.  Therefore  enough 
rounds  of  each  configuration  should  be  fired  to  obtain  a  statistical  base  for  that 
configuration. 

Although  the  individual  measurements  are  normally  presented  in  mils,  the  dispersion 
pattern  is  frequently  described  as  circular  error  probable,  CEP.  The  CEP  of  the  dispersion 
pattern  is  defined  as  the  radius  of  the  circle  in  which  50  percent  of  the  rounds  will  impact. 
Usually,  with  only  minimal  effort,  dispersion  measurements  can  be  obtained  in  the  ffee- 
flight  spark  range  at  the  same  time  the  aerodynamic  tests  are  ongoing.  Frequently  this  is 
not  accomplished  either  because  it  is  believed  that  the  data  would  not  be  representative  of 
the  actual  configuration  or  because  something  was  done  during  testing  in  order  to  excite 
the  angular  motion.  Nevertheless  the  test  engineer  should  note  that  this  data  can  be 
obtained  at  very  small  additional  expense  and  effort  if  it  is  required  and  appropriate  for 
his  particular  tests. 


Ill 
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VIII.  FUTURE  REQUIREMENTS 


It  is  always  dangerous  in  attempting  to  predict  what  the  future  might  bring  to  any  line 
of  endeavor.  However,  it  is  somewhat  less  hazardous  here  since  each  of  the  technologies 
discussed  have  efforts  already  on-going  and  some  have  been  in  consideration  for  several 
years.  Other  than  the  interesting  aerodynamic  issues  a  researcher  can  be  involved  with 
there  are  also  three  primary  support  technologies  which  are  always  of  interest  to  the 
ballistician.  These  three  areas  are  advanced  instrumentation  systems,  improved  launchers 
and  launch  techniques,  and  new  theoretical  and  trajectory  analysis  algorithms.  A  brief 
discussion  of  the  on-going  efforts  in  these  three  areas  will  be  presented  in  this  section. 

8.1  INSTRUMENTATION 

8.1.1  Flow  Field  Visualization  and  Measurements 

Over  the  years  the  study  of  basic  fluid  flow  physics  has  required  the  development  of 
flow  visualization  techniques.  These  have  included  smoke  generation,  shadowgraph  and 
schlieren  photography  along  with  relatively  complex  interferometry  techniques.  The 
aerodynamist  is  also  constantly  attempting  to  relate  changes  /  variations  in  the  measured 
aerodynamic  coefficients  and  derivatives  with  observable  phenomenon  occurring  in  the 
flow  about  the  body.  For  these  reasons  and  others,  flow  visualization  is  of  interest  to  the 
researcher.  Recently  with  the  development  of  computational  fluid  dynamics  (  CFD  ) 
codes  a  requirement  to  be  able  to  validate  these  codes  with  experimental  data  has  arisen. 
Although  surface  pressure  and  non-obtrusive  velocity  data  can  be  obtained  in  the  wind 
tunnel,  non-obtrusive  pressure  and  density  measurements  in  the  far  flow  fields  or  wake 
can  be  difficult  if  not  impossible  to  obtain. 

In  order  to  obtain  basic  non-obtrusive  measurements  an  effort  was  initiated  in  the 
early  1990’s  to  develop  a  multiview  holographic  interferometry  system,  and  the 
corresponding  data  analysis  algorithms  for  mapping  the  flow  field  density  *  These 
flow  field  density  maps  would  then  be  used  to  assist  in  validating  the  corresponding  CFD 
predictions  or  to  develop  advanced  CFD  turbulence  models,  etc.  A  schematic  of  the 
holocamera  developed  for  this  purpose  is  presented  in  Fig.  8.1  and  was  positioned  in 
the  up-range  end  of  the  ARF.  Several  basic  model  configurations  were  flown  through  this 
holocamera  and  6  different  views  obtained  (  see  interferograms  in  Fig.  8.2  ).  These 
images  were  then  passed  into  the  phase  shift  interferometry  algorithm  *  ^  for  analysis. 
Although  this  algorithm  was  successfully  demonstrated  using  generated  data  the 
existence  of  experimental  noise  ( primarily  Laser  speckle  )  in  the  images  caused  the  phase 
map  unwrapping  process  to  fail. 

It  is  believed  that  the  primary  source  of  the  Laser  speckle  comes  from  the  scatter 
plates  in  the  holocamera,  see  Fig.  8. 1 .  In  fact,  preliminary  bench  tests  *  ^  using  a  mirror  in 
place  of  a  scatter  plate  tended  to  confirm  this  belief  As  a  result,  an  effort  to  continue 
development  of  this  technique  is  currently  in  progress.  The  initial  approach  will  be  to 
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Figure  8.1  Three  Dimensional  Layout  of  the  Holocamera 
(  Fig.  2  from  Ref.  8.1  ) 

replace  one  of  the  scatter  plates  with  a  first  surface  mirror  and  obtain  a  holographic 
interferogram  of  a  free-flight  object.  This  interferogram  will  be  duplicated  5  times  and  the 
corresponding  6  images  passed  into  the  phase  shift  interferometry  algorithm  for  analysis 
and  evaluation.  Assuming  that  a  successful  result  will  be  achieved  the  remaining  scatter 
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plates  in  the  holocamera  will  be  replaced  with  similar  mirrors.  This  would  then  complete 
a  system  which  would  be  capable  of  determining  the  flow  field  densities  about  ffee-flight 
missile  configurations  and  provide  heretofore  unobtainable  data  for  analysis  and 
comparison  with  theoretical  predictions. 


Figure  8.2  Six  Simultaneous  Interferograms  of  a  Cone  Cylinder  Flare 
Configuration  at  Mach  3.  ( Fig.  13  from  Ref.  8.1) 

8.1.2  Electronic  Imaging 

As  was  discussed  in  Section  5.1,  also  see  Ref.  5.4,  it  has  been  planned  for  some  time 
to  convert  the  existing  film  based  system  in  the  ARF  to  an  electronic  imaging  system 
based  on  modem  charge  coupled  device  (CCD)  technology.  The  first  part  of  this 
conversion  was  the  procurement  of  the  Ektron  film  scanner  and  the  development  of  the 
associated  image  processing  system  (CADRA)  .  Now  that  the  image  processing  system 
is  operational  and  used  on  a  daily  basis  the  primary  task  still  requiring  completion  is  the 
procurement  and  installation  of  the  CCD  cameras.  The  cameras  chosen  for  this 
installation  are  EG&G,  Optoelectronics-Reticon,  high  resolution  devices.  It  is  planned  to 
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install  Ik  X  Ik  pixel  array  cameras  in  the  up-range  end  of  the  facility  and  2k  x  2k  devices 
in  the  down-range  section  of  the  ARF.  It  is  believed  that  the  higher  resolution  is  required 
down-range  for  two  reasons.  First  the  cross  sectional  dimensions  of  the  facility  are  larger 
down-range  and  secondly,  because  of  potential  projectile  dispersion,  the  area  to  be 
viewed  is  significantly  larger.  The  network  connecting  the  cameras  to  the  control  room, 
data  collection  and  analysis  system  has  been  designed  and  will  incorporate  an  ethemet 
interface  and  protocol.  A  schematic  of  the  electronic  shadowgraph  system  including  the 
existing  film  scanner  is  included  in  Fig.  8.3. 


FUUWflEAOER-SYSTBM  TAPE-BACKUP 


Figure  8.3  Schematic  of  Electronic  Shadowgraph  System 
(  Fig.  5  fi'om  Ref.  5.4  ) 
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Once  this  system  has  been  installed  the  complete  process  from  test  firing  to  having 
the  trajectories  ready  for  analysis  will  be  minutes  instead  of  days  as  is  presently  required. 
The  film  scanning  and  digitization  is  the  single  most  manpower  intensive  operation  at 
the  present  time  and  this  along  with  film  handling  and  processing  will  be  completely 
eliminated  with  the  installation  of  the  electronic  cameras.  However,  because  of  limited 
resources,  it  is  contemplated  that  completion  may  require  several  years  to  procure  the  1 00 
cameras  required  to  outfit  the  facility. 

8.1.3  On-Board  Instrumentation 

On-board  instrumentation  techniques  have  been  used  in  aerodynamic  testing  for  many 
years  especially  for  the  development  of  manned  vehicles  or  relatively  large  unmanned 
systems.  Objects  as  small  as  artillery  rounds  or  sounding  rockets  *  have  also  utilized 
on-board  instrumentation;  however,  such  packages  have  not  been  used  extensively  for 
test  items  in  the  20-30  mm  diameter  range.  This  is  primarily  due  to  the  relatively  small 
volumes  available  for  packaging  the  desired  instrumentation  and  the  problems  associated 
with  recovering  the  data.  Historically,  two  different  approaches  to  recovering  the  data 
have  been  used.  The  first  is  to  transmit  the  data  to  a  receiving  station  while  the  test  item  is 
in  flight.  This  approach  requires  the  model  to  not  only  carry  the  desired  sensor(s)  but  a 
miniaturized  transmitter  and  antenna  as  well.  The  second  approach  is  to  record  the  data 
on-board  and  to  recover  the  test  item  and  remove  the  stored  data.  This  approach  requires 
a  miniaturized  recording  device  that  must  survive  the  flight  termination  event  and  also  be 
recoverable.  Although  both  techniques  have  been  successfully  used  in  the  past  neither 
has  been  particularly  beneficial  to  free-flight  testing  of  subscale  items. 

This  may  change  in  the  very  near  future  with  the  development  of  a  family  of 
hardened  subminiature  sensor  and  telemetry  packages.  Such  a  family  of  on-board 
instrumentation  is  presently  being  developed  by  the  U.S  Army  *  It  is  hoped  that  this 
technology  demonstration  program  will  develop  instrumentation  for  the  continuous 
measurement  of  projectile  attitude,  pressure,  temperature,  accelerations  and  other  flight 
data  at  reasonable  costs.  If  this  program  is  successful  it  could  revolutionize  free-flight 
aerodynamic  testing  as  we  know  it  and  open  numerous  horizons  heretofore  unimaginable. 

8.2  LAUNCH  TECHNIQUES 

8.2.1  Launchers 

Ballistic  researchers  /  practitioners  are  always  interested  in  improving  the  existing 
launchers  and  developing  new  launch  techniques.  This  interest  stems  from  the  desire  to 
launch  larger  and  heavier  packages  at  higher  and  higher  velocities  or  to  subject  existing 
packages  to  reduced  in-bore  acceleration  loads.  There  is  always  the  need  to  expand  the 
existing  operational  envelopes  into  regions  unobtainable  in  the  past.  Also,  if  the  in-bore 
accelerations  can  be  reduced  while  maintaining  the  muzzle  velocity,  relatively  fragile 
model  /  sabot  packages  can  be  launched  and  still  survive  the  launch  loads.  In  general 
these  are  the  considerations  which  led  to  the  development  of  two-stage  light  gas  gims 
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for  use  in  various  laboratories  and  initially  spurred  the  interest  in  other  more  advanced 
concepts,  i.e.,  electromagnetic  launchers,  ram  accelerators,  etc. 

During  a  preliminary  study  into  the  possibility  of  weaponizing  a  two-stage  light 
gas  gun  a  rather  unique  two-stage  firing  cycle  was  discovered  *  In  contrast  to  the 
conventional  two-stage  light  gas  gun  this  new  launcher(  wave  gun )  employed  a  very  light 
piston  in  a  short  pump  tube  (low  volume)  containing  a  light  gas,  i.e.  hydrogen  or  helium, 
at  relatively  high  pressures.  These  conditions  permitted  the  piston  to  oscillate  in  the  pump 
tube  prior  to  the  rupture  of  the  diaphragm  and  the  resulting  acceleration  of  the  model  / 
sabot  package  down  the  launch  tube.  This  firing  cycle  potentially  can  minimize  the  in¬ 
bore  accelerations  for  a  given  muzzle  velocity.  Also  portions  of  the  existing  two-stage 
light  gas  gun  hardware  could  be  used  with  only  minor  changes.  Therefore  this  appears  to 
be  an  attractive  technique  for  use  in  free-flight  experimental  /  research  facilities.  A 
comparison  of  the  two  types  of  cycles  are  shown  in  Fig.  8.4. 


Conventional  Li^t  Gas  Gun 


Piston 


Wave  Gun 
Piston 


Projectile 


Projectile 


Heavy  pistcxL 

Low  charge  piessure. 

Large  pump  tube  volume. 


□0= 

Light  piston. 

Hi^  charge  pressure. 
Small  puny  tube  vohnne. 


Figure  8.4  Firing  Cycles  for  a  Conventional  Two-Stage  Light  Gas  Gim  and  a  Wave  Gun 
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A  recent  investigation  was  conducted  to  develop  a  mathematical  model  of  a  wave  gun 
and  validate  this  model  with  some  experimental  tests  *  That  investigation  was 
successful  leaving  two  questions  still  to  be  determined.  First,  does  the  wave  gun  firing 
cycle  illustrated  in  Fig.  8.4  possess  enough  advantages  over  the  conventional  two-stage 
light  gas  gun  to  justify  the  required  modifications.  Secondly,  what  are  the  optimum 
dimensions  of  the  necessary  hardware  modifications,  i.e.,  pump  tube  and  piston.  An 
effort  is  on-going  to  answer  these  questions  and  depending  on  the  results  the  existing 
two-stage  light  gas  gun  will  either  be  modified  or  maintained  in  its  present 
configuration. 

8.2.2  Sabot  Design 

As  mentioned  in  Section  3.3,  at  the  present  time  the  designer  does  not  have  access  to 
a  user  friendly  sabot  design  software  package.  There  are  several  existing  finite  element 
software  packages  which  have  been  successfully  used  for  point  designs,  however,  they 
could  not  be  described  as  user  friendly.  Typically  they  require  highly  trained  and 
experienced  analysts  and  have  been  developed  to  handle  a  multitude  of  diverse  and 
generalized  design  problems.  What  is  required  for  the  ballistics  researcher  is  an  easily 
used  code  and  design  process  developed  solely  for  sabots.  This  requirement  was 
recognized  and  an  on-going  effort  established  in  an  attempt  to  provide  the  designer 
with  such  tools,  see  Ref  8.9.  This  reference  describes  the  on-going  program  and  also 
contains  an  excellent  bibliography  of  sabot  technology.  Essentially  the  final  output 
of  the  effort  will  be  a  design  code  which  will  contain  a  finite  element  routine  specifically 
written  for  sabots  and  simple  closed  form  solution  routines  based  on  past  point  designs. 
The  approach  will  be  to  develop  several  basic  sabot  designs  using  these  codes  which  are 
expected  to  fail  under  certain  launch  conditions.  These  designs  will  then  be  built  and 
tested  at  conditions  in  and  aroimd  the  expected  failure  conditions.  The  results  of  these 
tests  will  provide  the  required  data  to  evaluate  the  design  codes. 

In  conjunction  with  the  effort  described  above  a  simple  mechanism  to  evaluate  the 
mechanical  properties  of  various  materials  commonly  used  in  sabot  construction  will 
also  be  constructed  and  evaluated,  see  sketch  below. 
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It  is  suspected  that  the  material  properties  of  some  of  these  materials,  i.e.,  polycarbonates, 
may  vary  significantly  from  batch  to  batch  or  with  different  manufacturers.  There  is  also 
considerable  concern  that  the  material  properties  may  change  over  time.  Hence  sabots 
constructed  in  the  past  and  that  were  successful  may  fail  at  identical  test  conditions  in  the 
future.  That  certainly  is  a  designer’s  nightmare  and  this  simple  testing  apparatus  may  be 
sufficient  to  guard  against  that  possibility. 

Undoubtedly  there  will  be  several  publications,  papers  and  reports,  documenting  the 
results  of  this  on-going  effort  in  the  relatively  near  future.  It  is  hoped  that  this  work  and 
the  associated  publications  will  take  some  of  the  art  out  of  sabot  design  and  put  it  on  a 
firm  scientific  basis. 

8.3  TRAJECTORY  ANALYSIS 

8.3.1  Aerodynamic  Modeling 

Considerable  progress  has  been  made  in  modeling  the  aerodynamic  coefficients  and 
derivatives  as  functions  of  angle  of  attack,  Mach  number,  and  roll  orientation  during  the 
past,  see  the  discussions  on  the  quasi-nonlinear  analysis  and  numerical  integration  in 
Sections  6.2  and  6.3  and  the  references  at  the  end  of  Section  VI.  All  of  these  existing 
techniques  assume  that  aerodyneimic  parameters  can  be  expanded  as  predetermined 
mathematical  functions,  normally  polynomials  (see  the  expansions  presented  in  Section 
2.3).  The  analysis  routines  then  determine  the  coefficients  of  these  assumed  functions 
such  that  the  theoretical  equations  of  motion  best  match  the  experimentally  measured 
motions.  For  a  great  majority  of  the  test  flights  this  process  works  exceptionally  well  and 
the  resulting  aerodynamics  adequately  describes  the  experimental  motion  profiles. 

However,  there  are  configurations  of  interest  where  these  predetermined 
mathematical  functions  do  not  adequately  describe  the  aerodynamics.  For  example  any 
configuration  which  possesses  small  angle-of-attack  ( -3  to  +  3  degs  )  nonlinearities 
cannot  be  adequately  modeled  with  polynomials  since  these  functions  inherently  describe 
nonlinearities  which  increase  with  increasing  angle-of-attack.  Therefore  any  polynomial 
which  would  adequately  describe  a  small  angle-of-attack  nonlinearity  would  greatly 
exaggerate  the  nonlinearities  at  the  higher  angles.  Also,  these  functions  can  pose  serious 
restrictions  for  almost  any  configuration  at  very  high  angles-of-attack  ( angles  greater  that 
30  degs).  Finally  this  existing  process  imposes  a  heavy  burden  on  the  analyst  in  selecting 
the  correct  form  of  the  predetermined  mathematical  functions  and  requires  considerable 
experience  in  recognizing  which  of  the  aerodynamic  expansions  are  likely  required 
to  successfully  describe  the  experimentally  measured  motions. 

Considering  the  above  comments  the  ffee-flight  trajectory  analyst  is  always  vigilant 
in  looking  for  advanced  mathematical  modeling  techniques  or  processes  which  could 
alleviate  the  necessity  of  assuming  a  known  form  for  the  various  nonlinearities.  Recently 
three  such  techniques  were  investigated  and  discussed  in  Ref.  5.2.  These  three  techniques 
were  identified  as  splines,  the  Mook  technique,  and  the  Mook-Taylor  series  technique. 
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The  splines  investigated  in  Ref.  5.2  utilized  continuous  piecemeal  linear  functions. 
However,  this  technique  would  not  be  restricted  to  linear  functions  and  higher  order 
polynomials  could  be  splined  together  in  a  similar  manner.  Essentially  this  process  would 
describe  the  nonlinear  aerodynamics  as  a  series  of  short  line  segments  where  the  analyst 
could  select  the  interval  for  the  knot  locations,  i.e.,  every  degree  or  so.  In  the  Mook 
technique  no  mathematical  model  for  the  aerodynamic  parameters  are  specified  and  the 
aerodynamic  coefficients  and  derivatives  are  computed  at  each  experimentally  measured 
data  point.  Therefore  this  technique  provides  the  aerodynamics  as  a  function  of  time  not 
angle-of-attack,  Mach  number,  or  roll  orientation.  In  the  Mook-Taylor  series  technique 
the  Mook  terms  represent  the  Taylor  series  corrections. 

All  of  the  above  data  correlation  techniques  have  been  coded  into  the  CADRA 
software  package  along  with  the  conventional  polynomial  expeinsions.  But,  other  than 
what  was  done  in  Ref  5.2,  a  detailed  and  systematic  study  of  the  various  options  have  not 
been  accomplished  to  date.  Presently  what  is  required  is  a  pressing  reason  and  need  to 
accomplish  such  a  detailed  investigation.  A  set  of  experimental  data  which  cannot  be 
adequately  analyzed  using  the  conventional  polynomial  expansions  would  provide  this 
reason  and  need. 

If  one  of  these  or  any  other  advanced  data  correlation  technique  proves  successful  in 
the  future  it  is  believed  that  it  would  be  incorporated  into  the  existing  ARFDAS  routine, 
see  Fig.  6.1,  instead  of  replacing  the  existing  system.  The  idea  being  that  the  trajectory 
analysis  would  proceed  exactly  the  same  as  now  using  the  polynomial  expansions  for  the 
aerodynamics.  However,  the  last  step  would  then  be  to  pass  these  results  into  the  new 
technique  as  the  initial  starting  point  for  the  advanced  data  correlation  routine  and 
accomplish  one  additional  level  of  analysis.  This  last  step  may  prove  not  to  be  necessary 
for  all  configurations  but  would  be  available  for  those  that  required  a  more  general 
aerodynamic  model. 

8.3.2  AUTOMATION 

At  the  present  time  free-flight  data  collection  and  trajectory  analysis  requires  a 
considerable  amoimt  of  hands-on  interaction  of  the  facility  operators  and  analysts  and  are 
manpower  intensive.  This  hands-on  interaction  starts  with  the  film  handling,  i.e.,  loading 
film  into  cassettes,  removing  exposed  film  from  the  cassettes,  film  processing,  and  film 
scanning  and  digitization  into  the  CADRA  software  system.  The  final  completion  and 
installation  of  the  electronic  shadowgraph  system  as  described  in  Section  8.1.2  will 
eliminate  these  manpower  intensive  operations.  In  fact  once  the  electronic  cameras  are 
installed  the  CCD  images  will  be  available  for  trajectory  computations  immediately  after 
the  test  flight  is  terminated.  The  CADRA  system  has  an  automatic  image  scanning 
routine  already  imbedded  in  the  software.  This  automatic  image  scanning  routine  has 
been  checked  out  using  ideally  generated  images  but  has  not  been  utilized  for 
experimentally  obtained  images.  This  automatic  routine  requires  that  the  reference  points 
in  the  images  be  somewhat  consistently  positioned  each  time  an  image  at  a  given 
shadowgraph  station  is  obtained.  At  the  present  time  the  film  cannot  be  manually 
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positioned  on  the  scanner  with  enough  precision  and  consistency  for  the  automatic 
routine  to  locate  the  reference  points.  However,  when  the  CCD  images  arrive  for 
scanning  they  will  be  consistently  positioned  and  the  automatic  routine  can  be  utilized. 
This  not  only  means  that  the  images  are  available  for  immediate  analysis  but  that  within 
minutes  after  the  test  flight  the  experimental  trajectory  will  be  available  for  analysis  and 
extraction  of  the  aerodynamic  coefficients  and  derivatives. 

In  the  past  few  years  considerable  effort  has  also  been  expended  in  automating  much 
of  the  trajectory  analysis  process,  see  for  example  Ref  8.10.  This  reference  describes  an 
expert  system  based  on  artificial  intelligence  (AI)  techniques.  It  is  believed  that  a  similar 
system  in  conjunction  with  the  automatic  image  processing  and  trajectory  computation 
routines  would  permit  real  time  trajectory  analysis  and  the  determination  of  the 
aerodynamic  coefficients  and  derivatives  for  many  if  not  most  of  the  test  configurations. 
Certainly  this  could  be  accomplished  for  most  configurations  within  a  hour  or  two  of  the 
test  flight  and  allow  the  user  and  test  engineer  to  consider  these  results  prior  to  launching 
the  next  test  model. 

Continuing  this  process  of  automating  the  data  collection  and  trajectory  analysis 
process  is  an  absolute  necessity  in  the  future.  The  benefits  which  will  be  derived  from 
this  process  include: 

(1)  the  elimination  of  potential  human  errors  which  presently  exist  in  the  existing 
manpower  intensive  operations 

(2)  reduced  man-hours  required,  i.e.,  film  handling,  processing,  scanning 

(3)  reduced  material  resources  required,  i.e.,  film,  processing  chemicals 

(4)  immediate  availability  of  aerodynamic  data  for  programmatic  decisions 

(5)  improved  facility  efficiency  ( more  tests  in  less  time) 

(6)  improved  user  perceptions  associated  with  perceived  outdated  processes  and 
techniques 

Automation  is  the  key  to  the  future  viability  of  free-flight  spark  ranges.  Such  facilities 
have  and  should  continue  to  provide  a  unique  role  in  aerodynamic  testing  but  like  any 
other  tool  they  must  be  maintained  and  improved  to  meet  future  requirements. 
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APPENDIX  1 

LOCATION  OF  INSTRUMENTATION  IN  THE 
AEROBALLISTIC  RESEARCH  FACILITY 
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WINDOW 

NUMBER 


ARF  INSTRUMENTATION  LOCATIONS 

xSnce  instrumentation  INSTRUMENTATION 


INTERFEROMETER  1 


S.G.  STA.  1 


10.0  *  S.G.  STA.  2.  T.C.  1 


S.G,  STA.  3 


20.0  I  INTERFEROMETER  2 


S.G.  STA.  4 


DIRECTS.  G.  STA. 


A.C.  CONTROL 


S.G.  STA.  5 


S.G.  STA.  6 


, _ i^ 

mam 


S.G.  STA.  7 


S.G.  STA.  8 


S.G.  STA.  9 


MULTI.  SPARK  1 


T.C.  2 


S.G.  STA.  10 


WINDOW  NOMINAL 
NUMBER  X-DISTANCE 


35  170.0  S.G.  STA.  16,  T.C.  3 


S.G.  STA.  17 


S.G.  STA.  18 


MULTI.  SPARKS 


S.G.  STA.  19 


ELATIVE  HUMIDIT 


S.G.  STA.  20 


S.G.  STA.  21 


S.G.  STA.  22 


S.G.  STA.  23  ,  A.C.  CON. 


36 

175.0 

37 

180.0 

38 

185.0 

39 

190.0 

40 

195.0 

41 

209.8 

42 

214.8 

219.8 

44 

224.8 

45 

229.8 

46 

234.8 

49 


50  254.8 


259.8  S.G.  STA.  24  ,  T.C.  4 


52  264.8 


53 

269.8 

S.G.  STA.  25 

54 

274.8 

S.G.  STA.  26 

110.0 


115.0 


125.0 


130.0 


135.0 


140.0 


145.0 


150.0 


160.0 


S.G.  STA.  11 


S.G.  STA.  12 


S.G.  STA.  13 


S.G.  STA.  14 


S.G.  STA.  15 


MULTI.  SPARK  2 


57 

289.8 

58 

294.8 

mumi 

61 

309.8 

63 

319.9 

mm 

65 

329.9 

S.G.  STA.  27 


S.G.  STA.  28 


S.G.  STA.  29 


S.G.  STA.  30 


WINDOW 

NUMBER 


ARF  INSTRUMENTATION  LOCATIONS 

x°dSce  instrumentation  INSTRUMENTATION 


WINDOW  NOMINAL 
NUMBER  X-DISTANCE 


339.9  T.c.  5 _ 

344.9  _ 

349.9  S.G.  STA.  31  ,  A.C.  CON. 

354.9  _ 

359.9  _ 

365.0  multi,  spark  4 

370.0 _ 

375.0 _ 

380.0  S.G.  STA.  32 

385.0 _ 

390.0 _ 

395.0  SQ-  STA.  33 

400.0 _ 

405.0 _ 

410.0  S.G.  STA.  34 

415.0 _ 

420.0 _ 

425.0  ® 

430.0 _ 

435.1  _ 

440.1  S.G.  STA.  36 

445.1  _ 

450.1  A.C.  CONTROL 

455.1  S.G.STA.37 

460.1  _ 

465.1  _ 

470.1  S.G.  STA.  38 

475.1  _ 

480.1  _ 

435.1  S.G.  STA.  39 

490.1  _ 

495.1  _ 

500.1  S.G.  STA.  40,  T.C.  6 


100 

505.1 

101 

510.1 

102 

I 

515.1 

103 

520.1 

104 

525.1 

105 

530.1 

106 

535.1 

107 

540.1 

108 

545.1 

109 

550.1 

110 

555.1 

111 

560.1 

112 

565.1 

113 

570.1 

114 

575.1 

115 

580.2 

116 

585.2 

117 

590.2 

118 

595.2 

119 

600.2 

120 

605.2 

121 

610.2 

122 

615.2 

123 

620.2 

124 

625.1 

125 

630.1 

126 

635.1 

127 

640.1 

128 

645.1 

129 

650.1 

130 

655.1 

131 

660.1 

NC 

)TE;  T.C. 

S.G.  STA.  41 

S.G.  STA.  42 

S.G.  STA.  43 
A.C.  CONTROL 

S.G.  STA.  44 

MULTI.  SPARK  5 

T.C.  7 

S.G.  STA.  45 

S.G.  STA.  46 

S.G.  STA.  47 

S.G.  STA.  48 

S.G.  STA.  49 

A.C.  CONTROL 
S.G.  STA.  50 


PH  STATION  # 
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APPENDIX  2 


FACILITY  CALIBRATION  TECHNIQUE 

( Extracted  from  Ref.  3.13) 


129 


A  2.1  BASIC  CONCEPT 


When  developing  the  calibration  technique  several  features  appear  advantageous.  It 
was  hoped  that  most  if  not  all  of  these  features  could  be  incorporated  into  the  final 
technique.  These  features  were: 

a.  no  surveying  should  be  required, 

b.  technique  to  be  accomplished  by  the  normal  operating  crew, 

c.  the  range  coordinates  of  all  reference  beads  and  spark  source  orifices  ( 100  of  each, 
two  at  each  of  the  50  shadowgraph  stations)  should  be  simultaneously  determined, 

d.  technique  should  require  a  minimum  of  range  downtime, 

e.  a  minimum  of  special  hardware  and  fixtures  should  be  required. 

Working  with  these  considerations  in  mind  the  technique  discussed  herein  was 
developed. 

Early  in  the  development  process  it  was  recognized  that  if  shadowgrams  of  points  in 
space  with  known  range  coordinates  could  be  obtained  that  these  would  represent  an  ideal 
source  of  range  calibration.  The  logic  being  that  once  these  shadowgrams  were  obtained 
the  range  coordinates  of  the  shadows  of  the  known  points  could  be  determined  with 
respect  to  the  reference  system  which  would  likely  contain  errors.  If  the  range  coordinates 
of  the  reference  system  beads  were  not  in  error  then  the  vectors  cormecting  the  shadows 
to  the  associated  known  points  in  space  would  intersect  at  the  spark  source  orifice.  This  is 
obvious  since  the  rays  casting  the  shadows  all  originated  at  the  orifice.  Therefore  if  the 
range  coordinates  of  a  reference  bead  were  in  error,  as  likely  they  would  be,  the  vectors 
connecting  the  shadows  to  the  associated  known  points  would  not  intersect.  It  was  also 
recognized  that  the  error  (other  than  film  reading)  in  the  shadow  of  the  known  point  is 
equal  to  the  error  in  the  range  coordinates  of  the  reference  system  bead.  Note  that  the 
shadow  coordinates  are  obtained  with  respect  to  the  reference  bead.  Therefore  if  the  range 
coordinate  of  the  reference  bead  are  varied  until  the  vectors,  defined  by  the  shadows  and 
the  associated  known  points  in  space,  intersect  then  these  new  coordinates  are  the  true 
range  coordinates  of  the  reference  bead  and  the  point  of  intersection  is  the  spark  source 
orifice  location.  Actually,  since  film  reading  errors  are  also  present  the  above  process 
amoimts  to  minimizing  a  circle  of  intersection  of  the  vectors  and  assuming  that  the  center 
of  this  circle  corresponds  to  the  spark  source  orifice. 

A  2.2  DETERMINATION  OF  KNOWN  POINTS 

The  above  discussion  illustrates  that,  using  this  approach,  the  problem  of  range 
calibration  reduces  to  a  problem  of  positioning  points  in  space  whose  range  coordinates 
are  known.  However,  this  also  can  be  a  serious  problem  especially  when  surveying  is 
undesirable  and  it  is  wished  to  calibrate  all  shadowgraph  stations  simultaneously.  The 
method  selected  for  obtaining  points  with  known  coordinates  consisted  of  a  Kevlar  wire 
with  beads  positioned  precisely  1.524  meters  (5  ft)  apart  strung  the  entire  length  of  the 
instrumented  range.  This  wire  was  ai^usted  such  that  one  bead  was  in  the  field  of  view  of 
both  cameras  at  each  of  the  50  shadowgraph  stations.  This  calibration  wire  had  been 
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constructed  while  under  3.63  Kgm  (8  lb)  of  tension  using  an  Invar  tape  calibrated  every 
1 .524  meters.  When  hanging  in  the  range  the  calibration  wire  was  fixed  at  the  up-range 
end  and  the  downrange  end  was  passed  over  a  near  ffictionless  pulley  with  the  same 
weight  attached.  The  calibration  wire  with  the  beads  attached,  the  calibrated  tape,  and  the 
support  system  at  both  ends  are  all  of  the  special  hardware  required. 

With  this  calibration  wire  in  place  shadowgrams  are  obtained  at  all  100  cameras  (two 
at  each  of  the  50  shadowgraph  locations).  The  calibration  wire  is  then  moved  to  another 
location  both  cross  range  and  vertically  and  the  process  of  obtaining  shadowgrams  at 
each  of  the  50  locations  repeated.  This  process  is  accomplished  for  several  different 
positions  of  the  calibration  wire.  Three  sets  of  shadowgrams  are  to  be  used  to  calibrate 
the  reference  bead  and  spark  source  orifice  locations  and  the  remaining  sets  of 
shadowgrams  will  be  used  to  get  a  statistical  base  for  determining  the  measurement 
precision.  The  positions  of  the  shadows  of  the  beads  on  the  calibration  wire  are  then  read 
and  numerically  coded  as  if  they  were  test  items  in  flight.  Using  the  resulting  measured 
range  coordinates  of  the  calibration  beads  the  actual  range  coordinates  of  the  calibration 
beads  are  determined  statistically  by  using  the  known  physical  state  of  the  calibration 
wire.  This  known  physical  state  consists  of  the  facts  listed  below: 

a.  the  cross-range  coordinates  (y)  plotted  vs  the  down-range  coordinates  (x)  define  a 
straight  line, 

b.  the  vertical-range  coordinates  (z)  plotted  vs  x  define  the  equation  of  a  catenary, 

c.  initially  the  arc  length  between  beads  on  the  calibration  wire  was  1 .524  meters. 

The  mathematical  treatment  as  discussed  in  the  following  sections  permits  the  calibration 
wire  to  stretch  or  contract  assuming  it  occurs  linearly  over  the  entire  length  of  the  range. 

A  2.2.1  CROSS-RANGE  COORDINATES 

The  actual  y  (cross-range)  coordinate  of  each  calibration  bead  is  determined 
statistically  by  least  squares  fitting  a  straight  line  to  the  measured  y  cuid  x  coordinate  data. 
These  measured  y  and  x  coordinates  were  obtained  using  a  previous  range  calibration  or 
initial  estimates  for  the  spark  source  orifice  and  reference  bead  locations. 

Once  this  fit  has  been  obtained,  the  theoretical  y  coordinate  corresponding  to  each 
measured  y  coordinate  represents  the  best  statistical  estimate  of  the  actual  cross-range 
coordinate  for  each  bead  location.  It  should  be  noted  that  when  the  calibration 
shadowgrams  are  obtained  some  care  should  be  exercised  to  ensure  that  the  calibration 
wire  is  nominally  parallel  to  the  range  reference  system.  By  insuring  that  the  wires  are 
nominally  parallel,  large  errors  in  the  measured  x  coordinates  can  be  tolerated  during  the 
fitting  process  and  still  only  insignificant  errors  in  the  theoretically  determined  y 
coordinates  will  result.  Also  ,  when  the  calibration  shadowgrams  are  obtained  the  heating 
/  air  conditioning  system  must  be  turned  off  in  order  minimize  the  effects  of  air  currents 
on  the  calibration  wire.  This  ensures  that  the  wire  hanging  under  its  own  weight  describes 
a  vertical  plane  and  the  straight  line  assumption  is  valid.  Since  the  technique  of  least 
squares  fitting  a  straight  line  to  data  points  is  common  place  it  is  not  elaborated  on  here. 
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A  2.2.2  VERTICAL-RANGE  COORDINATES 

Two  different  statistical  approaches  of  determining  the  actual  z  (vertical)  coordinate 
of  each  calibration  bead  were  investigated.  These  two  approaches  consisted  of  fitting  the 
measured  z  and  x  coordinates  by  either  assuming  that  z  was  a  parabolic  function  of  x,  or 
that  z  was  a  hyperbolic  cosine  function  of  x.  The  actual  equation  of  a  wire  hanging  under 
its  own  weight  ( the  shape  of  which  is  called  a  catenary)  is  a  hyperbolic  cosine  function 
but  the  fitting  process  is  somewhat  more  complex  than  the  parabola.  It  is  also  well  known 
that  a  parabola  is  an  excellent  approximation  to  the  equation  of  a  catenary  when  the  sag 
is  small.  Both  of  these  methods  were  therefore  evaluated  by  fitting  the  measured  z  and  x 
coordinates  of  a  typical  calibration  wire  and  comparing  results.  The  maximum^difference 
between  the  two  theoretical  curves  as  determined  from  these  fits  was  0.34x10  meters 
and  the  average  difference  was  0.14x10"^  meters.  Since  these  differences  are  a  couple 
orders  of  magnitude  smaller  that  our  measuring  capability  it  was  assumed  that  both 
methods  were  essentially  identical  and  the  simpler  parabolic  function  was  used  in  the 
calibration  routines. 

Again,  once  the  fit  has  been  obtained  the  theoretical  z  coordinate  corresponding  to 
each  measured  z  coordinate  represents  the  best  statistical  estimate  of  the  actual  vertical 
coordinate  of  each  of  the  calibration  beads.  The  technique  of  least  squares  fitting  a 
quadratic  (parabola)  is  also  common  place  and  will  not  be  discussed  here. 

Here  it  should  be  noted  that  the  primary  reason  why  Kevlar®  is  used  for  the 
calibration  wire  is  because  its  weight  is  minimal  for  the  strength  characteristics  it 
possesses.  The  total  sag  experienced  over  the  20 1  meters,  even  with  the  beads  in  place, 
only  amounts  to  0.40  meters.  This  permits  relatively  large  errors  in  the  measured  down- 
range  coordinates  without  significantly  affecting  the  statistical  estimate  of  the  actual  z 
coordinate  for  each  calibration  bead. 

A  2.2.3  DOWN-RANGE  COORDINATES 

The  actual  x  ( down-range )  coordinate  of  each  calibration  bead  is  also  determined 
statistically  by  a  least  squares  fitting  technique  to  the  measured  x  coordinates.  Since  this 
fitting  technique  is  somewhat  unique  it  will  be  described  in  some  detail  herein.  If  we 
initially  assume  that  the  first  measured  x  location  is  precise  ( this  assumption  will  be 
relaxed  later)  and  use  the  known  arc  lengths  between  the  calibration  beads  the  calibration 
bead  locations  are  computed  using  the  differential  relationship  shown  below: 

ds^  =  dx^  +  dz^  A2. 1 

where  s  =  known  arc  length. 

Also  we  know  from  the  previous  section  that  z  is  a  quadratic  function  of  x,  or: 

z  =  Cl  x^  +  C2  X  +  C3  A2.2 

and  the  constants  C]  ,  C2  ,  and  C3  have  previously  been  determined  during  the  fitting 
process.  Then  Eq.  A2.1  can  be  rewritten  as  shown: 
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X 


A2.3 


s  =  I  ( a  +  b  X  +  c  dx 

X 

where,  a  =  4Ci  ,  b  =  4CiC2  ,  c  =  C2"^+l 
Integrating  Eq.  A2.3  we  obtain, 

s  =  [(2ax+b)X‘^^ / 4a]  [(4ac-b^) /  8a][a‘^^ In +  x a‘^  +  b a'^^/ 2)]  | A2.4 

X 

0 

note  that ,  X  =  ax^  +  bx  +  c 

Since  everything  in  Eq.  A2.4  is  known  except  the  down-range  coordinate  (x) 
associated  with  each  known  arc  length  ,  Eq.  A2.4  is  solved  iteratively  for  each  of  these 
corresponding  values.  These  computed  values  of  x  will  likely  be  in  error  since  the  first  x 
coordinate  was  arbitrarily  set  equal  to  the  measured  x  coordinate  of  the  first  calibration 
bead.  Also,  the  calibration  wire  may  have  stretched  or  contracted  since  the  beads  were 
positioned  on  the  wire.  However,  the  down-range  distances  between  these  calibration 
beads  are  now  known  as  accurately  as  the  arc  lengths  were  known  when  the  wire  was 
constructed. 

In  order  to  understand  the  method  of  determining  the  actual  x  coordinate  of  each 
calibration  bead,  visualize  a  scale  containing  the  computed  x  coordinates  positioned 
alongside  the  measured  x  coordinates  of  each  calibration  bead.  Slide  this  scale  containing 
the  computed  x  coordinates  back  and  forth  and  permit  the  scale  to  stretch  or  contract  until 
the  scale  best  matches  the  measured  x  coordinates.  The  final  position  of  the  x  coordinates 
on  the  scale  now  represents  the  best  statistical  estimate  of  the  actual  down-range 
coordinate  at  each  of  the  measured  calibration  bead  locations. 

In  describing  this  mathematically  refer  to  the  sketch  below: 
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where, 


N  =  number  of  data  points 
Xin,  =  measured  x  coordinates 

x'ia  =  computed  x  coordinates  using  known  arc  length 
Xja  =  actual  X  coordinates  determined  statistically  from  least  squares  fit 
Sf  =  stretch  factor  between  the  x'jg  scale  and  the  Xj^  coordinates 
6x0  ~  shift  between  the  x'j^  and  Xj^  scales 


The  relation  relating  the  two  scales  can  be  written  as  follows: 

^ia  ~  ^  ia  ^xO  ^im 

If  we  define  a  delta  x  as, 


AXi  =  Xia-Xi, 


A2.5 

A2.6 


then  substituting  Eq.  A2.6  into  A2.5  and  rearranging  we  can  obtain  the  equation  below: 

AXi  =  x'ia  +  Sxo  +  Xim(  Sf  - 1 )  A2.7 

In  order  to  minimize  the  sum  of  AXj^ ,  the  following  relation  is  used, 

^/^6xoZ[x'ia  +  5xo  +  Xta(Sf-l)f  =0  A2.8 

i=l 

Equation  A2.8  can  be  simplified  to, 

6xo  =  -A-B(Sf- 1 )  A2.9 

where,  n  n 

A  =  (l/N)i:x'ia  ,  and  B=(l/N)Zxi„ 

i=I  i=l 

Now  if  we  substitute  Eq.  A2.9  into  A2.7  we  obtain, 

AXj  =  x'ia  -  A  +  B  -  Xjm  +  Sf  (Xin,  -  B)  A2.10 

2 

Again,  in  order  to  minimize  the  sum  of  AXi  , 

N 

^/<7SfZ[x'ia-A  +  B-Xin,  +  Sf(Xin,-B)f  =0  A2.11 

i=I 

Equation  A2.1 1  can  now  be  solved  for  the  stretch  factor  and  after  some  manipulation  the 
following  relation  obtained: 

Sf  =  [( NAB  -  E  Xi„  x'ia )  /  ( s  Xi„,'  -  N  B^ )]  +  1  A2.12 

i=l  i=l 

The  stretch  factor,  Sf ,  can  now  be  computed  directly  using  Eq.  A2.12.  Then  the  zero  shift, 
6x0,  is  computed  using  Eq.  A2.9,  and  finally  the  best  statistical  estimate  for  each  of  the 
down-range  coordinates  computed  using  Eq.  A2.5. 


This  completes  the  methodology  behind  the  determination  of  the  x,y,z,  coordinates 
of  the  calibration  beads.  We  can  now  assume  that  the  location  of  these  beads  are  precisely 
known.  Obviously  much  reliance  has  been  placed  on  the  fact  that  a  least  squares  fit  to 
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measured  data  is  a  better  representation  of  reality  than  any  one  measured  data  point.  This 
is  well  known  and  a  proven  fact  as  long  as  the  mathematical  model  used  to  describe  the 
measured  data  points  adequately  represents  the  true  physics  of  the  system. 

A2.3  SPARK  GAP  AND  REFERENCE  BEAD  LOCATIONS 

The  logic  behind  the  computation  of  the  spark  gap  and  reference  bead  locations  will 
be  illustrated  in  this  section.  The  mathematics  are  similar  to  the  vector  intersection 
routine  as  discussed  in  Section  5.2  with  two  exceptions.  The  first  is  that  here  we  are 
looking  for  the  intersection  of  three  vectors  instead  of  two  and  secondly  we  will  vary  one 
set  of  coordinates  (shadow  coordinates  of  the  known  beads  in  space).  The  mathematical 
details  are  presented  in  Ref.  3.13  and  not  repeated  here,  but  a  simplified  discussion  using 
a  planar  system  will  be  used  to  illustrate  the  process.  This  simplified  planar  system  is 
shown  in  the  sketch  below. 


Where,  X] ,  X2  ,  X3  =  known  points  in  space 

xis ,  X2s  ,  X35  =  shadow  locations  of  known  points 

^Ise  ’  ^2se  ’  ^3se  “  erroneous  shadow  locations 

Ax  =  error  in  reference  bead  location  ( same  error  in  shadow 

locations,  i.e.  Ax  =  Xjg  -  xjse ) 

As  illustrated  in  the  sketch  above,  if  the  location  of  the  reference  bead  is  in  error  by 
Ax  the  measured  shadow  locations  are  in  error  by  the  same  amount(  neglecting  film 
reading  error).  Also  illustrated  is  the  fact  that  the  vectors  defined  by  the  erroneous 
shadow  locations  and  the  known  points  in  space  do  not  intersect  at  a  point.  However,  if 
the  erroneous  reference  bead  location  is  varied  (each  shadow  location  is  varied  by  the 
same  amount)  until  the  vectors  do  intersect  at  a  point  then  the  final  location  of  the 
reference  bead  is  the  true  reference  bead  location  and  the  point  of  intersection  is  the 
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updated  spark  gap  location.  This  essentially  is  the  method  as  described  mathematically  in 
Ref.  3.13. 

A2A  RANGE  CALIBRATION  EVALUATION 

As  mentioned  previously  several  sets  of  shadowgrams  were  obtained  with  the 
calibration  wire  located  at  different  positions  in  the  range.  Three  sets  of  the  shadowgrams 
were  used  to  calibrate  the  spark  source  and  reference  bead  locations  using  the  method 
described  in  the  previous  section.  The  remaining  sets  of  shadowgrams  are  used  to 
establish  the  measurement  precision  and  to  evaluate  possible  station  biases.  The  beads  on 
these  shadowgrams  are  read  and  analyzed  as  if  they  represented  a  free-flight  test  item  in 
space.  The  measured  range  coordinates  of  calibration  beads  are  then  fit  using  the  straight 
line,  quadratic,  and  sliding-stretching  technique  as  previously  discussed.  The  differences 
between  the  measured  coordinates  and  the  statistically  determined  best  estimates  are  then 
analyzed  at  each  shadowgraph  station  to  determine  possible  measurement  biases.  These 
biases  are  then  recorded  and  subtracted  from  all  future  measurements  up  imtil  the  time  the 
facility  is  recalibrated  and  the  process  is  repeated.  An  estimate  of  the  position 
measurement  precision  is  determined  from  the  variation  in  the  measurements  at  each 
station. 

The  precision  of  the  attitude  measurements  is  obtained  by  comparing  the  theoretical 
local  slopes  of  the  calibration  wire  at  each  shadowgraph  station  with  the  measured  slopes 
obtained  from  the  shadowgraph  film.  The  theoretical  local  slopes  can  easily  be 
determined  from  the  curve  fits.  The  measured  slopes  are  determined  by  locating  two 
points  on  the  wire  at  each  shadowgraph  station  and  computing  the  corresponding 
attitudes. 

Using  the  techniques  discussed  above  the  measuring  capability  of  the  ART  is 
approximately  0.08  degrees  and  0.08  cm  respectively  for  the  orientation  and  position  of 
well  defined  points  in  space.  This  measurement  precision  has  also  been  confirmed  from 
the  analysis  of  numerous  free-flight  trajectories  as  measured  in  the  facility.  However,  it 
must  be  noted  that  the  measurement  precision  can  be  a  function  of  the  test  item’s  size  and 
geometry.  For  example  the  attitude  of  a  short  blunt  configuration  cannot  be  measured  as 
precisely  as  a  long  slender  configuration. 
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APPENDIX  3 


CURVE  FITTING  PROCEDURE  FOR  THE 
TRICYCLIC  EQUATION  OF  MOTION 


A  3.1  STATEMENT  OF  THE  PROBLEM 

The  problem  at  hand  is  to  fit  the  theoretical  equation  of  the  angular  motion  to  the 
experimentally  measured  data  points  and  determine  the  unknowns  such  that  they 
represent  a  best  fit  to  the  data.  The  theoretical  equation  of  motion  (Eq.  6.20)  is  repeated 
below, 

p  +  ia  =  Kp  exp[(^p  +  icjipOx]  +  Kn  exp[(pN  +  i<t>N')x]  +  K-p  exp(ipx)  A3 . 1 

and  the  unknowns  to  be  determined  are:  Kp,  pp,  (j)p',  K^,  Pn,  (j)N\  ^nd  K-p.  As  was  indicated 
in  Section  6.1.3  the  spin  rate,  p,  is  a  function  of  (})p'and  ([in'  and  therefore  is  not  considered 
as  an  unknown  in  Eq.  A3. 1 . 

A  3.2  LEAST  SQUARES  THEORY 

The  residual  is  defined  as  the  difference  between  the  fitted  angles,  Pp  +  iap,  and  the 
measured  angles.  Pm  +  iaM,  as  indicated  in  the  sketch  below: 


Fitted  Curve 


For  convenience  we  will  redefine  the  angles  as  §p  =  Pp  +  iap  and  =  Pm+  i(XM  then  the 
residual,  R^,  can  be  written  as  shown  below: 


Where  f(  x^)  is  the  tricyclic  equation  of  motion,  Eq.  A3.1  written  for  the  measured 
distance  x^,  or 

f(  Xm)  =  Kp  exp  (j)p  Xm  +  Kn  exp  x^  +  Kj  exp  Xm  A3.3 

where 

Kn  =  a  +  ib  ,  Kp  =  c  +  id  ,  Kj  =  f  +  ig 

<I>N  =  Fn  +  »  <^p  =  Fp  +  i<t>p'  .  <t'T=  ip 

and 

p  =  (V  I,)(  (^n'  +  <t)p') 
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The  values  of  K^,  Kp,  Kj,  (|)n,  and  (j)p  which  satisfy  the  least  squares  requirement  that 
the  sum  of  the  residuals  squared  be  a  minimum  can  be  found  by  a  simultaneous  solution 
of  the  set  of  equations  shown  below: 


^ERm  /  =0 

^IRm' 


7  ^  (|)m  =  0 


^  E  Rm  /  ^  Kp  =  0 
E  Rm^  /  ^  ct>p  =  0 


^ZRm"/  ^Kt  =0 


A3 .4 


This  set  of  equations  (ASA)  is  nonlinear  in  respect  to  the  unknowns  and  hence  cannot  be 
solved  directly.  Therefore,  an  interative  technique  must  be  used  for  the  fitting  process. 

A  3.3  THE  DIFFERENTIAL  CORRECTIONS  PROCEDURE 

The  differential  corrections  procedure  is  an  iterative  method  which  can  be  used  to  fit 
Eqs.  A3. 4  to  the  required  precision.  We  are  given  a  set  of  data  points  and  the  function  Eq. 
3.3  written  as. 


^  =  Kp  exp  (jip  X  +  K]vj  exp  x  +  K-p  exp  (j)  j  x 
Note  that  ^  =  p  +  ia  and  the  unknowns  (Kp,  (jip,  K^,  <1)^,  and  Kj)  are  defined  as, 
Kn  =  a  +  ib  =  (ao  +  Aa)  +  i(  bo  +  Ab)  =  K^o  +  AK^ 

Kp  =  c  +  id  =  (cq  +  Ac)  +  i(  do  +  Ad)  =  Kpo  +  AKp 
Kp  s  f  +  ig  =  (fo  +  Af)  +  i(  go  +  Ag)  =  Kpo  +  AKp 

<I>N  =  I^N  +  ^  M-NO  +  APn  i(‘t>No'  A(j)N')  =  <t>N0  A(t)N 

(jip  s  Pp  +  i(j)p'  =  Ppo  +  App  +  i((t)po'  +  A(t)p')  =  (jjpo  +  A(j)p 


Then  it  becomes  obvious  that, 
Kno  ~  ^  ibo 

Kpo  ~  Co  +  ido 
Kpo  ~  fo  lEo 

<t>N0  “  M’NO 
<t*P0  ~  M'PO  i^t’po’ 


AKn  =  Aa  +  iAb 
AKp  =  Ac  +  iAd 
AKp  =  Af  +  iAg 
A<t>N  ~  Apn  +  iA(|)f,j' 
A(|)p  =  App  +  iA(j)p' 


>  A3.5 


The  subscripted  zero  values  represent  the  first  approximation  to  the  unknowns  and  the 
delta  terms  are  the  corrections  added  after  each  iteration.  We  now  need  to  determine  these 
corrections  such  that  the  sum  of  the  residuals  squared  is  a  minimum.  Substituting  the 
quantities  above  into  the  residual  equation,  Eq.A3.2,  that  equation  can  be  written  as 
shown: 
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Rm  +  (Kpo  +AKp)  exp[((|)po  +  A(|)p)]  Xm+  (K^q  +AKN)exp[((i)No  +  A<|)n)]  x^ 

(^TO  +AK-r)  exp[(<j)7o  ■*"  A<|)i-)]  Xj^ 

Expanding  the  above  equation  with  Taylor’s  theorem  we  can  obtain, 

Rm  +  =  Kpo  exp  (t)po  Xm  +  Kno  exp  (t)No  Xm  +  Kto  exp  (|)to  Xm  +  AKp  ^  F  /  <^  Kpo 

+  AKn^  F  Id  Kno  +  AKt^  F  Kto  +  A(t)p^  F  Id  (|)po 
+  A<j)N  dVId  <t)No  +  (  higher  order  terms)  A3 .6 


where. 


F  —  Kpo  exp  (|)po  Xm  +  Kno  exp  (j)No  Xm  ■*"  Kjq  exp  (jjjo  x^ 

When  the  first  approximations  are  reasonably  good  the  A’s  will  be  small  and  the  higher 
order  terms  are  neglected. 

It  is  now  necessary  to  define  another  residual  as  the  difference  between  the  ordinates 
of  the  first  approximation  curve  and  the  data  points,  or 

Ri  =  ^M-F 

Substituting  this  relationship  into  Eq.  A3. 6  we  arrive  at  the  equation  shown  below: 

Rm  =  AKp^  F  Id  Kpo+  AKn<?  F  Id  Kno  +  F  Kxo  +  Acjip^  F  (jipo 


+  A(j)N F  / (T*  (t»No  -  Ri  A3. 7 

Equation  A3. 7  can  be  written  for  each  data  point  and  this  series  of  equations  are  known 
as  the  residual  equations.  Using  this  definition  for  R^  the  condition  for  the  sum  of  the 
squares  to  be  minimum  is, 

^7  Z  Rm^  /  ^  AKn  =  0  ,  ^  Z  Rm^  /  ^  AKp  =  0  ,  ^  Z  Rm^  /  ^  AKx  =  0 
^7  Z  Rm^  /  ^  A(t)N  =  0  ,  ^  E  Rm^  ^  A(|)p  =  0 
Substituting  Eq  A3 .7  into  one  of  the  above  relations  we  get, 

^  Z  Rm^  /  ^  AKn  =  ^  Id  AKnZ  [AKp^7  F  Id  Kpo+  AKn^  F  Id  Kno  +  AKt^7  F  Id  Kjo 
+  A(t)p<7  ¥  I d  (jipo  +  A(t)N<7  F  Id  (|)no  -  Rj  ]^  =  0 
Taking  the  indicated  partial  derivative  this  equation  can  be  written  as, 

^7  Z  Rm^/ ^7  AKp  =  2  ZAi  [AKp  Aj  +  AKn  Bj  +  AKt  Cj  +  A^p  D;  +  E;  -  Rj]  =  0  A3.8 
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where, 


Ai=  ^  F Kpo ,  Bi=  F Kno  ,  Ci=  F /<?  Kto  , 
Di=  F  / (^  (t)po ,  Ej  =  c?  F/^  (j)No 


Similarly  the  other  partial  derivatives  can  be  obtained,  or 

^  I  Rm^/  AKn  =  2  iBj  [AKp  Ai  +  AKn  B;  +  AKj  Cj  +  A(j)p  D,  +  A(t)N  Ej  -  RJ  =  0  A3 .9 
E  AKt  =  2  ECi  [AKp  Aj  +  AKn  Bj  +  AKt  Cj  +  A(l)p  D-,  +  A(j)N  E;  -  Rj]  =  0  A3. 10 
E  Rm^/^  A<l)p  =  2  EDj  [AKp  Ai  +  AKn Bi  +  AK-f  C;  +  A(t)p  Dj  +  A(t)N  Ej  -  Rj]  =  0  A3. 1 1 
E  Rm^/^  A(t)N  =  2  EEi  [AKp  Ai  +  AKn  Bi  +  AKj  C,  +  A(t)p  D;  +  A(t)N  Ei  -  Ri]  =  0  A3. 12 


The  summations  are  accomplished  over  the  complete  data  set  or  1  to  N.  These  five 
equations  (A3. 8  to  A3. 12)  are  usually  written  in  matrix  form  and  solved  for  the  five 
unknowns:  AKp,  AKn,  AKy,  A(j)p,  and  A(|)n.  Once  the  A’s  are  determined  they  are  added  to 
the  initial  approximations  and  the  process  repeated.  This  iterative  process  is  continued 
until  the  A’s  approach  zero  and  then  convergence  is  achieved. 
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